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SECTION  I 


INTRODUCTION  AND  REPORT  ARRANGEMENT 


This  report,  in  two  volumes,  describes  a computer  program  that 
has  been  developed  for  the  design  of  axial  compressors.  It  is  an  updated 
version  of  the  program  presented  in  References  1 and  2.  The  principal 
purpose  of  the  program  is  to  enable  a single  computer  program  to  deter- 
mine the  geometry  of  the  compressor  blading,  details  of  the  flow  within 
the  compressor,  and  the  design  point  performance  of  the  machine.  Some 
optional  calculation  routines  will  also  enable  effects  of  mixing  of  the  flow 
to  be  investigated.  The  program  consists  fundamentally  of  three  sections: 
two  alternative  means  of  determining  blade  geometry,  and  an  aerodynamic 
computation  for  the  flow  through  the  compressor. 

One  n ethod  of  determining  blade  geometry  involves  the  use  of  analytic 
meanlincs  for  the  blade  sections.  This  section  of  the  program  is  essentially 
the  computer  program  described  in  Reference  3,  modified  according  to 
Reference  5.  These  modifications  are  the  difference  between  this  section 
of  the  new  program  and  the  same  section  of  the  original  program  of  Refer- 
ences 1 and  2,  and  extend  the  capabilities  of  the  program  section  to  include 
the  optional  generation  of  "splitter"  blades  instead  of  conventional  "princi- 
pal" blades,  and  also  the  optional  output  of  data  compatible  with  the  NASTRAN 
stress  analysis  program.  Output  from  the  section  optionally  includes 
relative  flow  angles  through  the  blades.  These  may  be  used  as  partial  input 
data  to  the  aerodynamic  determination  section  of  the  program  which  will  then 
yield  the  flow  characteristics  for  the  blading.  An  alternative  input  form  to 
the  aerodynamic  section  is  the  specification  of  the  desired  angular  momentum 
distribution  through  the  blade.  Output  from  this  calculation  -v ill  include 
relative  flow  angles  through  the  blade.  These  may  be  used  as  input  to  the 
third  section  of  the  program,  the  generation  of  blades  having  arbitrary 
section  meanlines.  This  section  is  based  upon  the  computer  program 
described  in  Reference  4,  modified  according  to  Reference  5.  The  modifi- 
cations described  in  Reference  5 extend  the  program  section  capabilities 
as  for  the  analytic  meanline  section  mentioned  above,  and  further,  the 
method  of  generating  the  meanline  between  the  points  prescribed  by  the 
aerodynamic  section  calculation  has  been  improved.  For  both  blading 
methods,  flow  deviation  angles  link  the  blade  angles  to  relative  flow  angles. 

The  aerodynamic  section  was  written  anew  for  the  original  program 
of  References  1 and  2,  and  has  been  modified  somewhat  in  this  revised 
version.  The  principal  changes  are  the  addition  of  an  optional  alternative 
form  of  momentum  equation  for  the  meridional  velocity  profiles  and 
optional  output  of  pressure  force  data  in  a format  compatible  with  the 
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NASTRAN  stress  analysis  program.  The  aerodynamic  section  is  described 
in  detail  in  this  report,  and  sufficient  information  is  given  regarding  the 
blading  sections  to  enable  them  and  the  whole  program  to  be  used,  but  for 
finer  details  of  the  blading  calculations,  the  references  given  above  should 
be  consulted. 

The  following  three  sections  of  this  report  give  details  of  the  capab.li- 
ties  and  method  of  the  program,  the  input  data  requirements,  and  the 
resulting  output  data.  These  sections  should  give  the  casual  program  user 
all  information  necessary  to  use  the  program.  Section  V details  the  theory 
of  the  methods  used  in  the  aerodynamic  section.  This  will  enable  the 
aerodynamic  flow  model  to  be  evaluated,  and  the  various  program,  limitations 
ar.d  capabilities  to  be  more  fully  understood.  The  system  of  equations 
derived  is  necessarily  solved  numerically  and  iteratively,  and  the  various 
numerical  procedures  involved  are  described  in  Section  VI.  Details  of 
the  FORTRAN  programming,  including  implementation  of  the  program  on 
a particular  computing  system,  are  given  in  Section  VII. 

FORTRAN  listings  of  the  program  and  a sample  run  are  shown  in 
Volume  II  of  the  report. 
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SECTION  II 


SCOPE  AND  GENERAL  METHOD  OF  PROGRAM 

1.  OVERALL  PROGRAM 

In  one  run  of  the  program,  any  one  section  alone  may  be  used  or  the 
aerodynamic  section  may  be  used  in  conjunction  with  either  or  both  of  the 
blading  sections.  When  the  blading  sections  are  used  alone,  it  is  assumed 
that  one  blade  only  will  be  produced  in  each  run.  When  they  are  used  in 
conjunction  with  the  aerodynamic  section,  any  number  of  blades  may  be 
produced  by  either  method,  subject  only  to  the  limitations  imposed  by  the 
number  of  computing  stations  the  aerodynamic  section  can  handle.  When 
the  analytic  meanline  blade  section  is  used  ir.  conjunction  with  the  aero- 
dynamic section,  the  data  from  the  blade  calculations  that  are  carried  over 
to  the  aerodynamic  analysis  consists  essentially  of  the  blockage  caused  by 
the  blade  and  the  blade  angles.  Two  angles  are  involved,  the  "section 
angle"  defined  by  the  intersection  of  the  blade  mean  surface  with  a cylinder, 
and  the  "lean  angle"  defined  by  the  intersection  of  the  blade  mean  surface 
with  a computing  station.  Additional  data  describing  the  blade  perform- 
ance are  required  before  an  aerodynamic  analysis  can  be  performed, 
namely,  the  deviation  angles  and  a specification  of  the  losses.  This  is 
combined  with  the  data  from  the  blading  section  in  an  interface  routine  and 
then  the  completed  data  are  passed  to  the  aerodynamic  section.  When  the 
arbitrary  meanline  blade  section  is  used  following  an  aerodynamic  design 
calculation,  the  data  carried  over  to  the  blading  calculation  consist 
essentially  of  the  blade  section  angles  for  the  blade.  The  aerodynamic 
design  calculation  produces  the  relative  flow  angles,  which,  in  a manner 
similar  to  that  mentioned  above,  are  combined  with  the  flow  deviation 
angles  except  that  in  this  case  the  calculation  was  included  in  the  blade 
program  described  in  Reference  4. 

The  program  is  written  in  standard  FORTRAN  IV  and  should  be 
immediately  compatible  with  aL'  mediun  „o  large  computing  systems.  All 
development  running  was  done  on  a CDC  6600  system  incorporating 
GALCOMP  software  and  on-line  hardware  at  Wright- Patterson  Air  Force 
Base.  The  program  uses  overlay;  each  of  the  main  sections  uses  close  to 
130K  (octal)  of  central  memory. 

2.  AERODYNAMIC  SECTION 

The  compressor  flow  computed  in  the  aerodynamic  section  is 
assumed  to  be  axisymmetric  and  inviscid.  The  two  alternative  momentum 
equations  include  entropy  gradients  in  both  the  cross  streamwise  and 
streamwisc  directions,  and  also  tht  blade  forces.  The  fluid  properties  are 
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all  computed  in  a group  of  FUNCTION  SUBPROGRAMS.  As  shown  in  this 
report,  the  properties  are  computed  for  a perfect  gas.  This  may  be 
changed  by  replacing  these  routines.  The  streamline  curvature  method  of 
solution  is  employed  to  solve  the  system  of  equations.  In  this  method, 
a number  of  computing  stations  are  located  at  strategic  points  in  the  flow 
and  preferably,  but  not  necessarily,  close  to  orthogonal  to  the  local 
meridional  streamline  direction.  Typically,  several  stations  are  situated 
upstream  of  the  first  blade  row;  one  will  be  at  the  edge  of  each  blade  row 
and  several  will  be  downstream  of  the  last  blade  row.  Additionally, 
stations  may  be  placed  within  the  blade  rows,  and  this  is  normally  done 
when  detailed  blade  evaluations  are  being  performed  for  high  speed  blades. 
The  stations  may  be  defined  quite  freely;  they  need  not  be  radial  nor  need 
they  be  straight  lines.  A computing  mesh  is  formed  by  the  intersection  of 
these  (fixed)  stations  with  the  streamlines,  whose  locations  are  iteratively 
determined.  An  initial  estimate  is  made  of  the  streamline  locations.  The 
flow  within  the  compressor  is  computed  on  the  basis  of  this  estimate  and 
the  resulting  flow  distribution  enables  a new  estimate  to  be  made.  This 
procedure  is  repeated  until  the  estimated  streamline  pattern  is  correct 
(to  within  a given  tolerance).  Up  to  30  computing  stations  and  21  stream •* 
lines  may  be  used  to  describe  the  flow  field. 

Losses  occurring  for  each  blade  row  may  be  specified  by  means  of 
the  radial  distribution  of  relative  total  pressure  loss  coefficient,  isen- 
tropic  efficiencies,  or  entropy  increases.  Methods  of  specifying  the 
angular  momentum  fall  into  two  classes,  "design"  and  "analysis".  For 
the  design  cases,  the  radial  distribution  of  total  pressure,  total  enthalpy, 
absolute  angular  momentum,  or  absolute  whirl  velocity  may  be  given. 

For  the  analysis  cases,  the  relative  flow  angle  distribution  is  required  and 
this  may  be  specified  directly,  or  a blade  angle  and  deviation  angle  may  be 
given.  Relative  flow  angle  or  blade  angle  may  be  given  in  the  streamline 
direction,  or  as  would  occur  were  the  stream  surfaces  all  concentric 
cylinders.  An  off-design  analysis  option  involving  relative  total  pressure 
loss  coefficients  and  relative  flow  angles  is  included.  In  this,  a mimimum 
loss  coefficient  and  a reference  relative  outlet  flow  angle  are  given  and 
deviations  from  these  values  computed  as  a function  of  relative  inlet  angle 
and  its  reference  value.  In  order  to  produce  designs  in  which  the  blade 
row  losses  are  consistent  with  a loss  model  (rather  than  being  specified  in 
the  data  and  then  updated  manually  between  computer  runs),  a loss  re- 
estimation procedure  is  incorporated  which  may  be  used  to  continuously 
update  relative  total  pressure  loss  coefficient  values  as  the  computation 
proceeds.  A second  option  using  the  same  routine  is  to  only  use  it  to 
compute  and  display  new  loss  coefficients  after  the  main  solution  is 
computed.  Although  intended  as  a design  calculation  aid,  this  routine  may 
be  employed  regardless  of  whether  a design  or  analysis  method  of 
determining  angular  momentum  is  employed  for  the  blade  row. 
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Optionally,  the  blockage  due  to  boundary  layers  on  the  annulus  walls 
may  be  computed  from  a simple,  attached  turbulent  boundary  layer 
equation. 

Provision  is  made  in  the  procedure  to  incorporate  two  viscous  flow 
effects,  which  are  the  turbulent  mixing  of  the  flow  and  the  radial  transfer 
of  the  blade  wakes.  Used  in  conjunction  with  the  off-design  analysis 
option  mentioned  above,  this  should  enable  the  effect  of  these  phenomena 
on  multistage  compressor  performance  to  be  investigated. 

When  the  analytic  meanline  blade  section  is  used  to  create  blade 
geometry  prior  to  execution  of  the  aerodynamic  section,  an  interface 
routine  determines  the  flow  deviation  angles.  Then  the  relative  flow  angles 
required  by  the  aerodynamic  section  are  the  sum  of  the  blade  angles  and 
the  deviation  angles.  The  deviation  angles  are  calculated  basically 
according  to  the  method  proposed  in  Reference  6,  Equations  269  and  271 
In  this  method,  the  deviation  angles  are  composed  of  a reference  deviation 
for  a 1 9%  thick  urcambere-1  section,  modified  for  blade  type  and  actual 
thickness,  plus  a component  depending  upon  camber,  inlet  angle,  solidity, 
and  meanline  form.  Two  extensions  to  the  method  have  been  made.  First, 
the  second  component  has  been  irade  a function  also  of  point  of  maximum, 
camber,  and  second,  provisions  have  been  made  to  add  in  an  arbitrary 
additional  component  of  deviation  angle. 

Optionally,  output  from  the  aerodynamic  section  is  data  describing 
the  aerodynamic  loading  on  the  blade(s)  in  a format  compatible  with  the 
input  requirements  of  the  NASTRAN  stress  analysis  program.  A convention 
has  been  established  for  locating  the  elements  into  which  the  blade  is 
divided  that  has  been  followed  in  each  of  the  three  program  sections  where 
NASTRAN-compatible  output  is  concerned.  Reference  7 is  a users  manual 
for  the  NASTRAN  program. 

3.  ANALYTIC  MEANLINE  BLADE  SECTION 

As  mentioned  previously,  this  section  is  basically  the  program 
described  in  Reference  3,  modified  according  to  Reference  5.  The  blade 
is  defined  by  first  producing  a series  of  blade  sections,  one  on  each  of  a 
number  of  stream  surfaces  (surfaces  of  revolution)  through  the  blade,  and 
then  stacking  them  in  a prescribed  manner  to  produce  the  blade.  "Manu- 
facturing section"  definitions  in  Cartesian  coordinates  may  then  be 
produced  by  interpolation  at  a series  of  planes  perpendicular  to  the  stacking 
axis. 


Four  alternative  meaniines  are  available  for  the  blade  section.  (Any 
one  blade  uses  one  meanline  type  throughout.)  The  first  one  consists  of  a 
fourth  order  polynomial  in  which  the  coefficients  are  determined  by 


satisfying  the  desired  inlet  and  outlet  angles  and  also  specifying  the  ratios 
of  second  derivative  at  the  inlet  and  outlet  to  the  maximum  value  occurring 
on  the  section.  This  is  particularly  suited  to  blades  having  inlet  Mach 
numbers  in  the  range  of  approximately  0.  8 to  1.7,  where  a relatively 
straight  leading  edge  followed  by  increasing  camber  is  desirable.  The 
second  meanline  is  composed  of  two  exponential  functions  and  allows 
additional  degrees  of  freedom.  Similar  specifications  are  made  at  the  inlet 
and  outlet,  and  additionally,  the  location  of  the  point  where  the  two  mean- 
line equations  meet  and  the  blade  angle  at  that  point  are  specified.  This 
enables  "S"  sections  to  be  specified,  in  addition  to  sections  having  contin- 
uously positive  camber  or  a straight  leading  edge  region  and  a cambered 
aft  portion.  Thus,  a blade  may  be  produced  that  is  potentially  useful  in 
the  Mach  number  range  of  approximately  0.  8 to  over  2.  0.  The  third 
meanline  is  simply  a circular  arc  so  that  only  the  inlet  and  outlet  angles 
may  be  defined.  Blades  based  upon  this  meanline  have  been  successfully 
used  at  subsonic  and  low  transonic  Mach  numbers.  The  final  meanline  is 
composed  of  two  circular  arcs.  Both  the  inlet  and  outlet  angles  and  the 
location  of  the  point  where  the  two  meanline  segments  meet  and  the  blade 
angle  at  this  point  are  required  to  be  specified.  This  meanline  has 
similar  capabilities  to  the  exponential  meanline  described  above,  and  is 
currently  in  use  in  industry.  In  order  to  complete  the  blade  section 
specification  (with  any  of  the  above  mentioned  meanlines),  a thickness 
definition  is  required.  Two  thickness  specifications  are  used.  The  first 
consists  of  two  third  order  polynomials,  one  applying  from  the  leading 
edge  to  the  point  of  maximum  thickness,  and  the  second  to  the  remaining 
portion  of  the  section.  By  varying  the  location  of  the  point  of  maximum 
thickness,  and  the  thickness  at  the  section  edges  and  at  the  point  of 
maximum  thickness,  thickness  distributions  compatible  with  mechanical 
requirements  and  all  but  the  lowest  speed  aerodynamic  requirements  may 
easily  be  produced.  This  thickness  is  applied  to  all  meanlines  except  the 
(single)  circular  arc.  However,  this  limitation  can  be  overcome  by  using 
the  multiple  circu'ar  arc  meanline  with  the  -extent  of  one  arc  set  to  some 
insignificant  fraction  of  the  total  chord.  For  the  case  of  the  circular  arc 
meanline  only,  the  thickness  distribution  is  set  so  that  the  two  surfaces 
oroduced  are  also  circular  arcs,  producing  the  so-called  double  circular 
arc  section. 

The  stacking  ol  the  blade  ?s  achieved  by  locating  each  of  the  (two- 
dimensional)  sect. ons  on  the  appropriate  surface  of  revolution  to  form  a 
three-dimensional  blade  definition.  There  are  available  three  variations 
in  the  stacking  method.  The  first  is  to  stack  the  section  (that  is,  pass  the 
stacking  axis,  a radial  line)  through  the  two-dimensional  section  centroids, 
or  points  specified  relative  to  the  centroid.  This  will  normally  be  used  for 
rotor  blades  for  reasons  of  stress.  The  other  two  methods  are  to  stack 
the  blade  at  either  the  leading  or  trailing  edge. 
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The  interpolation  of  sections  in  Cartesian  coordinates,  and  there 
fore  convenient  for  the  manufacturing  process,  is  achieved  by  fitting  a 
spline  curve  through  points  on  the  blade  surfaces  and  noting  where  it 
passes  through  the  desired  manufacturing  planes. 


Up  to  80  points  may  be  used  to  define  each  blade  surface,  and  blade- 
sections  on  up  to  21  stream  surfaces  may  be  computed.  The  stream  surfaces 
may  be  defined  at  up  to  ten  stations  distributed  upstream,  within,  and  down- 
stream of  the  blade.  When  the  blade  program -section  is  used  in  conjunction 
with  the  aerodynamic  program -section,  the  streamsurfaces  and  stations 
are  common  to  the  two  calculations,  and  hence,  iteration  between  the  two 
program-sections  is  required  to  achieve  similarity  between  the  stream  - 
surfaces  that  are  specified  to  the  blading  pregram- section  and  those 
generated  in  the  ensuing  aerodynamic  analysis. 

A further  capability  of  the  analytic  meanline  program  section  is  the 
option  to  design  "splitter"  blades.  Splitter  blades  are  defined  as  blades 
having  the  same  trailing  edge  location  at  the  principal  blades,  a leading 
edge  location  somewhere  between  the  leading  and  trailing  edges  of  the 
principal  blades,  and  the  same  meanline  as  that  portion  of  the  principal 
blades  that  they  duplicate.  (Typically,  the  chord  of  splitter  blades  will  be 
about  one  half  that  of  the  principal  blades.)  Two  thickness  distributions 
may  be  used  in  conjunction  with  the  meanline  obtained  in  this  way;  the 
double -circular -arc  (DCA)  thicKiess  distribution,  and  the  ’ standard" 
double -polynom ial  thickness  distribution  described  above  for  the  princ.oal 
blades. 

Also,  output  data  may  be  produced  that  can  be  used  directly  as  input 
to  the  NASTRAN  stress  analysis  program.  Cards  are  created  that  describe 
the  grid  and  elements  that  the  blade  has  been  divided  into.  These  are 
geometric  properties  of  the  blade.  Cards  may  also  be  produced  that 
reflect  the  aerodynamic  loading  on  the  elements.  As  this  program  section 
is  basically  concerned  with  blade  geometric  properties,  extra  input  is 
required  to  generate  these  loading  figures.  These  may  be  derived  from  the 
output  of  a previous  run  of  the  aerodynamic  section.  Alternatively,  the 
aerodynamic  section  itself  may  be  directed  to  output  the  pressure-load 
data  in  NASTRAN-compatible  form. 

4.  ARBITRARY  MEANLINE  BLADE  SECTION 

As  mentioned  previously,  this  program  section  is  based  upon  the 
program  described  in  Reference  4 and  modified  according  to  Reference  5. 

As  for  the  Analytic  Meanline  case,  the  blade  is  defined  by  first  producing 
a series  of  blade  sections,  one  on  each  of  a number  of  streamsurfaces. 


7 


and  then  stacking  them  to  produce  the  blade.  "Manufacturing  section" 
definitions  in  Cartesian  coordinates  may  then  be  produced  by  interpolation 
at  a series  of  planes  perpendicular  to  the  stacking  axis. 

The  unique  feature  of  this  blading  method  is  the  procedure  used  to 
produce  the  section  meanlines.  For  each  section,  input  to  the  calculation 
includes  relative  flow  angles  at  a number  of  points  along  the  meridional 
chordline.  (These  may  either  be  input  directly  or  have  been  produced  by 
the  aerodynamic  section.)  The  incidence  a.ngle  at  the  leading  edge  is 
specified  and  hence  the  blade  section  inlet  angle  is  known.  The  trailing 
edge  deviation  angle  is  determined  from  an  estimate  of  the  cascade  solidity, 
because  the  actual  chord  is,  as  yet,  unknown.  Then,  a specification  of 
the  variation  of  deviation  angle  within  the  blade  from  the  incidence  angle  at 
the  leading  edge  to  the  trailing  edge  value  yields  the  blade  angle  at  each 
point  where  the  relative  flow  angle  was  given.  A spline-fit  is  made  of 
the  tangent  of  this  angle  distribution  versus  axially -projected  chord  length, 
and  the  resulting  piece-wise  cubic  analytically  integrated  to  yield  a piece- 
wise  quartic  meanline.  This  is  a departure  from  the  method  presented  in 
the  original  program  of  References  1 and  2.  In  that  case,  a piece-wise 
cubic  meanline  was  produced,  but  it  was  necessary  to  compute  meanlines 
for  a range  of  end-point  conditions  of  the  first  cubic  segment  in  order  to 
find  a satisfactory  solution,  and,  in  some  cases,  no  feasible  meanline  was 
ever  found.  The  new  approach  appears  to  produce  an  attractive  meanline 
in  all  cases  with  no  difficulty. 

The  length  of  the  resulting  meanline  enables  the  cascade  solidity 
estimate  to  be  revised,  and  some  iteration  is  normally  required  to  get  the 
trailing  edge  deviation  angle  consistent.  This  is  performed  automatically 
by  the  program. 

The  remainder  of  the  procedure  to  produce  the  final  blade  definition 
is  similar  to  that  described  for  the  Analytic  Meanline  Blade  Section.  The 
two-part  cubic  thickness  distribution  is  applied  to  each  meanline  to  define 
the  section  surfaces.  Stacking  and  interpolation  of  "manufacturing  sections" 
proceed  as  described  above. 

Splitter  blades  may  be  designed  by  this  program  section  just  as  des- 
cribed for  the  Analytic  Meanline  Section,  except  that  in  this  case  only  the 
thickness  distribution  comprising  of  two  third-order  polynomials  is 
available. 

Output  for  subsequent  use  with  the  NASTRAN  stress  analysis  program 
may  be  produced  just  as  for  the  Analytic  Meanline  Section. 


INPUT  DATA 

1.  DATA  FORMAT 

The  input  data  consist  of  an  initial  indication  of  the  number  of 
entries  that  are  to  be  made  to  each  of  three  program  sections,  and  then  a 


data-set  for  each  entry  to  each  section.  The  data  that  are  required  for  the 
interfacing  of  the  output  from  the  analytic  meanline  blade  section  to  the 
aerodynamic  section  are  included  in  the  data-set  for  the  analytic  meanline 
section.  Because  input  to  the  aerodynamic  and  arbitrary  meanline  sections 
may  be  generated  by  previous  execution  of  the  analytic  meanline  and  aero- 
dynamic sections,  respectively,  the  input  for  these  sections  that  are  to  be 
supplied  directly  by  the  user  varies.  This  is  indicated  in  the  chart  below 
by  giving  the  FORTRAN  variable  name  for  the  file  from  which  any  data 
are  taken  that  are  not  always  supplied  directly. 

LOG5  is  the  file  from  which  input  is  taken  that  may  be  generated  by 
the  analytic  meanline  section.  When  the  analytic  meanline  section  has  been 
directed  to  produce  data  for  the  aerodynamic  section  for  a particular  com- 
puting station,  LOG5  becomes  an  internally  generated  scratchfile.  Other- 
wise, LOG5  is  attached  to  the  standard  input  unit  and  the  user  supplies  the 
data.  LOG6  is  similarly  used  for  the  arbitrarv  meanline  section.  If  the 
aerodynamic  section  has  been  directed  to  produce  appropriate  data,  these 
will  be  stored  within  the  computer.  Otherwise,  the  user  supplies  the  data 
directly. 

Three  input  formats  are  used;  alphanumeric,  real,  and  integer.  Only 
one  type  of  data  occurs  on  any  one  card  (with  the  exception  noted  below), 
and  normal  FORTRAN  convention  indicates  which  quantities  are  integer. 

The  alphanumeric  format  is  used  for  the  four  title  cards.  Up  to  72  charac- 
ters may  be  used,  starting  in  column  1.  Real  numbers  are  punched  in 
fields  of  12  locations,  starting  in  cobimn  1 of  the  card,  up  to  six  numbers 
per  card.  Decimal  points  should  be  included  to  ensure  correct  interpre- 
tation of  the  data,  and  the  numbers  may  be  placed  anywhere  within  the 
allocated  fields.  Integers  are  punched  in  fields  of  three  locations,  starting 
in  column  1.  No  decimal  point  may  be  used,  and  the  numbers  must  occupy 
the  right-most  locations  of  the  allocated  fields.  The  one  exception  referred 
to  above  is  the  last  record  defined  for  the  aerodynamic  section.  It  consists 
of  three  real  numbers  input  in  the  usual  manner  in  the  first  36  columns, 
followed  by  two  integers  in  the  next  three  columns  each.  These  data  would 
normally  be  produced  by  a prior  run  of  the  program,  and  the  integers  are 
included  merely  to  check  that  the  cards  have  net  become  mis-ordered. 
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In  the  following  chart,  one  line  corresponds  to  one  card  (except 
where  more  than  one  line  is  required  to  complete  the  description  of  the 
card). 

TITLE1 

NANAL  NAERO  NARBIT 

The  following  data-set  is  input  to  the  analytic  meaniine  section, 
and  will  occur  NANAL  times.  The  last  record  in  this  set  is 
indicated  with  an  asterisk. 


(Con) 

(Con) 


TITLE2 


NLINES  NSTNS  NZ  NSPEC  NFOINT  N BLADE  ISTAK 
IPUNCH  ISECN  IFCORD  IFPLOT  IPRINT  ISPLIT  INAST 


IRLE  IRTE  NSIGN 


ZINNER  ZOUTER  SCALE  STACKX  PLTSZE 


KPTS  IFANGS 


XSTA  RSTA 


R B LA  FOR 


Occurs  KPTS  times 


Occurs  NLINES  times 


ZR  B1  B2  PP  QQ  RLE 


TC  TE  Z CORD  DELX  DELY 


S BS  - Only  if  ISECN  = 1 or  3 


Occurs 

NSTNS 

times 

Occurs 

NSPEC 

times 


RLES  TCS  TES  ZZS  PERSPT  - Occurs  NSPEC  times  if 

ISPLIT  ^ 1 


XI  X2  X3  XB  - Occurs  NLINES  times  for  each  station  where 
IFANGS  = 2 if  ISPLIT  :>  1 
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(Con) 

(Con) 


NRAD  NDPTS  NDATR  NSWITC  NLE  NTE 
XKSHPE  SPEED 

NOUT1  NOCJT2  NOUT3  -Refers  to  leading  edge  station. 
NR  NTERP  NMACH  NLOSS  NL1 
NL2  NEVAL  NCURVE  NLITER  NDEL 
NOUT1  NOUT2  N0UT3  NBLAD 
R XLOSS  -Occurs  NR  times 
RTE 

DM  DVFRAC  -Occurs  NDPTS  times 
RDTE  DELTAD  AC-Occurs  NDATR  times 


Occurs 

This 

for  each 

group 

station 

only 

within 

occurs 

blade  or 

if 

at  trailing 

NAERO 

edge 

= 1 or 
IPUNCH 

Occurs 

NRAD 

times 

= 1 

The  following  data-set  is  input  to  the  aerodynamic  section  and 
will  only  occur  if  NAERO  = 1.  The  last  record  in  this  set  is  indicated  with 
a double  asterisk. 

TITLE3 

CP  GASR  G EJ 

NSTNS  NSTRMS  NMAX  N FORCE  NBL  NCASE 
(Con)  NSPLIT  NSET1  NSET2  NREAD  NPUNCH  NPLOT 

(Con)  NPAGE  NTRANS  NMIX  NMANY  NSTPLT  NEQN 

NWHTCH  - Occurs  NMANY  times 
G EJ  SCLFAC  TOLNCE  VISK  SHAPE 
XSCALE  PSCALE  RLOW  PLOW  XMMAX  R CONST 
CONTR  CONMX 

FLOW  SPDFAC  - Occurs  NCASE  times 
NSPEC 

XSTN  RSTN  - Occurs  NSPEC  times 


Occurs 

NTNS 

tiu.es 


(LOG5) 

(Con) 

(Con) 

(LOG5) 

(LOG5) 

(Con) 

(LOG5) 


NDATA  NTERP  NDIMEN  NMACH 

DATAC  DATA1  DATA 2 DATA 3 -Occurs 

NDATA  times 

NDATA  NTERP  NDIMEN  NMACH  NWORK 
NLOSS  NL1  NL2  NEVAL  NCURVE  NLITER 
NDEL  NOUT1  NOUT2  NOUT3  NBLADE 
SPEED-If  NDATA  >0 


Inlet 

condition 

specification 


Occurs 

NDATA 

times 


DATAC  DATA1  DATA 2 DATA 3 DATA 4 
DATA  5 

DATA6  DATA 7 DATA8  DATA 9 
DELC  DELTA  - Occurs  NDEL  times 
W BLOCK  BBLOCK  BDIST  -Occurs  NSTNS  times 
NDIFF 

DIFF  FDHUB  FDMID  FDTIP  -Occurs  NDIFF 
NM  NRAD 
TERAD 

DM  WFRAC  -Occurs  NM  times 


For 

sta- 

tions 

2 

thru 

NSTNS 


Occurs 

NRAD 

times 


Occurs 

NSET1 

times 

Occurs 

NSET2 

times 


DELF(l)  DELF(2) DELF  (NSTRMS)  - if  NSPLIT  = 1 

or  NREAD  = 1 

R X XL  II  JJ  - Occurs  NSTRMS  times  for  NSTNS  stations 
if  NREAD  = 1 


The  following  (and  final)  data-set  is  input  to  the  arbitrary  mean- 
line section,  and  will  occur  NARBIT  times. 
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TITLE4 


NLINES  NSTNS  NZ  NSPEC  ISEGPT  NBLADE  ISTAK 
(Con)  IPUNCH  IFPLOT  IPRINT  ISPLIT  INAST 

ZINNER  ZOUTER  SCALE  STACKX  PLTSZE 
IRLE  IRTE  NR  A DEV  NINC  NSIGN  IFCA 
XKSHPE  SOL  TO  L 
NPTS 
RADEV 

SM  DEVCRV  -Occurs  NPTS  times 
RINC  XINC  DELDEV  - Occurs  NINC  times 
IFANGS 

(LOG6)  KPTS 

(LOG6)  XSTA  RSTA  - Occurs  KPTS  times 

(LOG6)  R AIRANG  BLAFOR  -Occurs  NUNES  I 

times  I 

ZR  YA  RLE  TC  TE  ZZ  Occurs 

NSPEC 

DELX  DELY  times 

RLES  TCS  TES  ZZS  PERSPT  -Occurs  NSPEC  times  if 

ISPLIT  s 1 

XI  X2  X3  XB  - Occurs  NLINES  times  for  each  station  where 
IFANGS  = 2 if  ISPLIT  5 1 


Occurs 

NSTNS 

times 


Occurs 

NRADEV 

times 


2.  DATA  ITEM  DEFINITIONS 

The  aerodynamic  section  may  be  used  with  any  self-consistent 
unit  system  and,  additionally,  a "linear  dimension  scaling  factor" 
(SCLFAC)  is  incorporated  into  the  input  so  that  some  commonly  used  but 
inconsistent  unit  systems  may  be  used.  This  is  principally  intended  to 
allow  the  use  of  inches  for  physical  dimensions  and  yet  retain  feet  for 
velocities.  The  basic  dimensions  used  in  the  data  are  length  (L),  time  (T), 
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and  force  (F).  Angles  are  expressed  in  degrees  (A),  and  temperatures  on 
an  absolute  temperature  scale  (D).  Heat  capacities  (H)  are  also  required. 
Some  possible  unit  systems  are  given  below,  together  with  the  corres- 
ponding value  of  SCLFAC. 


L 

T 

F 

D 

H 

SCLFAC 

Feet 

Seconds 

Pounds 

Deg.  Rankine 

BTU 

1.0 

Inches 

Seconds 

Pounds 

Deg.  Rankine 

BTU 

12.0 

Metres 

Seconds 

Kilograms 

Deg.  Kelvin 

CHU 

1.0 

Note  t!  ii  some  data  names  are  used  in  more  than  one  section;  care 
should  be  take;  to  consult  the  correct  sub-division  below  for  definitions. 


a.  Initial  Directives 


TITLE  1 


This  is  a title  card  for  the  run. 


NANAL  The  number  of  blades  to  be  generated  by  analytic  mean- 

line blade  section. 

NAERO  If  NAERG  = 0,  the  aerodynamic  section  will  not  be 

entered. 

If  NAERO  = 1,  the  aerodynamic  section  will  be  entered. 

NARBIT  The  number  of  blades  to  be  generated  by  the  arbitrary 

meanline  blade  section. 


b.  Analytic  Meanline  Blade  Section 

For  a more  detailed  discussion  of  the  input  to  this  section 
through  item  XB,  see  References  3 and  5.  For  this  section,  the  dimen- 
sioned input  is  either  in  degrees  (A)  or  in  length  (L). 


TITLE2 

NLINES 


NSTNS 


NZ 


A title  ca  -d  for  the  analytic  meanline  section  of  the 
program. 

Hie  number  of  stream  surfaces  which  are  defined,  and 
on  which  blade  sections  will  be  designed.  Must  satisfy 
2 £ NLINES  s 21. 

The  number  of  computing  stations  at  which  the  stream - 
surface  radii  are  specified.  Must  satisfy  3s  NSTNS  £ 10. 

The  number  of  constant-z  planes  on  which  manufacturing 
(Cartesian)  coordinates  for  the  blade  are  required.  Must 
satisfy  3 £ NZ  £ 1 5. 
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NSPEC 


The  number  of  radially  disposed  points  at  which  the 
parameters  of  the  blade  sections  are  specified.  Must 
satisfy  1 £ NSPEC  £21. 


NPOINT  The  number  of  points  that  will  be  generated  to  specify 

the  pressure  and  suction  surfaces  of  each  blade  section. 
Must  satisfy  2 £ NPOINT  £ 80.  Generally,  no  less  than 
30  should  be  used.  It  will  be  advantageous  to  specify  80 
points  when  precision  plots  of  the  sections  are  to  be 
produced  directly  by  the  program. 

NBLADE  The  number  of  blades  in  the  blade  row. 

ISTAK  If  ISTAK  = 0,  the  blade  will  be  stacked  at  the  leading 

edge. 

If  ISTAK  = 1 , the  blade  will  be  stacked  at  the  trailing 
edge. 

If  ISTAK  = 2,  the  blade  will  be  stacked  at,  or  offset 
from,  the  section  centroid. 

IPUNGH  If  the  data  for  subsequent  aerodynamic  analysis  are 

computed  (IFANGS  = 1),  they  will  also  be  punched  out  if 
IPUNCH  =1.  If  not  required,  set  IPUNCH  to  zero. 

ISECN  If  ISECN  = 0,  the  blade  will  be  constructed  using  the 

polynomial  camber  line  and  the  standard  (i.  e.  , double- 
cubic)  thickness  distribution. 

If  ISECN  = 1,  the  exponential  camber  line  and  the 
standard  thickness  distribution  will  be  used. 

If  ISECN  = 2,  the  circular  arc  camber  line  and  the  double  - 
circular-arc  thickness  distribution  will  be  used. 

If  ISECN  = 3,  the  multiple-circular-arc  meanline  and  the 
standard  thickness  distribution  will  be  used. 

IFCORD  If  IFCORD  = 0,  the  meridional  projections  of  the  stream  - 

surface  blade  section  chords  are  specified. 

If  IFCORD  = 1,  the  stream  surface  blade  section  chords 
are  specified. 
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IFPLOT 


IPRINT 


Where  CALCOMP  software  is  incorporated  into  the 
computing  system,  IFPLOT  specifies  the  creation  of 
precision  plots.  (Further  information  regarding  the 
requirements  for  this  are  given  in  the  section  entitled 
"Program  Operation  and  Structure.  ") 

If  IFPLOT  = 0,  no  plots  will  be  produced. 

If  IFPLOT  = 1,  a plot  of  the  streamsurface  sections  will 
be  produced.  All  NLINES  sections  are  shown  super- 
imposed. The  origin  for  each  section  plot  is  offset 
from  the  centroid  of  the  section  by  distances  specified 
by  DELX  and  DELY. 

If  IFPLOT  = 2,  a plot  of  th*  manufacturing  sections  will 
be  produced.  The  origin  is  the  blade  stacking  axis,  and 
all  NZ  sections  are  shown  superimposed. 

If  IFPLOT  = 3,  both  of  the  plots  described  for  IFPLOT 
= 1 and  2 will  be  produced. 

If  IFPLOT  = 4,  individual  plots  of  each  of  the  manu- 
facturing sections  will  be  produced.  The  axes  are 
rotated  clockwise  by  the  section  stagger  angle  for  each 
plot,  so  that  a large  scale  may  be  used. 

The  input  data  is  always  listed  by  the  program.  Details 
of  the  streamsurface  and  manufacturing  sections  are 
printed  as  prescribed  by  IPRINT. 

If  IPRINT  = 0,  details  of  streamsurface  and  manufacturing 
sections  are  printed. 

If  IPRINT  = 1,  details  of  streamsurface  sections  -e 
printed. 

If  IPRINT  = 2,  details  of  manufacturing  sections  are 
printed. 

If  IPRINT  = 3,  details  of  neither  streamsurface  nor 
manufacturing  sections  are  printed.  (The  interface  data 
for  use  with  the  aerodynamic  section  of  the  program  is 
still  displayed. ) 
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ISPLIT 


INAST 


IRLE 

IRTE 

NSIGN 


ZINNER, 

ZOUTER 

SCALE 


If  ISPLIT  ~ 0,  principal  (that  is,  conventional)  blades 
are  to  be  designed. 

If  ISPLIT  - 1,  splitter  blades  having  the  standard  thick- 
ness distribution  are  to  be  designed. 

If  ISPLIT  = 2,  splitter  blades  having  the  DGA  thickness 
distribution  are  to  be  designed. 

If  INAST  = 0,  no  NASTRAN-compatible  data  is  generated. 

If  |INAST|  = 3,  NASTRAN-compatible  data  is  generated 
using  three-point  averages. 

If  |INAST  | = 4,  NASTRAN-compatible  data  is  generated 
using  four-point  averages. 

If  INAST  is  positive,  both  geometric  and  pressure-load 
data  are  output. 

If  INAST  is  negative,  the  pressure-load  data  is  not  out- 
put, and  will  have  to  be  generated  elsewhere,  such  as 
by  the  aerodynamic  section. 

See  the  Output  Data  description  (Section  IV,  1)  for  further 
details. 

The  computing  station  number  at  the  blade  leading  edge. 

The  computing  station  number  at  the  blade  trailing  edge. 

Indicator  used  to  sign  blade  pressure  forces  according 
to  program  sign  conventions.  For  compressor  rotors, 
if  the  machine  rotates  clockwise  when  viewed  from  the 
front,  set  NSIGN  to  1;  otherwise,  set  NSIGN  to  -1. 

For  compressor  stators,  the  two  values  given  for 
NSIGN  are  reversed. 

The  NZ  manufacturing  sections  are  equi-spaced  between 
z equals  ZINNER  and  ZOUTER. 

When  precision  plots  are  produced,  SCALE  is  the  scale 
factor  employed. 


17 


STACKX 


PLTSZE 

KPTS 


IFANGS 


XSTA 


RSTA 

R 

BLAFOR 


This  is  the  axial  coordinate  of  the  stacking  axis  for  the 
blade,  relative  to  the  same  origin  as  used  for  the 
station  locations,  XSTA. 

The  size  (inches)  of  the  plotter  to  be  used  in  the  creation 
of  precision  plots. 

The  number  of  points  provided  to  specify  the  shape  of  a 
computing  station. 

If  KPTS  = 1,  the  computing  station  is  upright  and  linear. 

If  KPTS  = 2,  the  computing  station  is  linear  and  either 
upright  or  inclined. 

If  KPTS  >2,  a spline  curve  is  fit  through  the  points 
provided  to  specify  the  shape  of  the  station. 

If  IFANGS  = 0,  the  calculations  of  the  quantities  required 
for  aerodynamic  analysis  will  be  omitted  at  a particular 
computing  station. 

If  IFANGS  = 1,  these  calculations  will  be  performed  at 
that  station. 

An  array  of  KPTS  axial  coordinates  (relative  to  an  arbi- 
trary origin)  which,  together  with  RSTA,  specify  the 
shape  of  a particular  computing  station. 

An  array  of  KPTS  radii  which,  together  with  XSTA, 
specify  the  shape  of  a particular  computing  station. 

The  stream  surface  radii  at  NLINES  locations  at  each  of 
the  NSTNS  stations. 

The  difference  in  pressure  between  points  on  the  blade 
pressure  surface  and  suction  surface.  If  zero  is  input, 
the  output  generated  for  the  NASTRAN  program  will 
reflect  no  aerodynamic  loading.  The  pressure  difference 
should  be  in  units  of  force  per  length  squared,  where 
length  is  the  Sc.nie  unit  as  the  dimensions  of  the  blade 
under  consideration 


i8 


«K**N*03 


The  variation  of  properties  of  the  stream  surface  blade 
section  is  specified  as  a function  of  stream  surface 
number.  The  various  quantities  are  then  interpolated 
(or  extrapolated)  at  each  stream  surface.  The  stream  - 
surfaces  are  numbered  consecutively  from  the  inner- 
most outward,  starting  with  1. 0.  ZR  must  increase 
monotonically,  there  being  NSPEC  values  in  all. 

The  blade  irlet  angle. 

The  blade  outlet  angle. 

If  ISECN  = 0,  PP  is  the  ratio  of  the  second  derivative  of 
the  camber  line  at  the  leading  edge  to  its  maximum  value. 
Must  satisfy  -2.  0 < PP  < 1.0. 

If  ISECN  = 1,  PP  is  the  ratio  of  the  second  derivative  of 
the  camber  line  at  the  leading  edge  to  its  maximum  value 
forward  of  the  inflection  point.  Must  satisfy  0.  0 < PP  s 1.0. 

If  ISECN  = 2 or  3,  PP  is  superfluous. 

If  ISECN  = 0,  QQ  is  the  ratio  of  the  second  derivative  of 
the  camber  line  at  the  trailing  edge  to  its  maximum  value. 
Must  satisiy  0.  Ofi  QQ  £ 1.  0. 

If  ISECN  = 1 , QQ  is  the  ratio  of  the  second  derivative  of 
the  camber  line  at  the  trailing  edge  to  its  maximum  value 
rearward  of  the  inflection  point.  Must  satisfy  0.  0 < QQ 
£1.0. 

If  ISECN  = 2 or  3,  QQ  is  superfluous. 

The  ratio  of  blade  leading  edge  radius  to  chord. 

The  ratio  of  blade  maximum  thickness  to  chord. 

The  ratio  of  blade  trailing  edge  half-thickness  to  chord. 

If  ISECN  = 2,  TE  is  superfluous. 

The  location  of  the  blade  maximum  thickness,  as  a fraction 
of  camber  line  length  from  the  leading  edge. 


If  ISECN  = 2,  Z is  superfluous. 
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? 


CORD 


DELX, 

DELY 


S,  BS 


RLES 

TCS 

TES 

ZZS 

PERSPT 

XI,  X2,  X3 
XB 


If  IFCORD  = 0,  CORD  is  the  meridional  projection  of  the 
blade  chord. 

If  IFCORD  = 1,  CORD  is  the  blade  chord. 

The  stacking  axis  passes  through  the  stream  surface  blade 
sections,  offset  from  the  centroids,  leading,  or  tra.'ling 
edge  by  DELX  and  DELY  in  the  x and  y directions 
respectively. 

If  ISECN  = 1 or  3,  S and  BS  are  used  to  specify  the 
locations  of  the  inflection  point  (as  a fraction  of  the 
meridionally-projected  chord  length)  and  the  change  in 
camber  angle  from  the  leading  edge  to  the  inflection  point. 
If  the  absolute  value  of  the  angle  at  the  inflection  point 
is  larger  than  the  absolute  value  of  Bl,  BS  should  have 
the  same  sign  as  Bl,  otherwise,  Bl  and  BS  should  be 
of  opposite  sign. 

The  ratio  of  splitter  section  leading  edge  radius  to 
sp’itter  chord. 

The  ratio  of  splitter  section  maximum  thickness  to 
splitter  chord. 

The  ratio  of  splitter  section  trailing -edge  half-thicknes; 
to  splitter  chord.  If  ISPLIT  = 2,  this  item  is  super- 
fluous. 

The  location  of  the  splitter  section  maximum  thickness, 
as  a fraction  of  splitter  camber  line  length.  If  ISPLIT  = 

2,  this  item  is  superfluous. 

The  ratio  of  desired  splitter  axial  chord  to  main  blade 
axial  chord  on  the  particular  stream  surface. 

These  quantities  are  not  read.  They  may  be  left  blank. 

The  blockage  due  to  the  principal  blades  that  is  to  be 
added  to  that  computed  for  the  splitter  blades.  (Generally, 
cards  containing  these  data  in  the  correct  format  will  have 
been  produced  by  use  of  the  IFANGS  = 1 and  IPUNCH  = 1 
options  when  computing  the  principal  blades.  Take  care 
only  to  incorporate  those  cards  that  correspond  to 
locations  within  the  splitter  blade  envelope.) 
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NRAD 


NDPTS 

NDATR 


NSWITC 


NLE 

NTE 

XKSHPE 

SPEED 

NR 

NTERP 

NMACH 

NLOSS 

NL1 

NL2 

NEVAL 

NCURVE 

NL.ITER 

NDEL 

NOUT1 

NOUT2 

NOUT3 

NBLAD 

R 


The  number  of  radii  at  which  a distribution  of  the 
fraction  of  trailing  edge  deviation  is  input.  Must 
satisfy  1 ■&  NRAD  £ 5. 

1 .e  number  of  points  used  to  define  each  deviation  curve. 
Must  satisfy  1 £ NDPTS  £ 11. 

The  number  of  radii  at  which  an  additional  deviation  angle 
increment  and  the  point  of  maximum  camber  are 
specified.  Must  satisfy  1 £ NDATR  £ 21. 

If  NSWITC  = 1,  the  deviation  correlation  parameter  "m” 
for  the  NACA  (A  g)  meanline  is  used. 

If  NSWITC  = 2,  the  deviation  correlation  parameter  "m" 
for  double -circular-arc  blades  is  used. 

Station  number  at  leading  edge. 

Station  number  at  trailing  edge. 

The  blade  shape  correction  factor  in  the  deviation  rule. 
See  definition  for  Aerodynamic  Section. 

The  number  of  radii  where  a "loss"  is  input. 


See  definition  for  Aerodynamic  Section. 


Radius  at  which  loss  is  specified. 
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XLOSS 


RTE 


DM 


DVFRAC 

RDTE 

DELTAD 


AC 


Loss  description.  The  form  is  prescribed  by  NLOSS; 
see  aerodynamic  section. 

Radius  at  blade  trailing  edge  where  the  following  deviation 
fraction/ chord  curve  applies. 

If  NRAD  = 1,  it  has  no  significance.  Must  increase 
monotonically. 

The  location  on  the  meridional  chord  where  the  deviation 
fraction  is  given.  Expressed  as  a fraction  of  the 
meridional  chord  from  the  leading  edge.  Must  increase 
monotonically. 

Fraction  of  trailing  -edge  deviation  that  occurs  at  location 
DM. 


Radius  at  trailing  edge  where  additional  deviation  and 
point  of  maximum  camber  are  specified. 

Additional  deviation  angle  added  to  that  determined  by 
deviation  rule.  Input  positive  for  conventionally  positive 
deviation  for  both  rotors  and  stators. 

Fraction  of  blade  chord  from  leading  edge  where  maximum 
camber  occurs. 


c. 


TITLE3 

CP 

GASR 

G 


EJ 


NSTNS 


Aerodynamic  Section 

A title  card  for  the  aerodynamic  section  of  the  program. 

Specific  heat  at  constant  pressure.  An  input  value  of 
zero  will  be  reset  to  0.  24,  Units;  H/F/D. 

Gas  constant.  An  input  value  of  zero  will  be  reset  to 
53.32.  Units:  L/SCLFAC/D. 

Acceleration  due  to  gravity.  An  input  value  of  zero  will 
be  reset  to  32.  174.  Units:  L/SCLFAC/T/T. 

Joules  equivalent.  An  input  value  of  zero  will  be  reset 
to  778.  16,  Units:  LF/SCLFAC/H. 

Number  of  computing  stations.  Must  satisfy  3 s NSTNS 
5 30. 
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NSTRMS 

NMAX 


N FORCE 


NBL 


NCASE 


NSPLIT 


NSET1 


NSET2 


Nuriber  of  streamlines.  Must  satisfy  3 £ NSTRMS  s21. 

An  input  value  of  zero  will  be  reset  to  11. 

Maximum  number  of  passes  through  the  iterative  stream- 
line determination  procedure.  An  input  value  of  zero  will 
be  reset  to  40. 

The  first  NFORCE  passes  are  performed  with  arbitrary 
numbers  inserted  should  any  calculation  produce 
impossible  values.  Thereafter,  execution  will  cease, 
the  calculation  having  "failed".  An  input  value  of  zero 
will  be  reset  to  10. 

If  NBL  = 0,  the  annulus  wall  boundary  layer  blockage 
allowance  will  be  held  at  the  values  prescribed  by 
WBLOCK. 

If  NBL  = 1,  blockage  due  to  annulus  wall  boundary  layers 
will  be  recalculated  except  at  station  1.  VISK  and 
SHAPE  are  used  in  the  calculation. 

Number  of  speed/ flow  combinations  to  be  computed. 

Must  satisfy  1 s NCASE  £ 10,  except  that  an  input  value 
of  zero  will  be  reset  to  1. 

If  NSPLIT  = 0,  the  flow  distribution  between  the  stream- 
lines will  be  determined  by  the  program  so  that  roughly 
uniform  increments  of  computing  station  will  occur 
between  the  streamlines  at  station  1. 

If  NSPLIT  = 1,  the  flow  distribution  between  the  stream- 
lines is  read  in  (see  DELF). 

The  blade  loss  coefficient  re-evaluation  option  (specified 
by  NEVAL)  requires  loss  parameter/diffusion  factor 
data.  NSET1  sets  of  data  are  input,  the  set  numbers  being 
allocated  according  to  the  order  in  which  they  are  input. 

Up  to  4 sets  may  be  input  (see  NDIFF). 

When  NLOSS  - 4,  the  loss  coefficients  at  the  station  are 
determined  as  a fraction  of  the  value  at  the  trailing  edge. 
Then,  NSET2  sets  of  curves  are  input  to  define  this 
fraction  at  a function  of  radius  and  meridional  chord.  Up 
to  2 sets  may  be  input  (see  NM). 


NREAD 


NPUNCH 


NPLOT 


NPAGE 

NTRANS 


NMIX 


NMANY 


If  NREAD  = 0,  the  initial  streamline  pattern  estimate 
is  generated  by  the  program. 

If  NREAD  = 1,  the  initial  streamline  pattern  estimate  and 
also  the  DELF  values  are  read  in.  (See  DELF,  R,  X, 
XL.) 

If  NPUNCH  = 0,  no  action  is  taken. 

If  NPUNCH  = 1,  the  final  streamline  pattern  computed 
for  each  running  point  calculation  that  did  not  fail  will 
be  punched.  The  DELF  values  are  also  punched. 

If  NPLOT  = 0,  no  action  is  taken. 

If  NPLOT  = 1,  a CALCOMP  plot  of  the  static  pressure 
distributions  on  hub,  mid,  and  tip  streamlines  will  be 
made  for  each  point. 

If  NPLOT  = 2,  a CALCOMP  plot  of  the  final  streamline 
pattern  will  be  made  for  each  point. 

If  NPLOT  = 3,  the  plots  made  for  NPLOT  = 1 and  2 will 
both  be  made.  (See  XSCALE,  PSCALE,  RLOW,  PLOW.) 

The  maximum  number  of  lines  printed  per  page.  An 
input  value  of  zero  will  be  reset  to  80. 

If  NTRANS  = 0,  no  action  is  taken. 

If  NTRANS  = 1,  relative  total  pressure  loss  coefficients 
will  be  modified  to  account  for  radial  transfer  of  wakes. 
See  Section  V.  1 1 . 

If  NMIX  = 0,  no  action  is  taken. 

If  NMIX  = 1,  entropy,  angular  momentum,  and  total 
enthalpy  distributions  will  be  modified  to  account  for 
turbulent  mixing.  See  Section  V.  12. 

The  number  of  computing  stations  for  which  blade 
descriptive  data  is  being  generated  by  the  analytic  mean- 
line section. 
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NSTPLT 


NEQN 


i 

l 


If  NSTPLT  = 0,  no  action  is  taken. 

If  NSTPLT  = 1,  a line-printer  plot  of  the  changes  made 
to  the  midstreaxnline  ‘-t‘  coordinate  is  made  for  each 
computing  station.  If  more  than  59  passes  through  the 
iterative  procedure  have  been  made,  then  the  plots  will 
show  the  changes  for  the  last  59  passes.  The  graph 
should  decay  approximately  exponentially  towards  zero, 
indicating  that  the  streamline  locations  are  stabilizing. 
Decaying  oscillations  are  equally  acceptable,  but,  growing 
oscillations  show  the  need  for  heavier  damping  in  the 
streamline  relocation  calculations,  that  is,  a decrease 
in  R CONST. 

This  item  controls  the  selection  of  the  form  of  momentum 
equation  that  will  be  used  to  compute  the  meridional 
velocity  distribution  at  each  computing  station.  There  are 
two  basic  forms,  and  for  each  case,  one  may  select  not 
to  compute  the  tern.s  relating  to  blade  forces.  (See  also 
Section  V.  1 . ) 

If  NEQN  = 0,  the  momentum  equation  involves  the  differ- 
ential form  of  the  continuity  equations  and  hence  (l-Mrn) 
terms  in  the  denominator.  Streamwise  gradients  of 
entropy  and  angular  momentum  (blade  forces)  are  com- 
puted within  blades  and  at  the  blade  edges  (provided  data 
that  describe  the  blades  are  given).  Elsewhere,  stream  - 
wise  entropy  gradients  only  are  included  in  a sin. pier  form 
of  the  momentum  equation,  except  that  at  the  first  and  last 
computing  station,  all  streamwise  gradients  are  taken  to 
be  zero.  This  is  generally  the  preferred  option  when 
computing  stations  are  located  within  the  blade  rows. 

If  NEQN  = 1,  the  momentum  equation  form  is  similar  to 
that  used  when  NEQN  = 0,  but  angular  momentum  gradients 
(blade  force  terms)  are  nowhere  computed.  This  generally 
is  the  preferred  option  when  computing  stations  are  located 
at  the  blade  edges  only. 

If  NEQN  = 2,  the  momentum  equation  includes  an  explicit 
dVm/dm  term  instead  of  the  (I- Mm)  denominator  terms. 
All  streamwise  gradients  (including  blade  force  terms) 
are  computed  as  for  the  case  NEQN  = 0.  When  computing 
stations  are  located  within  the  blade  rows,  the  results 
will  generally  be  similar  to  those  obtained  with  NEQN  = 0, 
and  solutions  may  be  found  that  cannot  be  computed  with 
NEQN  = 0 due  to  high  meridional  Mach  numbers. 
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N WHICH 

SCLFAC 

TOLNCE 

VISK 

SHAPE 

XSCALE 

PSCALE 

RLOW 

PLOW 

XMMAX 


If  NEQN  = 3,  the  momentum  equation  is  similar  to  that 
used  when  NEQN  = 1,  but  (as  for  the  case  NEQN  = 1)  no 
angular  momentum  gradients  are  computed.  This  may  be 
used  when  computing  stations  are  located  only  at  the  blade 
edges  and  high  meridional  Mach  numbers  preclude  the 
use  of  NEQN  = 1. 

The  numbers  of  each  of  the  computing  stations  for  which 
blade  descriptive  data  is  being  generated  by  the  analytic 
meanline  section. 

Linear  dimension  scale  factor.  An  input  value  of  zero 
will  be  reset  to  12.  0. 

Basic  tolerance  in  iterative  calculation  scheme.  An  input 
value  of  zero  will  be  reset  to  0.  001.  (See  discussion  of 
tolerance  scheme  in  Section  VI. ) 

Kinematic  viscosity  of  gas  (for  annulus  wall  boundary 
layer  calculations).  An  input  value  of  zero  will  be  reset 
to  0.00018.  Units:  LL/SCLFAC/SCLFAC/T. 

Shape  factor  for  annulus  wall  boundary  layer  calculations. 
An  input  value  of  zero  will  be  reset  to  0.  7. 

The  scale  used  for  the  physical  dimensions  in  the  static 
pressure  and  streamline  plots  given  as  length  unit  per 
inch  of  plot. 

The  scale  used  to  plot  static  pressure  given  as  pressure 
unit  per  inch  of  plot. 

Minin. um  value  of  radius  shown  on  streamline  plots. 

Units:  L. 

Minimum  value  of  pressure  shown  on  static  pressure 
plots.  Units:  F/L/L. 

The  square  of  the  Mach  number  that  appears  in  the  equation 
for  the  streamline  relocation  relaxation  factor  is  limited 
to  be  not  greater  than  XMMAX.  Thus,  at  computing  sta- 
tions where  the  appropriate  Mach  number  is  high  enough 
for  the  limit  to  be  imposed,  a decrease  in  XMMAX 
corresponds  to  an  increase  in  damping.  If  a value  of  zero 
is  input,  it  is  reset  to  0.  6. 
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RCONST 


The  constant  in  the  equation  for  the  streamline  relocation 
relaxation  factor.  The  value  of  8.  0 that  the  analysis  yields 
is  often  too  high  for  stability.  If  zero  is  input,  it  is  reset 
to  6.  0. 


CONTR 

CONMX 


The  constant  in  the  blade  wake  radial  transfer  calculations. 

The  eddy  viscosity  for  the  turbulent  mixing  calculations. 
Units:  L2/SCLFAC2/T. 


FLOW 

SPDFAC 


Compressor  flow  rate.  Units:  F/T. 

The  speed  of  rotation  of  each  computing  station  is  SPDFAC 
times  SPEED  (I).  The  units  for  the  product  are  revolutions/ 
(60xT). 


NSPEC 


The  number  of  points  used  to  define  a computing  station. 
Must  satisfy  2 £ NSPEC  £ 21,  and  also  the  sum  of  NSPEC 
for  all  stations  £ 150.  If  2 points  are  used,  the  station  is 
a straight  line.  Otherwise,  a spline-curve  is  fitted 
through  the  given  points. 


XSTN,  RSTN 


The  axial  and  radial  coordinates,  respectively,  of  a point 
defining  a computing  station.  The  first  point  must  be  on 
the  hub  and  the  last  point  must  be  on  the  casing.  Units:  L. 


K I 


E i 


NDATA 


NTERP 


Number  of  points  defining  conditions  or  blade  geometry  at 
a computing  station.  Must  satisfy  0 £ NDATA  £ 21,  and 
also  the  sum  of  NDATA  for  ail  stations  £ 100. 

If  NTERP  = 0,  and  NDATA  a 3,  interpolation  of  the  data 
at  the  station  is  by  spline-fit. 


I ! 


I ! 


NDIMEN 


If  NTERP  = 1 (or  NDATA  £ 2),  interpolation  is  linear 
point-to-point. 

If  NDIMEN  = 0,  the  data  are  input  as  a function  of  radius. 


If  NDIMEN  = 1,  the  data  are  input  as  a function  of  radius 
normalized  with  respect  to  tip  radius. 


If  NDIMEN  = 2,  the  data  are  input  as  a function  of  distance 
along  the  computing  station  from  the  hub. 
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NMACH 


DATAC 


DATA1 
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If  NDIMEN  = 3„  the  data  are  input  as  a function  of 
distance  along  the  computing  station  normalized  with 
respect  to  the  total  computing  station  length. 

If  NMACH  = 0,  the  subsonic  solution  to  the  continuity 
equation  is  sought. 

If  NMACH  = 1,  the  supersonic  solution  to  the  continuity 
equation  is  sought.  This  should  only  be  used  at  stations 
where  the  relative  flow  angle  is  specified,  that  is, 

NWCRK  = 5,  6,  or  7. 

The  coordinate  on  the  computing  station,  defined  according 
to  NDIMEN,  where  the  following  data  items  apply.  Must 
increase  n.onotonically.  For  dimensional  cases,  units  are 
L. 


At  Station  1 and  if  NWCRK  = 1,  DATA1  is  total  pressure. 
Units:  F/L/L. 

If  NWORK  = 0 and  the  station  is  at  a blade  leading  edge, 
by  setting  NDATA  ^ 0,  the  blade  leading  edge  may  be 
described.  Then  DATA1  is  the  blade  angle  measured  in 
the  cylindrical  plane.  Generally  negative  for  a rotor, 
positive  for  a stator.  (Define  the  blade  lean  angle  (DATA3) 
also.)  Units:  A. 

If  NWORK  = 2,  DATA1  is  total  enthalpy.  Units:  H/F. 

If  NWORK  = 3,  DATA1  is  angular  momentum  (radius  tin.es 
absolute  whirl  velocity).  Units:  LL/SCLFAC/T. 

If  NWORK  = 4,  DATA1  is  absolute  whirl  velocity.  Units: 
L/SCLFAC/T. 

If  NWORK  = 5,  DATA1  is  blade  angle  measured  ip  the 
stream. surface  plane.  Generally  negative  for  a rotor, 
positive  for  a stator.  If  zero  deviation  is  input,  if  becomes 
the  relative  flow  angle.  Units:  A. 

If  NWORK  - 6,  DATA1  is  the  blade  angle  measured  in  the 
cylindrical  plane.  Generally  negative  for  a rotor,  positive 
for  a stator.  If  zero  deviation  is  input,  it  becomes,  after 
correction  for  stream  surface  orientation  and  station  lean 
angle,  the  relative  flow  angle.  Units;  A. 
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If  NWORK  = 7,  DATA1  is  the  reference  relative  outlet 
flow  angle  measured  in  the  stream  surface  plane.  Generally 
negative  for  a rotor,  positive  for  a stator.  Units:  A. 

At  Station  1,  DATA2  is  total  temperature.  Units:  D. 

If  NLOSS  = 1,  BATA2  is  the  relative  total  pressure  loss 
coefficient.  The  relative  total  pressure  loss  is  n»easured 
from  the  station  that  is  NL1  stations  removed  from  the 
current  station,  NL1  being  negative  to  indicate  an 
upstream  station.  The  relative  dynamic  head  is  determined 
NL2  stations  removed  from  the  current  station,  positive 
for  a downstream  station,  negative  for  an  upstream  station. 

If  NLOSS  = 2,  DATA2  is  the  isentropic  efficiency  of 
compression  relative  to  condition^  NL1  stations  removed, 
NL1  being  negative  to  indicate  an  upstream  station. 

If  NLOSS  = 3,  DATA2  is  the  entropy  rise  relative  to  the 
value  NL1  stations  removed,  NL1  being  negative  to 
indicate  an  upstream  station.  Units:  H/F/D. 

If  NLOSS  = 4,  DATA2  is  not  used,  but  a relative  total 
pressure  loss  coefficient  is  determined  from  the  trailing 
edge  value  and  curve  set  number  NCURVE  of  the  NSET2 
families  of  curves.  NL1  and  NL2  apply  as  for  NLOSS  = 1. 

If  NWORK  = 7,  DATA  2 is  the  reference  (minimum) 
relative  total  pressure  loss  coefficient.  NL1  and  NL2 
apply  as  for  NLOSS  = 1. 

The  blade  lean  angle  measured  from  the  projection  of  a 
radial  line  in  the  plane  of  the  computing  station,  positive 
when  the  innermost  portion  of  the  blade  precedes  the 
outermost  in  the  direction  of  rotor  rotation.  Units:  A. 

The  fraction  of  the  periphery  that  is  blocked  by  the  presence 
of  the  blades. 

Cascade  solidity.  When  a number  of  stations  are  used  to 
describe  the  flow  through  a blade,  values  are  only  required 
at  the  trailing  edge.  (They  are  used  in  the  loss  coefficient 
re-estimation  procedure,  and  to  evaluate  diffusion  factors 
for  the  output.) 
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DATA6 


DATA7 

DATA8 


DATA  9 


NWORK 


If  NWORK  = 5 or  6,  DATA6  is  the  deviation  angle 
measured  in  the  streamsurface  plane.  Generally  negative 
for  a rotor,  positive  for  a stator.  Units:  A. 

If  NWORK  = 7,  DATA6  is  reference  relative  inlet  angle, 
to  which  the  minimum  loss  coefficient  (DATA2)  and  the 
reference  relative  outlet  angle  (DATA7)  correspond. 
Measured  in  the  streamsurface  plane  and  generally 
negative  for  a rotor,  positive  for  a stator.  Units:  A. 

If  NWORK  = 7,  DATA?  is  the  rate  of  change  of  relative 
outlet  angle  with  relative  inlet  angle. 

If  NWORK  = 7,  DATA8  is  the  relative  inlet  angle  larger 
than  the  reference  value  at  which  the  loss  coefficient  attains 
twice  its  reference  value.  Measured  in  the  stream, surface 
plane.  Units:  A. 

If  NWORK  = 7,  DATA9  is  the  relative  inlet  angle  smaller 
than  the  reference  value  at  which  the  loss  coefficient  attains 
twice  its  reference  value.  Measured  in  the  streams  surface 
plane.  Units:  A. 

If  NWORK  = 0,  constant  entropy,  angular  momentum,  and 
total  enthalpy  exist  along  streamlines  from  the  previous 
station.  (If  NMIX  = 1,  the  distributions  will  be  modified.) 

If  NWORK  = 1,  the  total  pressure  distribution  at  the  com- 
puting station  is  specified.  Use  for  rotors  only. 

If  NWORK  = 2,  the  total  enthalpy  distribution  at  the  com- 
puting station  is  specified.  Use  for  rotors  only. 

If  NWORK  = 3,  the  absolute  angular  momentum  distribution 
at  the  ccruputing  station  is  specified. 

If  NWORK  = 4,  the  absolute  whirl  velocity  distribution  at 
the  computing  station  is  specified. 

If  NWORK  = 5,  the  relative  flow  angle  distribution  at  the 
station  is  specified  by  giving  blade  angles  and  deviation 
angles,  both  measured  in  the  streamsurface  plane. 

If  NWORK  = 6,  the  relative  flow  angle  distribution  at  the 
station  is  specified  by  giving  the  blade  angles  measured 
in  the  cylindrical  plane,  and  the  deviation  angles  measured 
in  the  streamsurface  plane. 
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NLOSS 


NL1 


NL2 


NEVAL 


If  NWORK  = 7,  the  relative  flow  angle  and  relative  total 
pressure  loss  coefficient  distributions  are  specified  by 
means  of  an  off-design  analysis  procedure.  "Reference", 
"stalling",  and  "choking"  relative  inlet  angles  are 
specified.  The  minimum  loss  coefficient  varies  para- 
bolically  with  the  relative  inlet  angle  so  that  it  is  twice 
the  minimum  value  at  the  "stalling"  or  "choking"  values. 

A maximum  value  of  0.  5 is  imposed.  "Reference" 
relative  outlet  angles  and  the  rate  of  change  of  outlet 
angle  with  inlet  angle  are  specified,  and  the  relative 
outlet  angle  varies  linearly  from  the  reference  value 
with  the  relative  inlet  angle.  NLOSS  should  be  set  to  zero. 

I‘  . 'LOSS  = 1,  the  relative  total  pressure  loss  coefficient 
distribution  is  specified. 

If  NLOSS  = 2,  the  isentropic  efficiency  (for  compression) 
distribution  is  specified. 

If  NLOSS  = 3,  the  entropy  rise  distribution  is  specified. 

If  NLOSS  = 4,  the  total  pressure  loss  coefficient  distribution 
is  specified  by  use  of  curve-set  NCURVE  of  the  NSET2 
families  of  curves  giving  the  fraction  of  final  (trailing 
edge)  loss  coefficient. 

The  station  from  which  the  loss  (in  whatever  form  NLOSS 
specifies)  is  measured,  is  NL1  stations  removed  from  the 
station  being  evaluated.  NL1  is  negative  to  indicate  an 
upstream  station. 

When  a relative  total  pressure  loss  coefficient  is  used  to 
specify  losses,  the  relative  dynamic  head  is  taken  NL2 
stations  removed  from  the  station  being  evaluated.  NL2 
may  be  positive,  zero,  or  negative;  a positive  value 
indicates  a downstream  station,  a negative  value  indicates 
an  upstream  station. 

If  NEVAL  = 0,  no  action  is  taken. 

If  NEVAL  >0,  curve-set  number  NEVAL  of  the  NSETI 
families  of  curve  giving  diffusion  loss  parameter  as  a 
function  of  diffusion  factor  will  be  used  to  re-estimate 
the  relative  total  pressure  loss  coefficient.  NLOSS  must 
be  1,  and  NL1  and  NL2  must  specify  the  leading  edge  of 
the  blade.  See  also  NDEL. 
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If  NEVAL  <0,  curve-set  number  | NEVAL  I is  usee? 

NEVAL  > 0,  except  that  the  re -estimation  is  only  ma-e 
after  the  overall  computation  is  completed  (with  the  input 
losses).  The  resulting  loss  coefficients  are  displayed  but 
not  incorporated  into  the  overall  calculation.  See  also 
NDEL. 

When  NLOSS  = 4,  curve-set  NCURVE  of  the  NSET2 
families  of  curves,  specifying  the  fraction  of  trailing- 
edge  loss  coefficient  as  a function  of  meridional  chord 
is  used. 

When  NEVAL  > 0,  up  to  NLITER  re-estimations  of  the  loss 
coefficient  will  be  n.ade  at  a given  station  during  any  one 
pass  through  the  overall  iterative  procedure.  Less  than 
NLITER  re-estimations  will  be  made  if  the  ve'ocity 
profile  is  unchanged  by  re-estimating  the  loss  coefficients. 
(See  discussion  of  tolerance  scheme  in  Section  VI.) 

When  NEVAL  = 0,  set  NDEL  to  0.  When  NEVAL  f 0.  and 
NDEL  > 0,  a component  of  the  re-estimated  loss  coefficient 
is  a shock  loss.  The  relative  inlet  Mach  number  is 
expanded  (or  com  pressed)  through  a Prandtl -Meyer 
expansion  on  the  suction  surface,  and  NDEL  is  the  number 
of  points  at  which  the  Prandtl-Meyer  angle  is  given.  If 
NDEL  - 0,  the  shock  loss  is  set  to  zero.  Must  satisfy 
0 s NDEL  £ 21,  and  also  the  sum  of  NDEL  tor  all  stations 
£ 100. 

If  NOUT1  = 0,  no  action  is  taken. 

If  NOUT1  = 1,  cards  will  be  punched  that  may  be  incorporated 
into  the  data  deck  for  a subsequent  run  of  the  analytic  mean- 
line section.  These  are  the  (NLINES  times  NSTNS)  cards 
specifying  R. 

If  NOUT2  = 0,  no  action  is  taken. 

If  NOUT2  = 1,  output  records  are  created  that  may  be  used 
as  input  to  the  arbitrary  meanline  section.  These  are 
records  specifying  KPTS,  XSTA,  RSTA,  R,  and  AIRANG. 

If  NARBIT  = 0,  the  records  will  be  issued  ?.s  punch  cards; 

If  NARBIT  = 1,  the  records  will  be  stored  within  the 
computer  on  file  LOG6. 
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NOUT3  This  data  :iem  controls  the  generation  of  NASTRAN- 

compatibli  pressure  difference  output  for  use  in  a 
subseque-*  blade  stress  analysis.  For  details  of  the 
triangular  mesh  that  is  used,  see  the  Output  Description 
in  Section  IV.  2.  c. 
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If  NOUT3  = 1,  the  station  is  at  a blade  leading  edge. 

If  NOUT3  = 2,  the  station  is  at  a blade  trailing  edge. 

If  NOUT3  = 3,  the  station  is  at  the  trailing  edge  of  one 
blade,  and  at  the  leading  edge  of  another. 

If  NOUT3  = 0,  the  station  may  be  between  blade  rows, 
or  within  a blade  row  for  which  output  is  required, 
depending  upon  the  use  of  NOUT3  ^ 0 elsewhere.  See 
also  description  of  NBLADE  below. 


NBLADE  This  item  is  used  in  determining  the  pressure  difference 

across  the  blade.  The  number  of  blades  is  j NBLADE  J. 

If  NBLADE  is  positive,  ’three-point  averaging"  is  used 
to  determine  the  pressure  difference  across  each  blade 
element.  If  NBLADE  is  negative,  "four-point  averaging" 
is  used.  (See  the  Output  Description  in  Section  IV.  2.  c.) 

If  NBLADE  is  input  as  zero,  a value  of  +10  is  used.  At 
a leading  edge,  the  value  for  the  following  station  is  used: 
elsewhere  the  value  at  a station  applies  to  the  interv,  1 
upstream  of  the  station.  Thus  by  varying  the  «ign  of 
NBLADE,  the  averaging  method  used  for  the  pressure  forces 
may  be  varied  for  different  axial  segments  of  a blade  row. 

SPEED  This  card  is  omitted  if  NDATA  = 0.  The  speed  of  rotation 

of  the  blade.  At  a blade  leading  edge,  it  should  be  set  to 
zero.  The  product  SPDFAC  tim.es  SPEED  has  units  of 
revolutions/(  T x 60). 

DELC  The  coordinate  at  which  Prandtl-Meyer  expansion  angles 

are  given.  It  defines  the  angle  as  a function  of  the 
dimensions  of  the  leading  edge  station,  in  the  manner 
specified  by  NDLMEN  for  the  current,  that  is  trailing 
edge  station.  Must  increase  monotonic=»lly.  For  dimen- 
sional cases,  units  are  L. 
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DELTA 


W BLOCK 


BBLOCK, 

BDIST 


NDIFF 

DIFF 

FDHUB 

FDMID 

FDTIP 

NM 

NRAD 


The  Prandtl-Meyer  expansion  angles.  A positive  value 
implies  expansion.  If  blade  angles  are  given  at  the  leading 
edge,  the  incidence  angles  are  added  to  the  value  specified 
by  DELTA.  Units:  A.  (Blade  angles  are  measured  in 
the  cylindrical  plane. ) 

A blockage  factor  that  is  incorporated  into  the  continuity 
equation  to  account  for  annulus  wall  boundary  layers.  It 
is  expressed  as  the  fraction  of  total  area  at  the  computing 
station  that  is  blocked.  If  NBL  = 1,  values  (except  at 
Station  1)  are  revised  during  computation,  involving  data 
items  VISK  and  SHAPE. 

A blockage  factor  is  incorporated  into  the  continuity 
equation  that  may  be  used  to  account  for  blade  wakes  or 
other  effects.  It  varies  linearly  with  distance  along  the 
computing  station.  EBLOCK  is  the  value  at  mid-station 
(expressed  as  the  fraction  of  the  periphery  blocked),  and 
BDIST  is  the  ratio  of  the  value  on  the  hub  to  the  mid- 
value. 

When  NSET1>  0,  there  are  NDIFF  points  defining  loss 
diffusion  parameter  as  a function  of  diffusion  factor. 

Must  satisfy  1 s NDIFF  £ 15. 

The  diffusion  factor  at  which  loss  parameters  are  specified. 
Must  increase  monotonically. 

Diffusion  loss  parameter  at  10  per  cent  of  the  radial  blade 
height. 

Diffusion  loss  parameter  at  50  per  cent  of  the  radial  blade 
height. 

Diffusion  loss  parameter  at  90  per  cent  of  the  radial  blade 
height. 

When  NSET2>  0,  there  are  NM  points  defining  the  fraction 
of  trailing  edge  loss  coefficient  as  a function  of  meridional 
chord.  Must  satisfy  1 s NM  ^ 11. 

The  number  of  radial  locations  where  NM  loss  fraction/ 
chord  points  are  given.  Must  satisfy  1 s NRAD  £ 5. 


*i 


PIWJ 


E » 


f1?* 


1^^ 


TERAD 
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The  fraction  of  radial  blade  height  at  the  trailing  edge 
where  the  following  loss  fraction/ chord  curve  applies. 
If  NRAD  = 1,  it  has  no  significance. 


assess* 


DM  The  location  on  the  meridional  chord  where  the  loss 

fraction  is  given.  Expressed  as  a fraction  of  meridional 
chord  from  the  leading  edge.  Must  increase  mono- 
tonically. 


WFRAC  Fraction  of  trailing  edge  loss  coefficient  that  occurs  at 

location  DM. 


DELF 


R 


X 

XL 


The  fraction  of  tha  total  flow  that  is  to  occur  between  the 
hub  and  each  streamline.  The  hub  and  casing  are  included, 
so  that  the  first  value  must  be  0.  0,  and  the  last  (NSTRM) 
value  must  be  1.0. 

Estimated  streamline  radius.  (These  data  are  input  from 
hub  to  tip  for  the  first  station,  from  hub  to  tip  for  the 
second  station,  and  so  on. ) Units:  L. 

Estimated  axial  coordinate  at  intersection  of  streamline 
with  computing  station.  Units:  L. 

Estimated  distance  along  computing  station  from  hub  to 
intersection  of  streamline  with  computing  station. 

Units:  L. 


II,  JJ  Station  end  streamline  number.  These  are  merely  read 

in  and  printed  out  to  give  a check  or  the  order  of  the  cards, 
presumed  to  have  been  punched  on  p evious  run. 

d.  Arbitrary  Meanline  Bla  Section 

For  a more  detailed  discussion  of  the  input  tc  this  section,  see 
Reference  4.  In  this  section,  the  dimensioned  input  is  either  in  degrees 
(A),  or  length  (L). 


TITLE4 

N LINES 

NSTNS 

NZ 


A title  card  for  the  arbitrary  meanline  section  of  the 
program. 

As  for  Analytic  Meanline  Section. 

As  for  Analytic  Meanlinc  Section. 

As  for  Analytic  Meanlinc  Section. 
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NSPEC 

ISEGPT 


N BLADE 

ISTAK 

IPUNCH 

IFPLOT 

IPRINT 

INAST 

ZINNER 

ZOUTER 

SCALE 

STACKX 

PLTSZE 

IRLE 

IRTE 
NR A DEV 

NINC 

NSIGN 


As  for  Analytic  Meanline  Section. 

The  number  of  points  to  define  each  surface  in  each  blade 
segment,  including  the  end-points  common  to  the  adjacent 
segments.  Must  satisfy  2^  ISEGPT  and  also  (ISEGPT-1) 
times  the  number  of  segments  £ 80.  (The  number  of 
segments  is  IRTE  - IRLE  + 1.) 

As  for  Analytic  Meanline  Section. 

As  for  Analytic  Meanline  Section. 

As  for  Analytic  Meanline  Section. 

As  for  Analytic  Meanline  Section. 

As  for  Analytic  Meanline  Section. 

As  for  Analytic  Meanline  Section. 

As  for  Analytic  Meanline  Section. 

As  for  Analytic  Meanline  Section. 

As  for  Analytic  Meanline  Section. 

As  for  Analytic  Meanline  Section. 

As  for  Analytic  Meanline  Section. 

Station  number  at  blade  leading  edge.  Note  that  the  first 
station  specified  ir.  .ne  arbitrary  meanline  section  is 
number  1. 

Station  number  at  blade  trailing  edge. 

The  number  of  radii  at  which  a distribution  of  the  fraction 
of  trailing  edge  deviation  is  inp'.i..  Must  satisfy 
1 £ NRADEV  s 5. 

The  number  of  radii  at  which  the  incidence  angle  is 
specified.  Must  satisfy  1 £ NINC  £ 5. 

This  is  used  to  establish  the  sign  convention.  Normally, 

+ 1 for  stators  and  -1  for  rotors. 
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IFCA 

XKSHPE 

SOLTOL 

NPTS 

RADEV 

SM 

DEVCRV 

RINC 

XINC 

DELDEV 

I FANGS 

XPTS 

XSTA 

RSTA 

R 


* 
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If  IFCA  = 1,  the  deviation  correlation  parameter  "m" 
for  the  NACA  (A  ) meanline  is  used. 

If  IFCA  - 2.  the  deviation  correlation  parameter  "m" 
for  double-circular-arc  blades  is  used. 

The  blade  shape  correction  factor  in  the  deviation  rule. 

Solidity  tolerance  in  iteration  for  deviation  angle. 

The  number  of  points  defining  the  deviation  curve. 

Radius  at  the  blade  tiailing-edge  where  the  following 
deviation  fraction/ chord  curve  applies.  If  NRADEV  = 1, 
it  has  no  significance.  Must  increase  monotonically. 

The  location  on  the  meridional  chord  where  the  deviation 
fraction  is  given.  Expressed  as  a fraction  of  the  meridional 
chord  from  the  leading  edge.  Must  increase  monotonically. 

Fraction  of  trailing-edge  deviation  that  occurs  at  location 
SM. 


Radius  at  which  incidence  angle  and  additional  deviation 
angle  is  given. 

Incidence  angle  at  radius  RINC.  Input  positive  for  conven 
tionaily-positive  incidence  for  both  rotors  and  stators. 

Additional  deviation  angle  added  to  that  determined  by 
deviation  rule.  Input  positive  for  conventionally  positive 
deviation  for  both  rotors  and  stators.  Applied  on  the 
stream  surface  that  passes  through  the  leading  edge 
station  at  radius  RINC. 

As  for  Analytic  Meanline  Section. 

As  for  Analytic  Meanline  Section. 

As  for  Analytic  Meanline  Section. 

As  for  Analytic  Meanline  Section. 

Streamsurface  radius. 
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AIRANG 


Relative  flow  angle  at  radius  R.  Normally  negative  for 
rotors,  positive  for  stators. 

BLAFOR  As  for  Analytic  Meanline  Section. 


ZR 


YA 


RLE 


The  variations  of  properties  of  the  stream  surface  blade 
sections  are  specified  as  a function  of  stream  surface 
number.  The  stream  surfaces  are  numbered  consecutively 
from  the  innermost  outward.  Must  increase  monotonically. 

The  chord  is  multiplied  by  the  factor  (1-YA)  when  the 
solidity  is  determined,  which  is  then  used  in  the  deviation 
angle  calculations. 

The  ratio  of  blade  leading-edge  radius  to  the  chord. 


TC 


The  ratio  of  blade  maximum  thickness  to  chord. 


TE 

ZZ 

DELX, 

DELY 

XI,  X2 
X3,  XB 


The  ratio  of  blade  trailing-edge  half-thickness  to  chord. 

The  point  of  blade  maximum  thickness,  as  a fraction  of 
camber  line  length  from  the  leading  edge. 

As  for  Analytic  Meanline  Section. 

As  for  Analytic  Meanline  Section. 
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SECTION  IV 


OUTPUT  DATA 

1.  ANALYTIC  MEANLINE  SECTION 

Printed  output  may  be  considered  to  consist  of  four  sections;  a print- 
out of  the  input  data,  details  of  the  blade  sections  on  each  streamsurface,  a 
listing  of  quantities  required  for  aerodynamic  analysis,  and  details  of  the 
manufacturing  sections  determined  on  the  constant-z  planes.  These  are 
briefly  described  below.  In  the  explanation  which  follows,  parenthetical 
statements  are  understood  to  refer  to  the  particular  case  of  the  double- 
circular-arc blade  (ISECN  = 2). 

The  input  data  printout  includes  all  quantities  read  in,  and  is  self- 
explanato  ry. 

Details  of  the  streamsurface  blade  sections  are  printed  if  IPRINT  = 

0 or  1.  Listed  first  are  the  parameters  defining  the  blade  section.  These 
are  interpolated  at  the  streamsurface  from  the  tables  read  in.  Then  follow 
details  of  the  blade  section  in  "normalized"  form.  The  blade  section  geometry 
is  given  for  the  section  specified,  except  that  the  meridional  projection  of 
the  chord  is  unity.  For  this  section  of  the  output,  the  coordinate  origin  is 
the  blade  leading  edge.  The  following  quantities  are  given:  blade  chord; 
stagger  angle;  can  ber  angle;  section  area;  location  of  the  centroid  of  the 
section;  second  momen  t of  area  of  the  section  about  the  centroid;  orienta- 
tion of  the  principal  axes;  and  the  principal  second  moments  of  area  of  the 
section  about  the  centroid.  Then  are  listed  the  coordinates  of  the  camber 
line,  the  camber  line  angle,  the  section  thickness,  and  the  coordinates  of 
the  blade  surfaces.  NPOINT  values  are  given. 

A lineprinter  plot  of  the  normalized  section  follows.  The  scales  for 
the  plot  are  arranged  so  that  the  section  just  fills  the  page,  so  that  the 
scales  will  generally  differ  from  one  plot  to  another.  "Dimensional"  details 
of  the  blade  section  are  given  next.  The  normalized  data  given  previously  is 
scaled  to  give  a blade  section  as  defined  by  IFCORD  and  CORD.  For  this 
section  of  the  output,  the  coordinates  are  with  respect  to  the  blade  stacking 
axis.  The  following  quantities  are  given:  blade  chord;  radius  and  location 
of  center  of  leading  (and  trailing)  edge(s);  section  area,  the  second  moments 
of  area  of  the  section  about  the  centroid  and  the  principal  second  moments 
of  area  of  the  section  about  the  centroid.  The  coordinates  of  NPOINT  points 
on  the  blade  surfaces  are  then  listed,  followed  by  the  coordinates  of  31 
points  distributed  at  (roughly)  six  degree  intervals  around  the  leading  (and 
trailing)  edges.  Finally,  the  coordinates  of  the  blade  surfaces  and  points 
around  the  leading  (and  trailing)  edge(s)  is  (are)  shown  in  Cartesian  form. 
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The  quantities  required  for  aerodynamic  analysis  are  printed  at  all 
computing;  stations  specified  by  the  IFANGS  parameter.  The  radius,  blade 
section  angle,  blade  lean  angle,  blade  blockage,  and  relative  angular 
location  of  the  camber  line  are  printed  at  each  stream  surface  intersection 
with  the  particular  computing  station.  The  blade  section  angle  is  measured 
in  the  cylindrical  plane,  and  the  blade  lean  angle  is  measured  in  the  constant- 
axial-coordinate  plane. 

Details  of  the  n.anufacturing  sections  are  printed  if  IPRINT  = 0 or  2. 

At  each  value  of  z specified  by  ZINNER,  ZOUTER,  and  NZ,  section 
properties  and  coordinates  are  given.  The  origin  for  the  coordinates  is  the 
blade  stacking  axis.  The  following  quantities  are  given:  section  area;  the 
location  of  the  centroid  of  the  section;  the  second  moments  of  area  of  the 
section  about  the  centroid;  the  principal  second  moments  of  area  of  the 
section  about  the  centroid;  the  orientation  of  the  principal  axes;  and  the 
section  torsional  constant.  Then  the  coordinates  of  NPOINT  points  on  the 
blade  section  surfaces  are  listed,  followed  by  31  points  around  the  leading 
(and  trailing)  edge(s). 

If  NAERO  = 1,  the  additional  input  and  output  required  for,  and 
generated  by,  the  interface  are  also  printed.  (Apart  from  the  input  data 
printout,  this  is  the  only  printed  output  when  IPRINT  = 3.) 

If  IPUNCH  = 1 and  NAERO  = 0,  the  program  punches  the  quantities 
required  for  aerodynamic  analysis,  together  with  identifying  indices  denoting 
station  number  and  streamsurface  number,  on  cards  in  the  following  format: 

5 fields  each  of  12  locations  for  the  quantities  themselves,  followed  by  2 
fields  each  of  3 locations  for  the  indices. 

If  INAST  ^0,  cards  are  punched  that  may  be  used  as  input  for  the 
NASTRAN  stress  analysis  program.  For  the  purpose  of  stress  analysis, 
the  blade  is  divided  into  a number  of  triangular  elements,  each  defined  by 
three  grid  points.  The  intersections  between  computing  stations  and 
stream  surfaces  are  used  as  the  grid  points,  and  the  grid  points  and  element 
numbering  scheme  adopted  through  Program  UD0300  is  illustrated  in 
Fig  are  1. 

The  NASTRAN  input  data  format  includes  cards  identified  by  the  codes 
GRID,  CTRIA2,  PTRLA2,  and  PLOAD2.  The  data  are  fully  described  in 
Reference  7,  but  briefly,  the  GRID  cards  each  define  a grid  point  number 
and  give  the  coordinates  at  the  grid  point,  the  CTRLA2  cards  each  define  an 
element  in  terms  of  the  three  appropriate  grid  points  (by  number,  and  in  a 
significant  order),  the  PTRLA2  cards  each  give  an  average  blade  thickness 
for  an  element,  and  the  PLOAD2  cards  each  give  an  average  pressure 
loading  for  an  element.  If  ] INAST  | = 3,  the  thicknesses  and  pressures  are 
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Figure  1.  NASTRAN  Grid  Point  and  Element  Numbering 
Scheme. 


obtained  by  averaging  the  values  at  the  three  grid  points  that  comprise  the 
blade  elemient.  Because  of  the  different  weighting  towards  the  "front" 
and  "rear"  of  the  elements,  the  averages  computed  do  not  always  vary 
sensibly,  especially  near  the  blade  edges,  and  therefore  an  alternative 
averaging  procedure  is  also  available.  When  IlNAST  | = 4,  the  averages 
are  formed  from  values  at  the  four  grid  points  which  are  rather  obviously 
related  to  the  blade  element.  (This  is  indicated  in  Figure  1.)  Note  that 
as  detailed  in  tne  input  data  description,  if  INAST  is  negative,  the  PLOAD2 
cards  are  not  punched. 


Precision  plots  are  produced  if  IFPLOT  = 1,  2,  3,  or  4 as  described 
under  the  definition  of  IFPLOT  given  previously. 
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2. 


AERODYNAMIC  SECTION 


a.  Regular  Printed  Output 

The  input  data  are  first  printed  out  in  its  entirety,  and  the  results 
for  each  running  point  follow.  The  output  is  generally  self-explanatory  and 
definitions  are  given  here  for  some  derived  quantities.  Tabular  output  is 
generally  not  started  on  a page  unless  it  can  be  completed  on  the  same  page, 
according  to  the  maximum  number  of  lines  permitted  by  the  input  variable 
NPAGE. 


The  results  of  each  running  point  are  given  under  a heading 
giving  the  running  point  number.  .Any  diagnostics  generated  during  the 
calculation  will  appear  first  under  the  heading.  (Diagnostics  are  described 
in  the  following  section.)  Then,  a station-by- station  print  out  follows  for 
each  station  through  to  the  last  station,  or  to  the  station  where  the  calcu- 
lation failed,  if  this  occurred.  One  or  more  diagnostics  will  indicate  the 
reason  for  the  failure,  in  this  event.  Included  in  the  meshpoint  coordinate 
data  is  the  distance  along  the  computing  station  from  the  hub  to  the  inter- 
ception of  the  streamline  with  the  station  (L),  and  the  station  lean  angle 
(GAMA).  Where  the  radius  of  curvature  of  a streamline  is  shown  as  zero, 
the  streamline  has  no  curvature.  The  whirl  angle  is  defined  by 


tan  o( 


Vo 

v* 


(i) 


For  stations  within  a blade,  or  at  a blade  trailing  edge,  a rel- 
ative total  pressure  loss  coefficient  is  shown.  The  loss  of  relative  total 
pressure  is  computed  from  the  station  defined  by  the  input  variable  NL1.  If 
a loss  coefficient  was  used  in  the  input  for  the  station  (NLOSS  = 1 or  4,  or 
NWORK  = 7),  the  input  variable  NL2  defines  the  station  where  the  norn.alizing 
relative  dynamic  head  is  taken;  otherwise,  it  is  taken  at  the  station  defined 
by  NL1.  If  the  cascade  solidity  is  given  as  anything  but  zero,  it  is  used  in 
the  determination  of  diffusion  factors.  The  following  definition  is  used: 


Vir  Ve„  v,; 

V„  l<r  V, 


Inlet  conditions  (subscript  1)  are  taken  from  the  station 
defined  by  the  input  variable  NL1. 
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The  last  term  in  Equation  2 is  multiplied  by  -1  if  the  blade 
speed  is  greater  than  zero,  or  the  blade  speed  is  zero  and  the  preceding 
rotating  blade  row  has  negative  rotation.  This  is  necessary  because 
relative  whirl  angles  are  (generally)  negative  for  rotor  blades  and  for 
stator  blades  that  follow  a rotor  having  "negative"  wheel  speed.  Incidence 
and  deviation  angles  are  treated  in  the  same  way.  so  that  positive  and 
negative  values  have  their  conventional  significance  for  all  blades. 

If  annulus  wall  boundary  layer  computations  were  made  (NBL 
= 1),  details  are  shown  foi  each  station.  Then,  an  overall  result  is  given, 
including  a statement  of  the  number  of  passes  that  have  been  performed  and 
whether  the  calculation  is  converged,  unconverged,  or  failed.  When  the 
calculation  is  unconverged,  the  number  of  mesh  points  where  the  meridional 
velocity  component  has  not  remained  constant  to  within  the  specified  toler- 
ance (TOLNCE)  on  the  last  two  passes  is  shown  as  IVFAIL.  Similarity,  the 
number  of  streamtubes,  defined  by  the  hub  and  each  streamline  in  turn, 
where  the  fraction  of  the  flow  is  not  within  the  same  tolerance  of  the  target 
value  is  shown  as  IFFAIL.  If  these  numbers  are  small,  say  less  than  10% 
of  the  maximum  possible  values,  the  results  may  generally  be  used.  Other- 
wise, the  computation  should  be  rerun,  either  for  a greater  number  of 
passes,  or  with  modified  relaxation  factor  constants.  If  the  punch- 
streamline  option  (NPUNCH  =1)  was  used,  the  streamline  pattern  generated 
in  the  first  run  may  be  input  into  the  second  run  (with  NREAD  =1),  restarting 
the  calculation  essentially  where  it  was  terminated  and  thus  effecting  a 
saving  in  computer  time.  The  default  option  relaxation  constants  will  gener- 
ally be  satisfactory  but  may  need  modification  for  some  cases.  If  insufficient 
damping  is  specified  by  the  constants,  the  streamlines  generated  will  tend 
to  oscillate  and  this  may  be  detected  by  observing  a relatively  small  radius 
of  curvature  for  the  mid-passage  streamline  that  also  changes  sign  from  one 
station  to  the  next.  This  may  be  corrected  by  rerunning  the  problem  (from 
scratch)  with  a lower  value  input  for  RCONST,  say,  of  4.  0 instead  of  6.  0.  When 
the  damping  is  excessive,  the  velocities  will  tend  to  remain  constant  while 
the  streamlines  will  not  adjust  rapidly  to  the  correct  locations.  This  will  be 
indicated  by  a small  IVFAIL  and  a relatively  large  IFFAIL.  For  optimum 
program  performance,  RCONST  should  be  increased,  and  the  stream’ine 
pattern  generated  thus  far  could  be  used  as  a starting  point.  The  second 
constant  XMMAX  (the  maximum  value  of  the  square  of  Mach  number  used 
in  the  relaxation  factor)  is  incorporated  so  tha*  high  subsonic  or  super- 
sonic cases  the  damping  does  not  decrease  i . ceptably.  The  default 
value  of  0.6  may  be  too  low  for  rapid  program  convergence  in  some  such 
cases. 


If  the  punching  of  blade  pressure  load  data  for  the  NASTRAN 
program  is  specified  (by  the  input  variable  NOUT3),a  self-explanatory 
printout  is  also  made.  The  blade  clement  numbering  scheme  is  the  same  as 
that  incorporated  into  both  blading  sections  of  the  program,  and  illustrated 
in  Figure  '. 
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If  the  loss  coefficient  re-estimation  routine  has  been  used  for 
any  bladerow(s)  (NEVAL  ^ 0),  a printout  summarizing  the  computations 
made  will  follow.  A heading  indicating  whether  the  re-estimation  was 
incorporated  into  the  overall  iterative  procedure  or  whether  it  was  merely 
made  "after  the  event"  is  first  printed.  Then  follows  a self-explanatory 
tabulation  of  various  quantities  involved  in  the  redetermination  of  the  loss 
coefficient  on  each  streamline. 

b.  Diagnostic  Printed  Output 

The  various  diagnostic  messages  that  may  be  produced  by  the 
aerodynamic  section  of  the  program  are  all  shown.  Where  a computed  value 
will  occur,  "x"  is  shown  here. 

JOB  STOPPED  - TOO  MUCH  INPUT  DATA 

The  above  message  will  occur  if  the  sun.  of  NSPEC  or  NDATA 
or  NDEL  for  all  stations  is  above  the  permitted  limit.  Execution  ceases. 

STATIC  ENTHALPY  BELOW  LIMIT  AT  xxx.  xxxxxExxx 

The  output  routine  (subroutine  UD0311)  calculates  static 
enthalpy  at  each  meshpoint  when  computing  the  various  output  parameters  and 
this  message  will  occur  if  a value  below  the  limit  (HMIN)  occurs.  The  limiting 
value  will  be  used,  and  the  results  printed  become  correspondingly  arbitrary. 
HMIN  is  set  in  the  Program  UD03AR  and  should  be  maintained  at  some 
positive  value  well  below  anv  value  that  will  be  validly  encountered  in 
calculation. 


PASSxxx  STA TIONxxx  STREAMLINExxx  PRANDTL -MEYER 
FUNCTION  NOT  CONVERGED  - USE  INLET  MACH  NO 

The  loss  coefficient  re-estimation  procedure  involves  iteratively 
solving  for  the  Mach  number  in  the  Prandtl-Meyer  function.  If  the  calculation 
does  not  converge  in  20  attempts,  the  above  message  is  printed,  and  as 
indicated,  the  Mach  number  following  the  expansion  (or  compression)  is 
assumed  to  equal  the  inlet  value.  (The  routine  only  prints  output  following 
the  completion  of  all  computations  and  printing  of  the  station -by-station 
output  data. ) 

PASSxxx  STA  TIONxxx  ITERA  TIONxxx  STREAMLINExxx 
MERIDIONAL  VELOCITY  UNCONVERGED  VM  = xx.  xxxxxxExx 
VM(OLD)  = xx.  xxxxxxExx 
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For  "analysis"  cases,  that  is  at  stations  where  relative  flow 
angle  is  specified,  the  calculation  of  meridional  velocity  proceeds 
iteratively  at  each  meshpoint  from  the  mid-streamline  to  the  case  and  then 
to  the  hub.  The  variable  LPMAX  (set  to  10  in  Subroutines  UD0308  and 
UD03Z6)  limits  the  maximum  number  of  iterations  that  may  be  made  at  a 
streamline  without  the  velocity  being  converged  before  the  calculation 
proceeds  to  the  next  streamline.  The  above  message  will  occur  if  all  iter- 
ations are  used  without  achieving  convergence,  and  the  pass  number  is 
greater  than  NFORCE.  Convergence  is  here  defined  as  occurring  when  the 
velocity  repeats  to  within  TOLNCE/5.0,  applied  nondimensionally.  No 
other  program  action  occurs. 

PASSxxx  STATIONxxx  MOMENTUM  AND/OR  CONTINUITY 
UNCONVERGED  W/W  SPEC  = xx.  xxxxx  VM/VM  (OLD)  HUB  = 
xx.  xxxxxMID=xx.  xxxxx  TIP  = xx.  xxxxx 

If, following  completion  of  all  ITMAX  iterations  permitted 
for  the  flow  rate  or  meridional  velocity,  the  simultaneous  solution  of  the 
momentum  and  continuity  equations  profile  is  unconverged,  and  the  pass 
number  is  greater  than  NFORCE,  the  above  message  occurs.  Here  con- 
verged means  that  the  flow  rate  equals  the  specified  value,  and  the 
j meridional  velocity  repeats,  to  within  TOLNCE/5.0,  applied  nondimension- 

; ally.  If  loss  coefficient  re-estimation  is  specified  (NEVAL>  0),  an 

j additional  iteration  is  involved,  and  the  tolerance  is  halved.  No  further 

I program  action  occurs. 

j PASSxxx  STATIONxxx  VM  PROFILE  NOT  CONVERGED  WITH 

LOSS  RECALC  VM  NEW/ VM  PREV  HUB  = xx.  xxxxxx  MID  = 
xx.  xxxxxx  CASE  = xx.  xxxxxx 

When  loss  re-estimation  is  specified  (NEVAL>  0),  up  to 
NLITER  solutions  to  the  momentum  and  continuity  equations  are  completed, 
each  with  a revised  loss  coefficier.  variation.  If,  when  the  pass  number  is 
greater  than  NFORCE,  the  velocity  profile  is  not  converged  after  the 
NLITER  cycles  of  calculation  have  been  performed,  the  above  message  is 
issued.  For  convergence,  the  n.eridional  velocities  niust  repeat  to  within 
TOLNCE/5.0,  applied  nondimensionally.  No  further  program  action  occurs. 

A further  check  on  the  convergence  of  this  procedure  is  to 
compare  the  loss  coefficients  used  on  the  final  pass  of  calculation,  and  thus 
shown  in  the  station-by-station  results,  with  those  shown  in  the  output  from 
the  loss  coefficient  re -estimation  routine,  which  are  computed  from  the 
final  velocities,  etc. 
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PASSxxx  STATIONxxx  ITERATIONxxx  STREAMTUBExxx  STATIC 
ENTHALPY  BELOW  LIMIT  IN  MOMENTUM  EQUATION  AT 
xxx.  xxxxxExxx 

The  static  enthalpy  is  calculated  (to  find  the  static  temperature) 
during  computation  of  the  "design"  case  momentum  equation,  that  is,  when 
whirl  velocity  is  specified.  If  a value  lower  than  HMIN  (see  discussion  of 
second  diagnostic  message)  is  produced,  the  limiting  value  is  Inserted.  If 
this  occurs  when  IPASS  > NFORCE,  the  above  message  is  printed.  If  this 
occurs  on  the  final  iteration,  the  calculation  is  deemed  to  have  failed, 
calculation  ceases,  and  results  are  printed  out  through  to  this  station. 

PASSxxx  STATIONxxx  ITERATIONxxx  STREAMTUBExxx  LOOPxxx 
STATIC  H IN  MOMENTUM  EQUN.  BELOW  LIMIT  AT  xxx.  xxxxxExxx 

This  corresponds  to  the  previous  message,  but  for  the 
"analysis"  case.  For  failure,  it  must  occur  on  the  final  iteration  and  loop. 

PASSxxx  STATIONxxx  ITERATIONxxx  STREAMTUBExxx 
MERIDIONAL  MACH  NUMBER  ABOVE  LIMIT  AT  xxx.  xxxxxExx 

When  Subroutine  UD0308  is  selected  (NEQN  = 0 or  1),  the 
meridional  Mach  number  is  calculated  during  computation  of  the  design 
momentum  equation,  and  a maximum  value  of  0.  99  is  permitted.  If  a 
higher  value  is  calculated,  the  limiting  value  is  inserted.  If  this  occurs 
when  IPASS  > NFORCE,  the  above  message  is  printed.  If  this  occurs  on 
the  final  iteration,  the  calculation  is  deemed  to  have  failed,  calculation 
ceases,  and  results  are  printed  through  to  this  station. 

PASSxxx  STATIONxxx  ITERAT'.ONxxx  STREAMTUBExxx  LOOPxxx 
MERIDIONAL  MACH  NUMBER  ABOVE  LIMIT  AT  xxx.  xxxxxExxx 

This  corresponds  to  the  previous  message,  but  for  the  "analysis" 
case.  For  failure,  it  must  occur  at  the  final  iteration  and  loop. 

PASSxxx  STATIONxxx  ITERATIONxxx  STREAMTUBExxx 
MOMENTUM  EQUATION  EXPONENT  ABOVE  LIMIT  AT  xxx.  xxxxxExxx 

An  exponentiation  is  performed  during  the  computation  of  the 
design  case  momentum  equation,  and  the  maximum  value  of  the  exponent  is 
limited  to  88.0.  If  this  substitution  is  required  when  IPASS  > NFORCE, 
the  above  message  is  printed.  If  it  occurs  on  the  final  iteration,  the  calcu- 
lation is  deemed  to  have  failed,  calculation  ceases,  and  results  are  printed 
through  to  this  station. 
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PASSxxx  STATIC) Nxxx  ITERATIONSxxx  STREAMLINExxx 
(MERIDIONAL  VELOCITY)  SQUARED  BELOW  LIMIT  AT 
xxx.  xxxxxExxx. 

If  a meridional  velocity,  squared,  of  less  than  1.  0 is  calcu- 
lated during  computation  of  the  design-case  momentum  equation,  this  limit 
is  imposed.  If  this  occurs  when  IPASS ENFORCE,  the  above  message  is 
printed.  If  this  occurs  on  the  final  iteration,  the  calculation  is  deemed  to 
have  failed,  calculation  ceases,  and  results  are  printed  out  through  to  this 
station. 

PASSxxx  STATIONxxx  ITERATIONxxx  STREAMLINExxx  LOOPxxx 
(MERIDIONAL  VELOCITY)  SQUARED  BELOW  LIMIT  AT 
xxx.  xxxxxExxx. 

This  corresponds  to  the  previous  message,  but  for  the 
"analysis"  case.  For  failure,  it  must  occur  on  the  last  iteration  and  loop. 

PASSxxx  STATIONxxx  ITERATIONxxx  STREAMTUBExxx 
STATIC  ENTHALPY  BELOW  LIMIT  IN  CONTINUITY  EQUATION 
AT  xxx.  xxxxxExxx. 

The  static  enthalpy  is  calculated  during  computation  of  the 
continuity  equation.  If  a value  lower  than  HMIN  (see  discussion  of  second 
diagnostic  message)  is  produced,  the  limiting  value  is  imposed.  If  this 
occurs  when  IPASS>NFGRCE,  the  above  message  is  printed.  If  this 
occurs  on  the  final  iteration,  the  calculation  is  deemed  to  have  failed, 
calculation  ceases,  and  results  are  printed  out  through  to  this  station. 

PASSxxx  STATIONxx  ITERATIONxxx  STREAMLINExxx 
MERIDIONAL  VELOCITY  BELOW  LIMIT  IN  CONTINUITY  AT 
xxx.  xxxxxExxx. 

If  a meridional  velocity  of  less  than  1. 0 is  calculated  when  the 
velocity  profile  is  incremented  by  the  amount  estimated  to  be  required  to 
satisfy  continuity,  this  limit  is  imposed.  If  this  occurs  when  IPASS  > 
NFOPCE,  the  above  message  is  printed.  If  this  occurs  on  the  final  iteration, 
the  calculation  is  deemed  to  have  failed,  calculation  ceases,  and  results  are 
printed  through  to  this  station. 

PASSxxx  STATIONxxx  ITERATIONxxx  OTHER  CONTINUITY 
EQUATION  BRANCH  REQUIRED 

If  when  IPASS>NFORCE.  a velocity  profile  is  produced  that 
corresponds  to  a subsonic  solution  to  the  continuity  equation  when  a super- 
sonic solution  is  required,  or  vice  versa,  the  above  message  is  printed.  If 
this  occurs  on  the  final  iteration,  failure  is  deemed  .o  have  occurred,  calcu- 
lation ceases,  and  results  are  printed  out  through  to  this  station. 
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PASSxxx  STATIONxxx  ITERATIONxxx  STREAMLINExxx 
MERIDIONAL  VELOCITY  GREATER  THAN  TWICE  MID  VALUE 

During  integration  of  the  "design"  momentum  equations, 
no  meridional  velocity  is  permitted  to  be  greater  than  twice  the  value  on 
the  mid-streamline.  If  this  occurs  when  IPASS>NFORCE,  the  above 
iriessage  is  printed.  If  this  occurs  on  the  final  iteration,  the  calculation 
is  deemed  to  have  failed,  calculation  ceases,  and  results  are  printed 
through  to  this  station.  In  the  event  that  this  limit  interferes  with  a valid 
velocity  profile,  the  constants  that  appear  on  cards  $08$.  272,  $08$.  279, 

$26$.  229,  and  $26$.  236  may  be  modified  accordingly.  Note  that  as  the 
calculation  is  at  this  point  working  with  the  square  of  the  meridional 
velocity,  the  constant  for  a limit  of  2.  0 tin.es  the  mid-streamline  value, 
for  instance,  appears  as  4.0. 

PASSxxx  STATIONxxx  ITERATIONxxx  STREAMLINExxx 
LOOPxxx  MERIDIONAL  VELOCITY  ABOVE  LIMIT  xxxxxExx 
LIMIT  - xxxxxExx. 

During  integration  of  the  "analysis"  momentum  equations,  no 
meridional  velocity  is  permitted  ?.->  be  greater  than  three  times  the  value 
on  the  mid-streamline.  If  this  o c .rs  when  IPASS>NFORCE,  the  above- 
message  is  printed.  If  this  occurs  on  the  final  loop  of  the  final  iteration, 
the  calculation  is  deemed  to  have  failed,  calculation  ceases,  and  results 
are  printed  through  to  this  station.  In  the  event  that  the  limit  interferes 
with  a valid  velocity  profile,  the  constants  that  appear  on  cards 
$08$.  398,  $08$.  409,  $26$.  323,  $26$.  334,  and  $26$.  329  may  be  modified 
accordingly.  In  each  case  except  that  of  the  last  card  noted,  the  program 
is  working  with  meridional  velocity  squared,  so  that  a lin.it  of,  for  instance, 
3.0  times  the  mid- streamline  value  appears  as  9.  0. 

PAiSxxx  STATIONxxx  STREAMLINExxx  LIMITING  MERIDIONAL 
VELOCITY  SQUARED  = xxxxxExx. 

In  the  Subroutine  UD0308  (NEQN=  0 or  1),  a maximum 
permissable  meridional  velocity  (equal  to  the  speed  of  sound)  is  established 
for  each  streamline  at  the  beginning  of  each  pass.  The  calculation  yields 
the  square  of  the  velocity,  and  if  a value  of  less  than  1. 0 is  obtained,  a 
value  of  6250000.0  is  superimposed  (which  corresponds  to  a meridional 
velocity  of  2500.0).  If  this  occurs  when  IPASS  >NFORCE,  the  above  message 
is  printed,  and  the  calculation  is  deemed  to  have  failed.  Calculation  ceases 
after  the  station  computations  are  made,  and  results  are  printed  through 
to  this  station. 
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PASSxxx  STATIONxxx  ITERATIONxxx  STREAMLINExxx 
MERIDIONAL  VELOCITY  ABOVE  SOUND  S-^EED  VM  = 
xxxx. xx  A = xxxx. xx. 


In  Subroutine  UD0308  (NEQN  = 0 or  1),  no  meridional  velocity 
is  permitted  to  be  larger  than  the  speed  of  sound.  The  above  message  will 
occur  if  this  limit  is  violated  during  integration  of  the  "design"  n. omentum 
when  IPASS  > NFORCE.  If  the  limit  is  violated  at  any  point  when  IPASS> 
NFORCE  and  on  the  last  permitted  iteration  (last  permitted  loop  also  in 
the  case  of  the  "analysis"  momentum  equation),  the  calculation  is  deemed 
to  have  failed.  Calculation  ceases,  and  the  results  are  printed  through  to 
this  station. 


MIXING  CALCULATION  FAILURE  NO.  n 

The  above  message  occurs  when  flow  mixing  calculations  are 
specified,  and  the  computation  fails.  The  overall  calculation  is  halted,  and 
results  are  printed  through  to  the  station  that  is  the  upstream  boundary  for 
the  mixing  interval  in  which  the  failure  occurred.  The  integer  n takes  on 
different  values  to  indicate  the  specific  problems  as  follows. 

n = 1 In  solving  for  the  static  pressure  distribution  at  the  upstream 

boundary  of  each  mixing  s*ep,  the  average  static  enthalpy  is 
determined  in  each  streamtube  (defined  by  an  adjacent  pair 
of  streamlines).  This  failure  indicates  that  a value  less  than 
HMIN  was  determined. 

n = Z Calculation  of  the  static  pressure  distribution  at  the  upstream 

boundary  of  the  mixing  step  is  iterative.  This  failure  indicates 
that  the  procedure  was  not  converged  after  10  iterations. 

n = 3 The  static  enthalpy  on  each  streamline  at  the  mixing  step 

upstream  boundary  is  determined  from  the  static  pressure 
and  entropy  there.  This  failure  indicates  that  a value  less 
than  HMIN  was  determined. 

n = 4 The  axial  velocity  distribution  at  the  mixing  step  upstream 

boundary  is  determined  from  the  total  enthalpy,  static  enthalpy, 
and  tangential  velocity  distributions.  This  failure  indicates 
that  a value  less  than  VMIN  was  determined. 

n = 5 In  solving  for  the  otatic  pressure  distribution  at  the  downstream 

boundary  of  each  mixing  step,  the  average  static  enthalpy  is 
determined  in  each  streamtube  (defined  by  an  adjacent  pair 
of  streamlines).  This  failure  indicates  that  a value  less  than 
HMIN  was  determined. 
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n = 6 Calculation  of  the  static  pressure  distribution  at  the  down- 

stream boundary  of  the  mixing  step  is  iterative.  This  failure 
indicates  that  the  procedure  was  not  converged  after  10 
iterations. 

n = 7 The  static  enthalpy  distribution  at  the  mixing  step  downstream 

boundary  is  found  from  the  total  enthalpy,  axial  velocity,  and 
tangential  velocity  distributions..  This  failure  indicates  that 
a value  less  than  HMIN  was  determined. 


In  order  to  satisfy  continuity,  the  static  pressure  level  at  the 
mixing  step  downstream  boundary  is  iteratively  determined. 

This  failure  indicates  that  after  15  attempts,  the  procedure 
was  unconverged. 

Punched  Card  Output 

Four  output  options  may  result  in  punched  cards  being  produced 
by  the  aerodynamic  section  of  the  program.  Use  of  the  input  item  NOUT3 
gives  "I*LOAD2  - Cards"  punched  in  a format  compatible  with  the  NASTRAN 
stress  analysis  program.  For  the  purposes  of  stress  analysis,  the  blade 
is  taker,  to  be  composed  of  a number  of  triangular  elements.  Two  such 
elements  are  formed  by  the  quadrilateral  defined  by  two  adjacent  stream  lines 
and  two  adjacent  computing  stations.  The  way  that  each  quadrilateral  is 
divided  into  twi  triangles,  and  the  element  numbering  scheme  that  is  used, 
are  illustrated  in  Figure  1.  The  pressure  difference  for  each  element  is 
given  by  an  average  of  either  three  or  four  values  at  surrounding  mesh- 
points.  The  pressure  difference  at  each  meshpoint  is  computed  from  the 
equation 


Vm  jL(fVe)^ 


(3) 


and  as  follows.  At  the  blade  leading  edge  (NOUT3  = 1 or  3),  a forward 
difference  is  used  to  determine  the  meridional  gradients.  At  the  blade 
trailing  edge  (NOUT3  = 2 or  3),  the  pressure  difference  is  taken  to  be  zero. 
At  stations  with  the  bladerow  (NOUT3  = 0,  following  a leading  edge),  mean 
central  differences  are  used  to  determine  the  meridional  gradients.  When 
the  input  item  NBLADE  is  positive  (or  zero)  for  a particular  blade  axial 
segment,  then  three-point  averaging  is  used.  For  instance,  for  element 
number  1 in  Figure  1,  pressure  differences  at  grid  points  1,  6,  and  7 would 
be  used.  If  NBLADE  is  negative,  four-point  averaging  is  used.  For 
instance  for  element  number  1 in  Figure  1,  pressure  differences  at  grid 
points  1,  2,  6,  and  7 would  be  used.  The  same  average  would  also  apply  to 
element  number  2. 
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If  NPUNCH  = 1,  cards  will  be  punched  that  may  be  used  as 
input  when  the  NREAD  = 1 option  is  used.  These  comprise  one  or  more  cards 
defining  the  DELF  values,  and  then  (NSTRMS  times  NSTNS)  cards  giving  the 
streamline  coordinates.  These  may  be  incorporated  directly  into  a sub- 
sequent input  data  deck.  Note  that  one  set  will  be  punched  for  every  running 
point  that  does  not  fail. 

For  any  station  for  which  NOUT1  = 1,  NSTRMS  cards  giving 
the  streamline  radius  (R)  will  be  punched  that  then  may  be  used  in  a 
subsequent  input  data  deck  for  the  analytic  meanline  blade  section.  A set 
will  be  punched  for  each  running  point  that  does  not  fail,  immediately 
following  the  NPUNCH  = i cards  i.  r the  running  point,  if  this  option  is  also 
specified. 

The  input  item  NOUT2  controls  the  final  possible  card  output. 

If  NARBIT  ^ 0,  this  output  is  stored  within  the  computer  system  on  file 
LOG6;  otherwise,  it  will  be  placed  on  the  punch  file.  For  any  station  where 
NOUT2  £ 0,  data  are  output  that  give  KPTS,  XSTA,  and  RSTA  (that  define 
the  station  geometry),  and  R and  AIRANG  (that  define  the  stream  surface 
radii  and  relative  air  angles)  for  possible  subsequent  use  in  the  arbitary 
meanline  section  of  the  progrm.  (Note  that  BLAFOR  is  left  blank  (zero); 
the  first  punch  option  described  above  will  give  PLOAD2  -cards  for  sub- 
sequent use  with  the  NASTRAN  program  if  desired.) 

d.  Precision  Plot  Output 

Plots  Gf  the  static  pressure  distribution  and  final  streamline 
pattern  may  be  made  for  each  running  point  according  to  the  input  variable 
NPLOT.  The  program  assumes  the  availability  of  CALCOMP  software  to 
generate  instructions  for  a CALCOMP  plotter.  Generally,  an  11 -inch 
plotter  would  be  used,  and  the  input  variables  XSCALE,  PSCALE,  RLOW, 
and  PLOW  should  be  set  accordingly. 

3.  ARBITRARY  MEANLINE  SECTION 

Printed  output  from  this  section  of  the  program  may  be  considered  to 
consist  of  four  sections;  a printout  of  the  input  data;  details  of  the  blade 
sections  on  each  stream  surface;  a listing  of  quantities  required  for  aero- 
dynamic analysis;  and  details  of  the  manufacturing  sections  determined  on 
the  constant-z  planes.  These  are  briefly  described  below. 

The  input  data  printout  includes  all  quantities  read  in,  and  is  self- 
explanatory. 
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Details  of  the  streaun surface  blade  sections  are  printed  if  IPRINT  = 

0 or  1.  The  iterative  determination  of  a consistent  deviation  angle  and 
cascade  solidity  is  recorded  by  the  successive  values  of  deviation  and 
solidity  that  are  produced.  Next  the  camberline  that  meets  the  angle/ 
meridional  coordinate  specification  thus  created  is  detailed,  camberline 
coordinates,  first  and  second  derivatives,  and  radius  of  curvature  being 
shown.  Also  shown  are  the  camberline  angles  and  meridional  coordinates 
for  which  the  camberline  was  computed.  Then  follow  details  of  the  blade 
section  in  "normalized"  form.  The  blade  section  geometry  is  given  for 
the  section  specified,  except  that  the  meridional  projection  of  the  chord  is 
unity.  For  this  section  of  the  output,  the  coordinate  origin  is  the  blade 
leading  edge.  The  following  quantities  are  given:  blade  chord;  stagger  angle; 
camber  angle;  section  area;  location  of  centroid  of  the  section;  second  n.oments 
of  area  of  the  section  about  the  centroid;  orientation  of  the  principal  axes; 
and  the  principal  second  moments  of  area  of  the  section  about  the  centroid. 

Then  are  listed  the  coordinates  of  the  camber  line,  the  camber  line  angle, 
the  section  thickness,  and  the  coordinates  of  the  blade  surfaces.  A line- 
printer  plot  of  the  normalized  section  follows.  The  scales  for  the  plot  are 
arranged  so  that  the  section  just  fills  the  page,  so  that  the  scales  will 
generally  differ  from  one  plot  to  another.  "Dimensional"  details  of  the  blade 
section  are  given  next.  The  normalized  data  given  previously  are  scaled 
to  give  the  proper  blade  section.  For  this  section  of  the  output,  the 
coordinates  are  with  respect  to  the  blade  stacking  axis.  The  following  quan- 
tities are  given:  blade  chord;  radius  and  location  of  center  of  the  leading 
edge;  section  area;  the  second  moments  of  area  of  the  section  about  the 
centroid.  The  cooidinates  of  points  on  the  blade  surfaces  are  then  listed, 
followed  by  the  coordinates  of  31  points  distributed  at  six  degree  intervals 
around  the  leading  edge.  Finally,  the  coordinates  of  the  blade  surfaces  and 
points  around  the  leading  edge  are  shown  in  Cartesian  form. 


| The  quantities  required  for  aerodynamic  analysis  are  printed  at  all 

E computing  stations  specified  by  the  IFANGS  parameter.  The  radius,  blade 

[section  angle,  blade  lean  angle,  blade  olockage,  and  relative  angular  location 
at  the  camber  line  are  printed  at  each  streamsurface  intersection  with  the 
particular  computing  section. 


Details  of  the  manufacturing  sections  are  printed  if  IPRINT  = 0 or  2. 

At  each  value  of  z specified  by  ZINNER,  ZOUTER,  and  NZ,  section 
properties  and  coordinates  are  given.  The  origin  for  the  coordinates  is 
the  blade  stacking  axis.  The  following  quantities  are  given:  section  area; 
the  location  of  the  centroid  of  the  section;  the  second  moments  of  area  of 
the  section  about  the  centroid;  the  principal  second  moments  of  area  of  the 
section  about  the  centroid;  the  orientation  of  the  principal  axes;  and  the 
section  torsional  constant.  Then  the  coordinates  of  points  on  the  blade 
section  surfaces  are  listed,  followed  by  31  points  around  the  leading  edge. 
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Precision  plots  are  produced  if  IFPLOT  =1,  2,  3,  or  4 as  described 
under  the  definition  of  IFPLOT  gi'-e>-  previously. 


If  IPUNCH  = 1,  the  program  punches  the  quantities  required  for  aero- 
dynamic analysis,  together  with  identifying  indices  denoting  station  number 
and  streamsurface  number,  on  cards  in  the  following  format;  5 fields  each 
of  12  locations  for  the  quantities  themselves,  followed  by  2 fields  each  of  3 
locations  for  the  indices. 


If  INAST  ^ 0,  punched  cards  are  output  for  use  with  the  NASTRAN 
stress  analysis  program.  See  the  description  given  earlier  for  the  corres- 
ponding output  from  the  analytic  meanline  section  of  the  program  for 
details. 
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SECTION  V 

THEORY  OF  AERODYNAMIC  ANALYSIS 


1.  MOMENTUM  EQUATION 

The  meridional  velocity  distribution  along  each  computing  station  is 
found  by  integrating  a momentum  equation.  (The  continuity  equation  supplies 
the  necessary  constant  of  integration.)  Two  different,  basic  momentum 
equations  are  included  in  the  program,  and  for  each  of  them  there  are  forms 
that  correspond  to  flow  in  the  bladed  and  nonbladed  regions  of  the  compressors 
Thus  there  are  four  momentum  equations  in  all,  and  selection  of  the  one  (or 
two)  appropriate  to  an  application  is  controlled  by  the  input  variable  NEQN. 
Each  of  the  four  momentum  equations  may  be  written  (and  solved)  in  two 
ways,  depending  upon  whether  the  calculation  is  of  the  "analysis"  or  "design" 
type.  In  the  analysis  case,  the  relative  flow  angle  ( |?)  is  fixed,  whereas  in 
the  design  case,  it  is  the  absolute  angular  momentum  (t*Vs  ) that  is  fixed. 

One  application  of  the  program  will  usually  involve  two  forms  of  momentum 
equation  (for  the  bladed  and  nonbladed  regions  of  the  compressor),  and  some 
analysis-type  and  some  design-type  calculations,  so-called,  regardless  of 
whether  the  computation  is  actually  a compressor  design  or  analysis 
exercise. 


The  starting  point  for  the  development  of  the  various  equations  is 
presented  by  Wennerstrom  in  Reference  8,  Equation  14,  and  reproduced 
here.  (Our  sign  convention  for  the  station  lean  angle  ( y ) is  the  reverse 
of  his. ) 
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mw  is  a force  acting  in  the  meridional  plane  and  in  the  (projected)  stream- 
line direction.  Fn  is  a force  also  in  the  meridional  plane  and  normal  to 
Fm  . Equation  4 is  our  first  basic  momentum  equation.  In  Reference  8, 
Wennerstrom  defines  the  forces  Fm  and  Fn  in  the  context  of  the  bladed 
regions  of  a turbomachine.  His  Equations  3 and  19  define  F„  thus: 
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Fn  is  defined  by  Wennerstrom' s Equations  25  and  27  as  follows: 

F„  = (s«ps«i  7 ^fr  H-  s«P  w*  3T*  j|  ) x 

S**»  C.  s«c  (M)  +-  S«*  p tos  £ Tan  +-^ ) ^ 


(6) 


The  first  basic  equation,  adapted  for  the  bladed  regions  of  the  compressor, 
is  given  by  combining  Equations  4,  5,  and  6.  If  the  input  variable  NEQN  = 
2,  this  equation  is  used  by  the  computer  program  at  computing  stations  that 
are  at  the  edges  of,  or  within,  blade  rows.  (At  leading  edges,  if  the  blade 
geometry  is  not  given,  the  equation  adapted  for  nonbladed  regions  is  used.) 

_ The  second  equation  that  we  obtain  is  given  by  considering  the  forces 
•Vn  and  Fn  for  nonbladed  regions  of  the  compressor.  Our  result  differs 
slightly  from  that  given  by  Wennerstrom  in  Reference  8,  but  he  has  agreed 
that  the  treatment  here  is  correct.  In  the  nonbladed  regions,  there  can  be 
no  force  normal  to  the  streamline  direction,  nor  can  there  by  any  change 
in  angular  momentum  along  a streamline.  Then  the  forces  become  as 
follows: 


(7) 


(8) 


The  first  basic  equation,  adapted  for  the  nonbladed  regions  of  the  compressor, 
is  given  by  combining  Equations  4,  7,  and  8.  If  the  input  variable  NEQN  = 2, 
this  equation  is  used  by  the  computer  program  at  computing  stations  that 
are  neither  at  the  edges  of,  nor  within,  bladerows.  If  the  input  variable 
NEQN  = 3,  this  equation  is  used  at  all  computing  stations. 


55 


For  the  second  pair  of  equations,  the  basic  momentum  equation 

(Equation  4)  is  rewritten.  The  gradient  *s  eliminated  through  use  of 

the  continuity  equation.  The  rules  of  partial  differentiation  together  with 
the  geometric  relationship 
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The  tern,  ^ £ must  now  be  eliminated,  and  the  continuity  equation  is  used 
for  this.  The  continuity  equation  may  be  written  (for  axisymmetric  flow) 


This  may  be  expanded  to  give 
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Assuming  the  fluid  to  be  a perfect  gas,  the  meridional  density  gradient 
may  be  replaced  using 
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Note  that  the  possibility  of  streamwise  gradients  of  entropy  and  angular 
momentum  is  included.  Equations  10,  12,  and  13  are  combined  to  give 
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Equation  14  is  used  to  eliminate  ''am  from  Equation  4 to  give  the  second 
basic  equation.  For  the  bladed  regions.  Equations  5 and  6 again  give  Fm 
and  Fh  , and  the  second  basic  equation  adapted  for  the  bladed  regions  of  the 
compressor  is  obtained  thus: 

Vm  4 -V-1  ^Wecesfr 

ST  “ rc  cosW+vju-nN‘)  r dH 


+ | (</+^)  (tos  (3  + p-j£»  ) — 


If  the  input  variable  NEQN  = 0,  this  equation  is  used  by  the  computer 
program  at  computing  stations  that  are  at  the  edges  of,  or  within,  blade 
rows.  (At  leading  edges,  if  the  blade  geometry  is  not  given,  the  equation 
adapted  for  nonbladed  regions,  Equation  16,  is  used.) 

For  the  nonbladed  regions,  the  streamwise  gradient  of  angular 
momentum  in  Equation  14  becomes  zero,  and  Fr*  and  F«  are  given  by 
Equations  7 and  8.  Hence  the  second  basic  equation  adapted  for  the  non- 
bladed regions  is  obtained  thus: 
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V IF  t C®s  + 2r)(l  — H-^)  tr  ST 

i T c*r  a7+dS  .tjtdS  l-)l*0f-i)  v;mw)  ; 
+ 3 5I_3J'rf3l+34  B,  -pni3 rrrc- 

$ - V^j^){sA/0+H«UV)+^) 


(16) 


If  the  input  variable  NEQN  = 0,  this  equation  is  used  by  the  computer 
program  at  computing  stations  that  are  neither  at  the  edges  of,  nor  within, 
bladerows.  If  the  input  variable  NEQN  = 1,  this  equation  is  used  at  all 
computing  stations. 


The  remaining  manipulations  are  those  to  write  the  equations  in  forms 
convenient  for  design  or  analysis -type  calculations.  For  analysis  calculations, 
when  the  relative  flow  angle  ( ^)  is  fixed,  it  is  convenient  to  use  a coordinate 
system  rotating  at  blade  speed.  Then  the  rothalpy  (I)  is  constant  along 
streamlines  through  any  one  bladerow.  Replacing  W®  by  Vw  fa"  £ in 
Equation  4,  and  using  Equations  5 and  6 for  the  forces  Fm  and  rn  , the 
analysis  form  of  the  first  basic  equation  adapted  for  the  bladed  regions  is 
obtained  thus: 


s)  \i  \ n y , -rdl  nTfcK  ^ (fW 


The  analysis  form  of  Equation  4 for  nonbladed  regions  (which  will  only  be 
used  when  blade  force  terms  are  deliberately  omitted  by  setting  NEQN  to 
3),  is  given  by  replacing  W®  by  Vrtu  tftnp  , and  using  Equations  7 and  8 
for  the  forces  and  F«  . The  result  is; 


Vml  Up 

7 " r X 

C (18) 
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For  the  nonbladed  regions,  the  analysis  form  of  Equation  16  is: 

V*  f - v4* 

- 2t»V».tP,pl«,iV4  «|3  ^ ''  + * ,22) 

(This  equation  will  only  be  used  when  biaJe  force  terms  are  deliberately 
omitted  by  setting  NEQN  to  1.  ) 

The  design  form  of  Equation  15  (in  stationary  coordinates,  and  for  the 
bla.'^d  regions)  is; 

v.f -,y{f  -jjtir  *?  , % i 

-wMi - 4 v‘ 


+ : 


| ^y*)  djf  _ ^4)^(1  -H.1+  ^ ) 

i cos(f*-S)  <** 


Finally,  for  the  unbladed  regions,  where  Equation  16  pertains,  the  design 
form  is: 

\l  clV*  TdH  -rid?  Ve  cifrW  UcK  cJfL+t]* 

V-  n = 3T5J  _331cQ  ~ T ST 


(24) 


I - kx(Y-i^  Jtl  f 

l-hw‘  l-H4(  ‘WH)  ** 

-**(^s)(.^0+H<>)  + -tfe) 
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2. 


CONTINUITY  EQUATION 


The  continuity  equation  is  used  at  each  computing  station,  and  serves 
to  fix  the  constant  of  integration  in  the  momentum  equation.  Tne  equation  is: 

y = ^ dW 

(25) 

= ^ to  Vm  cos  (</+&)  2n r ( \ -£>) ctE. 


In  order  to  determine  on  which  branch  of  the  continuity  equation  a 
velocity  profile  lies,  and  to  be  able  to  estimate  by  how  much  to  change  the 
velocity  level  to  satisfy  continuity,  the  rate  of  change  of  flow  with  mid- 
radius meridional  velocity  is  required.  If  a perfect  gas  is  assumed,  then, 
for  the  design  case  (angular  momentum  fixed)  it  is  convenient  to  write: 


U)-  k)T 


(' 


v^v; 

23K,T 


(26) 


Using  this  in  Equation  25  and  differentiating  with  respect  to  midradius 
meridional  velocity,  the  result  is: 


dW  f Mm  !-H„l  jy 

J dVwmii  V m 


It  is  assumed  that: 

jVw  _ V* 


(28) 


Then  the  final  result  is; 
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(29) 


For  the  analysis  case,  the  relative  flow  angle  and  relative  total 
density  are  constant,  and  we  write: 

I 


U)  = 


V»  (l+Wp)  s M 

23T&T*  1 


(30) 


Using  this  in  Equation  25,  and  differentiating  with  respect  to  midradius 
meridional  velocity,  the  result  is: 


dw_  fai. 


(31) 


For  analysis  type  calculations,  it  is  often  important  to  evaluate  Equation 
31  accurately,  and  therefore  Equation  28  is  not  assumed.  Instead,  the 
meridional  velocity  at  any  point  on  a computing  station  is  represented  by: 


Yw(-f)  = Vmj+ 

Hence: 


(32) 


(33) 


is  derived  from  Equation  21,  22,  17,  or  18, 


d (dim) 

The  gradient  dwell' 


depending  upon  whether  the  input  variable  NEQN  = 0,  1,  2,  or  3,  respect- 
ively. In  the  cases  where  Equations  17  or  21  apply,  the  streamwise 
gradient  of  angular  momentum  is  neglected.  Thus  from  Equation  21,  we 
obtain:  _ 


(dVm]  _ cos'B 

dUJJj- 


3^7  - A 


V«‘ 


r cU 


When  Equation  22  applies,  the  result  is  the  same  except  that  the  coefficient 
of  in  Equation  34  is  replaced  by 


V*1 


From  Equation  17,  we  obtain 


£(f)=  ^ 

~ ^ ^ ^ ^1  ( ^ ( ^ P ~ to*»  1.  S**  {3  CCS  p)| 


(35) 


When  Equation  18  applies,  the  result  is  the  same  except  that  the  coefficient 
of  in  Equation  ^5  is  replaced  by: 


S'w»  ((f  + V ) ) 


So  for  analysis-type  calculations,  the  gradient  ' *s  g'ven 

by  Equations  31  and  33,  with  Equations  34  or  35  modified  as  appropriate. 


3. 


FLUID  PROPERTIES 


As  mentioned  previously,  the  fluid  properties  are  computed  in  the 
program  for  a perfect  gas.  The  following  basic  relationships  are  used: 


dh  = Ns  + Vpty 

)i  = Cft 

t>  = 

By  arbitrarily  defining  the  entropy  as  zero  for  a pressure  and 
temperature  of  1.0,  the  following  relationships  ma;  be  derived: 

k = Cp  «x|>  ( S/Cp  + fy’SCp  j’) 

S = Cf  Uje  (V<3>)  - ^ 


(36) 

(37) 

(38) 


(39) 


(40) 


(41) 


The  three  previous  equations  are  used  to  relate  enthalpies,  pressure, 
and  entropies.  The  following  two  relationships  are  also  required: 


(42) 


(43' 
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4.  RELATIVE  TOTAL  PRESSURE  LOSS  COEFFICIENTS 

Relative  total  pressure  loss  coefficients  may  be  used  in  the  input 
to  specify,  indirectly,  the  entropy  rise  in  the  flow,  and  are  calculated  for 
the  output  data.  When  the  dynamic  head  specified  to  normalize  the  loss 
coefficient  occurs  at  some  location  different  to  that  where  the  total  pressure 
is  sought,  the  following  definition  applies  : 

Pr-Pr-io(Pr,-|>4)  <«> 

This  gives  the  actual  relative  total  pressure  as  the  ideal  value  (com- 
puted from  the  station  defined  by  the  input  variable  NL1)  less  the  pressure 
loss  which  is  the  product  of  the  loss  coefficient  and  the  relative  dynamic 
head  (computed  at  the  station  specified  by  the  input  variable  NL2).  The 
isentropic  relative  total  pressure  P/  is  given  by  the  entropy  at  the  station 
defined  by  NL1  and  the  relative  total  enth  Hoy.  Basically,  the  same  defini- 
tion applies  when  the  normalizing  dynamic  head  ,s  specified  to  occur  at  the 
location  where  the  total  pressure  is  sought  (input  variable  NL2=0),  but  it 
is  generally  then  more  convenient  to  write  the  equation  in  the  form 

Pr  = T+  l(i’  - i^>r ! 


5.  ENTHALPY  AND  ROTHALPY,  AND  RELATED  DEFINITIONS 

The  Euler  turbomachine  equation  relates  changes  in  angular  n. omen- 
tum to  change  in  total  enthalpy  thus 

= Li^e-r^o,)  (46) 

Absolute  and  relative  whirl  velocities  are  related  by 

Vg  = W©  4 W 

Rothalpy  is  related  to  total  enthalpy  by 

I = H - ^rVe/^j 


(47) 

(48) 
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Relative  total  enthalpy  is  given  by 


Hy=H-(v;-(Ve-^)l)/^T 

also 

Hr,-  Hr,  = *>'  I 

Static  enthalpy  is  given  by 

li  ~ H - V 

or 

Ui+(M-vw 


6.  ANNULUS  WALL  BOUNDARY  LAYERS 

The  program  inchides  provision  to  compute  blockage  due  to  annulus 
wall  boundary  layers  by  a simple  attached  turbulent  boundary  layer  method. 
The  resulting  blockage  is  applied  as  a uniformly  distributed  factor  in  the 
continuity  equation.  (It  is  a component  of  the  term  B in  Equation  25;  the 
other  component  is  the  "wake  blockage",  specified  in  the  .nout  data  by 
B BLOCK  and  BDIST. ) 

The  boundary  layer  momentum  thickness  (on  each  wall)  is  calculated 

from 

^ r . 6.8 

Q = Vw34  \)°,fc(C,  + °.0(^V^cimJ  (53) 


(49) 


(50) 


(51) 


(52) 


Schlichting  gives  this  equation  (Reference  9,  Equation  22.8)  for  a 
turbulent  boundary  layer  on  a flat  plate  with  pressure  gradient,  except  that 
the  total  velocity  and  tlow  length  have  been  replaced  by  the  meridional 
velocity  component  and  meridional  distance  after  the  n.anner  of  Jansen, 
Reference  10.  The  equation  is  strictly  valid  for  a constant  shape  factor 
(ratio  of  displacement  to  momentum  thickness)  of  l.  4.  In  the  program, 
any  constant  value  may  be  used,  and,  for  multistage  machines,  a value 
around  0.7  seems  appropriate. 
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The  constant  Ct  is  determined  from  the  specified  blockage  at 
the  inlet  station  which  is  assumed  to  occur  as  boundary  layers  of  equal 
thickness  on  each  wall.  Then  the  integral  is  made  from  the  inlet  station 
to  each  station  in  the  flow.  If  a hub  does  not  exist  in  the  machine  at  the 
inlet,  all  the  inlet  blockage  is  assigned  to  the  casing.  Integration  on  the 
hub  then  commences  at  the  "bullet  nose"  where  the  boundary  layer  is 
assumed  to  have  zero  thickness. 

7.  LOSS  COEFFICIENT  RE-ESTIMATION 

The  loss  coefficient  re-estimation  procedure,  if  enacted,  consists 
of  determining  a diffusion  loss  coefficient  and  a shock  loss  coefficient,  if 
applicable,  and  summing  these  to  give  the  total  loss  coefficient.  This  is 
done  for  each  streamline. 

The  diffusion  loss  coefficient  is  found  by  first  calculating  the  diffusion 
factor  from  Equation  2.  The  corresponding  loss  parameter  is  then  inter- 
polated from  each  of  the  three  loss  parameter/ diffusion  factor  curves 
specified  to  be  used  for  the  blade  by  |NEVAL  |.  These  three  curves  give 
loss  parameters  for  10,  50,  a:.d  90  per  cent  of  the  radial  blade  height 
(measured  at  the  blade  trailing  edge),  and  the  value  for  the  particular 
streamline  radius  is  interpolated  from  the  three  values.  Then,  the  diffusion 
loss  coefficient  is  given  by 

/o  p ss  id  j>«ro  7.  (T  S?C  |?  ( 54) 


The  shock  loss  coefficient,  if  calculated  (NDEL  £ 0),  is  found  from 
an  average  Mach  number  at  the  presumed  strong  shock  location.  A suction- 
surface  Mach  number  is  found  by  expanding  (or  con. pressing)  the  relative 
inlet  Mach  number  through  an  expansion  angle  using  the  Prandtl-Meyer 
function. 


(55) 


The  expansion  angle  is  the  sum  of  the  value  specified  in  the  input 
(DELTA)  and  the  incidence  angle.  If  the  blade  geometry  is  not  given  at  the 
blade  leading  edge  (NDATA  =0),  the  incidence  is  assumed  to  be  zero.  If  the 
relative  inlet  Mach  is  less  than  unity,  the  expansion  is  calculated  for  an 
inlet  value  of  unity.  The  shock  Mach  number  is  then  taken  as  the  mean 
of  the  relative  inlet  Mach  number  (or  unity)  and  the  suction  surface  Mach 
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number.  If  the  relative  inlet  Mach  number  is  less  than  unity,  the  shock 
Mach  number  is  then  multiplied  by  the  inlet  value.  Should  this  final 
value  be  less  than  unity,  the  shock  loss  coefficient  is  set  to  zero.  Other- 
wise, it  is  calculated  from  the  pressure  loss  across  a normal  shock,  using 


y+i  y _ 


(■  + 


>-i  _ i 


As  mentioned  above,  the  final  loss  coefficient  is 


JO  sr 


8.  OFF-DESIGN  ANALYSIS  SCHEME  (NWORK  =7) 

The  loss  coefficient  is  defined  in  terms  of  a minimum  value, 

(input  variable  DATA2),  a reference  (or  minimum  loss)  relative  inlet  angle, 
, (DATA6),  and  "stalling"  and  "choking"  inlet  angles,  pic  and 
(DATA8  and  DATA 9).  The  loss  coefficient  at  any  inlet  angle  ^ is  then  given 
by 


id  = id, 


Mg 


ft  ~ ft  1 l£F 


PI 


where  , °r  Fu  is  used  depending  upon  whether  ft  is  greater  or  less 
than  fti t£f.  A maximum  value  of  0.  5 is  imposed  on  the  resulting  loss 
coefficient.  The  relative  outlet  flow  angle  varies  linearly  from  the  reference 
value  fttfcF(DATAl)  according  to 

Plt£F  + (fi-  1^5^  (59) 

The  gradient  is  given  by  the  input  variable  DATA7. 
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9. 


BLADE  GEOMETRY 


For  the  purpose  of  determining  relative  flow  angles,  deviation  angles, 
and  incidence  angles,  blade  angles  are  required  to  be  known  on  the  stream  - 
surface  plane. 

When  the  input  variable  NWORK  equals  6,  or  blade  leading  edge 
geometry  is  given,  NWORK  = 0 and  NDATA  >0,  the  angle  given  in  the  input 
(as  DATA1)  is  measured  on  the  cylindrical  plane.  The  corresponding  value 
on  the  stream  surface  is  found  from 


-tap-  ,60, 


This  equation  is  given  in  Reference  11.  The  angle  is  related  to 
the  input  blade  lean  angle  by 


10.  INTEGRATED  PERFORMANCE 

Machine  performance  is  shown  in  the  output  in  terms  of  pressure 
ratio  and  isentropic  efficiency  on  each  streamline,  and  also  as  "integrated" 
figures  for  the  whole  flow. 

In  order  to  determine  the  integrated  values,  the  following  definitions 
are  made.  At  any  station,  the  mean  total  enthalpy  is  given  by 

P = ■€  H.£w/W  <«) 


The  corresponding  mean  total  pressure  is  given  by 

p = xP.Sw/w 


(63) 
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From  such  values  for  the  inlet  station,  a mean  inlet  entropy  ( S\  ) 
is  found.  An  isentropic  outlet  total  enthalpy  ( is  then  found  as  a function 

of  the  outlet  mean  total  pressure  ( ) and  . The  integrated  isentropic 

efficiency  is  then  given  by 


The  inlet  and  outlet  mean  total  pressures  may  also  be  used  to 
determine  an  integrated  total  pressure  ratio. 

11.  RADIAL  TRANSFER  OF  BLADE  WAKES 

A phenomenon  of  the  flow  through  typical  compressor  rotor  blade 
rows  is  the  high  losses  recorded  for  blade  sections  near  the  tip.  This  is 
illustrated  in  Figure  2.3  of  Reference  6 where,  for  a given  diffusion  factor, 
rotor  tip  sections  are  shown  to  have  a loss  parameter  from  two  to  four 
times  larger  than  other  blade  sections.  In  compressor  design  and  analysis 
procedures,  the  traditional  method  has  been  to  use  different  lo s =>  parameter 
/diffusion  factor  curves  for  different  blade  section  locations.  The  analysis 
in  this  subsection  is  an  attempt  to  account  analytically  for  the  phenomenon. 

It  is  postulated  that  the  result  seen  comes  about  not  from  differences  in  the 
boundary  layer  development  mechanism , but  from  radial  transfers  of  blade 
boundary  layer  and  wake  flows,  that  is,  secondary  flows.  These  are  caused 
by  an  imbalance  of  the  radial  forces  on  the  slow -moving  fluid  particles,  and 
the  in. balance  is  the  primary  target  of  the  following  analysis.  Because  the 
resultant  radial  transfers  are  generally  much  greater  for  rotors  than  stators, 
the  effect  is  usually  only  easily  discernab.le  there,  but  the  analysis  may  be 
applied  uniformly  to  rotors  and  stators.  If  the  radial  transfers  are  properly 
accounted  for,  the  same  wake  generation  model  (loss  coefficient  correlation) 
may  be  used  for  all  blade  sections,  provided  it  can  be  assumed  that  boundary 
layer  development  on  any  blade  section  is  affected  at  most  to  a second  order 
only  by  the  radial  transfers. 

The  approach  followed  is  to  determine  the  ratio  of  radial  shear  stress 
to  axial  shear  stress  on  an  element  of  the  blade  surface  boundary  layers. 

It  is  then  postulated  that  the  extent  of  the  radial  movement  of  the  boundary 
layer/wake  material  is  proportional  to  the  force  ratio.  The  radial  force  is 
found  by  considering  the  static  pressure  and  centrifugal  forces  acting  upon 
a small  volume  of  blade  surface  boundary  layer.  This  is  then  related  to  the 
shear  stress  at  the  blade  surface  by  the  boundary  layer  thickness.  When 
the  stress  ratio  is  formed,  the  boundary  layer  thickness  disappears  from 
the  final  result. 
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The  radial  stress  is  found  from  consideration  of  the  forces  acting 
on  a small  volume  of  the  boundary  layer,  which  show  the  net  radial  force 
to  be 


r>  1 I 

<Y  = «-  r - f ft 


(65) 


The  small  volume  is  in  contact  with  the  blade  surface  and,  therefore, 
acts  at  the  surface  to  produce  a shear  stress.  If  the  boundary  layer  thick- 
ness is  taken  to  be  the  conventional  displacement  thickness,  then  the  shear 
stress  is  given  by 


% = H) 


(66) 


The  axial  stress  is  found  as  follows.  For  a turbulent  boundary  layer 
without  pressure  gradient,  Schlichting  (Reference  9,  Equation  21.9)  gives 
the  momentum  thickness  as 


6 = -ottxOV) 


-Vs 


(67) 


and  the  shear  stress  is 


/V 


/)!  ^ 
y)  I/® 


(68) 


At  the  blade  trailing  edge  then,  the  axial  stress  is 


T = , LL  4 - 


(69) 


The  ratio  of  the  stresses  is 


r " % 

' % 


- It) 

ffii  1 1 


/ 


(70) 
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K. 


If  the  shape  factor  is  constant,  we  can  write 


(7’) 


It  is  now  assumed  that  the  tangent  of  the  angle  of  the  boundary  layer/wake 
movement  is  proportional  to  F#  . Thus  a suitable  axial  dimension  will 
yield  the  dimensional  radial  movement.  For  this  dimension,  the  blade 
chord  is  taken.  Hej.  e: 


F 


(72) 


For  our  purposes,  the  free  stream  velocity  is  the  relative  flow  velocity. 
The  constant  K remains  to  be  determined  from  experimental  data.  In 
conjunction  with  the  method  of  utilizing  F that  is  described  below,  a value 
of  2.  0 appears  to  be  approximately  correct. 


Equation  72  gives  a radial  movement  distance  at  any  radius,  but 
some  further  assumptions  are  required  before  this  gives  a loss  profile 
modified  for  the  effects  of  radial  transfers.  Relative  total  pressure  loss 
coefficient  is  approximately  directly  proportional  to  the  thickness  of  the 
wake  leaving  the  blade,  and  therefore  the  radial  movement  F is  directly 
applied  to  a two-dimensional  loss  coefficient  distribution  to  produce  the 
corresponding  three-dimensional  distribution.  Inputs  to  the  numerical 
procedure  are  loss  coefficients  presumed  to  be  based  upon  blade  element 
characteristics  (and  therefore  designated  two-dimensional)  for  each  stream 
line,  and  F-values  computed  for  each  streamline.  These  two  quantities 
are  spline-fit  interpolated  at  150  points  uniformly  distributed  between  hub 
and  casing.  Each  pair  of  values  is  taken  to  apply  to  an  interval  that  lies 
about  the  corresponding  radial  location.  Thus  at  radius  Fj  , there  is  a loss 
coefficient  and  a shift  Fj  that  apply  to  interval  of  width  bj  , where 

j = (Tj+i  “ V'  (73) 


In  the  case  where  _j  - 1 or  150,  bj  will  be  half  the  size  of  the  other  interval 
The  loss  coefficient  from  interval  bj  is  redistributed  according  to  Fj  . It 
is  assumed  that  the  new  distribution  is  a uniform  division  of  th<-  loss  over 
the  distance  ( |Fj|  ),  so  that  in  that  region  the  new  loss  coefficient 
attributable  to  till,  is 


(74) 


The  distributed  loss  coefficient  is  seen  at  all  of  the  150  radii  that  the 
movement  Fj  crosses,  including  the  starting  point  ITj  . In  the  event  that 
the  movement  Fj  prescribes  a shift  outside  of  the  compressor  annulus,  Fj 
is  reduced  so  that  the  movement  just  reaches  the  annulus  wall.  Then 
is  proportionally  larger  than  it  would  otherwise  have  been.  This  gives  a 
"piling  up"  effect  near  a wall  towards  which  there  is  boundary  layer  move- 
ment. The  total  loss  coefficient  at  any  radius  is  given  by  summing  all  of  the 
10 3>  components  that  apply  to  that  radius.  The  last  step  is  to  determine  loss 
coefficients  for  each  of  the  streamlines  that  are  being  used  in  the  basic 
flowfield  calculation,  and  for  which  the  original  two-dimensional  loss 
coefficients  were  given.  To  accomplish  this,  an  eighth-order  polynomial 
least-squares  fit  of  the  150(c^values  is  made,  and  then  loss  coefficients 
at  the  desired  radii  are  taken  from  the  polynomial  representation.  The 
least-squares  fit  was  introduced  into  the  procedure  because  the  150-point 
curve  exhibits  some  scatter  and  so  cannot  be  used  directly  as  input  to  a 
collating-curve  interpolation  scheme. 

The  procedure  described  above  gives  results  that  are,  for  engineering 
purposes,  independent  of  the  number  of  streamlines  used  to  describe  the 
flowfield.  The  selection  of  the  number  of  points  (150)  at  which  the  three- 
dimensional  loss  coefficients  are  derived,  and  the  use  of  an  eighth -order 
polynomial  fit  for  the  three  dimensional  loss  coefficient  distribution  do  have 
son.e  effect  on  the  results.  As  the  number  of  points  used  to  derive  the  loss 
coefficient  distribution  is  reduced,  the  scatter  of  the  resulting  curve 
increases,  and  effectively  less  radial  transfer  occurs.  This  is  because 
fewer  of  the  shifts  will  be  large  enough  to  enter  an  adjacent  accounting  area. 
For  larger  numbers  of  points  this  effect  diminishes,  and  there  is  generally 
little  difference  between  the  results  for  100,  150,  and  200  points.  An 
eighth-order  polynomial  fit  normally  retains  all  the  essential  features  of 
the  150-point  distribution,  without  being  unduly  influenced  by  the  scatter. 

12.  TURBULENT  MIXING  OF  THE  AXISYMMETRIC  FLOW 

Assumptions  that  are  usually  made  in  turbomachine  flowfield  analys:s 
are  that  the  flow  is  inviscid  and  nonconducting,  so  that  there  is  no  transfer 
of  mass  or  energy  across  streamlines  of  the  flow.  The  momentum  equations 
derived  earlier  (in  Section  V.  1)  incorporate  these  assumptions  as  well  as 
the  assumption  that  the  flow  ic  axisymmetric.  In  reality,  the  flow  is 
viscous  and  turbulent  in  nature, and  the  effects  of  this  are  incorporated 
into  the  computation  scheme. 
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The  momentum  equations  are  obtained  by  making  the  well-known 
assumption  that  the  flow  may  be  described  by  the  Navier-Stokes  equations 
for  laminar  flow,  with  the  molecular  kinematic  viscosity  replaced  by  an 
appropriate  eddy  viscosity.  (See,  for  example,  Reference  9.)  For  the 
purpose  of  analyzing  the  turbulence,  the  flow  is  assumed  to  be  contained 
by  frictionless  walls.  Thus  the  result  will  be  the  effects  of  mixing  of  the 
flow  without  regard  for  annulus  wall  boundary  layers.  In  fact,  the  entire 
flow  is  assumed  to  be  of  a boundary  layer  nature.  Thus  the  momentum 
equations,  for  axisymmetr ic,  steady  flow,  and  incorporating  the  boundary 
layer  type  flow  assumption  .which  eliminates  second  derivatives  of  velocity 
in  the  streamwise  direction,  are  as  follows: 


It  is  now  found  convenient  to  consider  the  flow  to  be  contained  by  a cylin- 
drical annulus.  By  making  this  assumption,  it  follows  that  the  radial 
velocity  component  must  everywhere  be  very  small  compared  to  the  other 
two  velocity  components,  so  that  Equations  75,  76,  and  77  may  be  written: 


(78) 


(79) 
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(80) 


It  is  interesting  to  note  that  Equation  78  is  the  Simple  Radial  Equilibrium 
Equation  of  turbomachinery,  whereas  Equations  79  and  80  retain  viscous 
terms.  TLiS  justifies  a feature  of  the  overall  method  that  is  introduced 
later:  the  assumption  that  even  when  turbulent  mixing  effects  are  incor- 
porated into  computation,  the  ir.viscid  momentum  equations  derived  in 
Section  V.  1 may  be  used  to  compute  the  velocity  profiles  at  each  computing 
station. 


The  boundary  condition  of  frictionless  walls  may  be  expressed  as: 


(81) 


In  order  to  determine  the  effect  of  turoulence  upon  the  total  enthalpy, 
an  energy  equation  is  required.  This  is  obtained  from  the  energy  equation 
for  laminar  flow,  neglecting  conduction  in  the  streamwise  direction.  Eddy 
viscosity  is  again  substituted  for  molecular  viscosity,  and  the  turbulent 
Prandtl  number  is  taken  to  be  unity.  Then  the  energy  equation  (for  a 
cylindrical  annulus)  is: 


£ 
v 5V 


j I- 1) 


(82) 


The  adiabatic  wall  boundary  condition  may  be  expressed  as 


V 


= o 


at  the  walls 


(This  assumes  that 
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(83) 


Continuity  is  satisfied  at  each  axial  location  by  use  of  the  relation: 
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(84) 


Equations  78,  79,  80,  82,  and  84  with  the  boundary  conditions  given 
by  Equations  81  and  83  describe  the  selected  model  of  turbulent  flow,  and 
are  used  to  determine  changes  in  the  distributions  of  entropy,  total  enthalpy, 
and  angular  momentum.  The  numerical  procedure  that  is  used  is  described 
in  Section  VI.  10.  The  ecdy  viscosity  is  presumed  to  be  constant  throughout 
the  flowfield.  This  may  be  partially  justified  by  considering  that  in  the 
typical  situation  where  an  eddy  viscosity  model  is  used,  (such  as  the  mixing 
of  a free  jet,  or  the  shear  layer  between  two  parallel  flows),  there  is  often 
a step  change  in  vel<  city  at  the  upstream  boundary,  whereas  in  the  case 
of  interest  here  the  gradients  in  velocity  are  normally  similar  in  magnitude 
throughout  the  flow.  Application  of  the  constant  eddy  viscosity  model  to 
cases  with  step  changes  in  conditions  at  the  upstream  boundary  could  be 
expected  to  yield  somewhat  inaccurate  profiles  of  conditions  downstream, 
although  the  basic  trends  would  be  reproduced. 
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SECTION  VI 


NUMERICAL  PROCEDURES  FOR  AERODYNAMIC  ANALYSIS 


1.  INTEGRATION  OF  DESIGN  MOMENTUM  EQUATION 


The  design-case  momentum  equation  (Equation  19,  20,  23,  or  24, 
as  appropriate)  is  integrated  as  a lineai,  first  order  equation  in  Vtn* 
by  assuming  that  only  Vm  is  unknown.  Then  the  integrating  factor 
technique  may  be  applied  to  give: 


v»; = v»;  (-A  f 


where,  in  the  case  of  Equation  19, 
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For  Equation  20,  A remains  unchanged,  and  the  last  two  terms  in  S 
are  replaced  by: 

Su,(^4iS) 


3 d F* 


For  Equation  23,  the  following  definitions  apply: 
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For  Equation  24,  fa  remains  unchanged,  and  the  last  two  terms  in  B are 
replaced  by: 
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Equation  85  fails  if  fa  - O,  but  then  the  following  relation  applies 


V.,'-  ».,■* 


In  the  computer  program,  the  switch  from  Equation  85  to  Equation  90  is 
made  when  I A I < io 

Equation  85  (or  Equation  90)  is  applied  to  the  interval  between  each 
adjacent  pair  of  streamlines  in  turn.  A single  value  of  the  parameters  fl 
and  6 is  determined  for  each  interval.  Gradients  in  the  direction  of  the 
computing  station  are  found  using: 


dF  _ M, 
dl  V 1, 


Gradients  in  the  streamline  direction  are  found  from,  (for  Station  I), 


4 = f j ll 

dm  \ Wn,  - mr  - w;-,  1 ' 
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At  the  first  and  last  computing  stations,  and  during  the  lirst  pass  of 
calculation,  the  streamwise  gradients  are  assumed  to  be  zero.  A value  is 
determined  from  Equation  92  for  each  streamline,  and  then  an  average  value 
is  used  in  each  streomtube.  This  procedure  is  followed  for  each  quantity 
that  is  given  on  a streamline  rather  than  for  a streamtube.  The  streamline 
slope  angle  ^ is  given,  at  Station  I,  by 


Xi*,  — Xi 
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At  the  first  (last)  computing  station, 


* 


is  given  by 


+ I y. 


(94) 


The  radius  of  curvature  of  the  streamlines  is  assumed  to  be  zero  at 
the  first  and  last  stations  and  elsewhere  is  determined  from: 

inVj-l ii \k 
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The  n.ethod  of  solution  described  above  is  based  upon  the  assumption 
that  the  parameters  A and  B are  fixed.  In  fact,  they  are  dependent  to  a 
varying  degree  upon  the  meridional  velocity,  and  so  the  determination  of 
the  meridional  velocity  is  iterative.  meridional  velocity  profile  is 
estimated,  the  parameters  A and  B are  determined,  and  a new  velocity 
profile  is  determined  from  the  mid-streamline  velocity  and  the  A and  B 
parameters.  Some  instability  in  this  procedure  was  observed  on  occasion, 
and  therefore  a relaxation  factor  R has  been  introduced  into  the  procedure 
for  the  new  estimate  of  the  meridional  velocities,  thus; 
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where 
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The  nondimensional  change  in  V*  in  Equation  97  is  computed  as  an  - 
coordinate  weighted  average,  and  R is  not  permitted  to  be  less  than  0.  1. 

Some  limits  are  set  on  the  values  of  Vrtl  that  may  be  recorded  for 
each  streamline.  A lower  limit  of  1.0,  and  an  upper  limit  of  twice  the 
value  on  the  n.id-streamline  are  imposed,  if  necessary.  For  the  case  of 
Subroutine  UD0308  only,  there  is  also  an  upper  limit  equal  to  the  local 
speed  of  sound. 

When  a new  meridional  velocity  profile  has  been  estimated  by  the 
above  procedure,  the  entropy  distribution  is  recalculated  if  losses  are 
specified  in  such  a way  that  the  entropy  is  a function  of  meridional  velocity. 
Then  the  continuity  summations  described  separately  below  are  made  and, 
if  continuity  is  not  satisfied,  or  the  resulting  meridional  velocity  profile 
differs  from  that  initially  assumed,  a new  iteration  is  started.  The 
maximum  number  of  permitted  iterations  is  controllea  by  the  variable 
ITMAX,  which  is  set  (separately)  to  20  in  both  Subroutine  UD0308  and 
Subroutine  UD0326. 

2.  INTEGRATION  OF  ANALYSIS  MOMENTUM  EQUATION 

The  analysis  momentum  equation  (Equation  17,  18,  21,  or  22,  as 
appropriate)  is  nonlinear  in  V**i  , and  therefore  a different  scheme  to  that 
described  above  for  the  design  case  equation  is  used.  All  the  gradients  and 
streamline  characteristics  are  computed  as  for  the  design  case.  Using  an 
estimate  of  the  average  meridional  velocity  in  a given  streamtube,  the 
gradient  dV-'/j-t  is  computed  from  either  Equation  17,  18,  21,  or  22. 
This  gradient  is  then  used  to  step  to  the  next  streamline  using: 


V-  V.,‘ . 


(98) 


The  velocity  determined  by  Equation  98  is  not  permicted  to  be  less  ihan  1.0 
nor  greater  than  thtee  tin.es  the  mid-streamline  value.  Because  the  gradient 
depends  upon  the  meridional  velocity,  tne  velocity  oetermination  by  chis 
means  is  iterative,  and  up  to  10  loops  are  allowed. 

When  a new  meridional  velocity  profile  has  oeen  established  by  the 
above  procedure,  the  entropy  distribution  is  recalculated  if  losses  are 
specified  in  such  a way  that  the  entropy  is  a function  of  meridional  velocity. 
Then  the  continuity  summations  described  sepavately  below  are  made  and, 
if  continuity  is  not  satisfied,  or  the  resulting  meridional  velocity  profile 
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differs  from  that  initially  assumed,  a new  iteration  is  started.  The  maxi- 
mum number  of  permitted  iterations  is  controlled  by  the  variable  ITMAX, 
which  is  set  (separately)  to  20  in  both  Subroutine  UD0308  and  Subroutine 
UD0326. 

3.  SATISFACTION  OF  CONTINUITY  EQUATION 

The  continuity  equation  (Equation  25)  is  required  to  be  satisfied  at 
each  computing  station  simultaneously  with  the  appropriate  momentum 
equation.  Accordingly,  the  following  procedure  is  incorporated  into  the  up 
to  20  iterations  that  may  be  performed  to  satisfy  the  momentum  equation. 


The  flow  rate  is  determined  from  Equation  25  using  trapezoidal 
integiation,  and  the  rate  of  change  of  flow  with  mid-streamline  meridional 
velocity  is  similarly  calculated  from  Equation  29  or  33,  as  appropriate. 
Several  courses  of  action  may  follow. 


dw/JV 


If  the  solution  type  is  that  desired  (that  is  " /dVt»  is  positive  when 
a subsonic  solution  is  specified  or  negative  when  a supersonic  solution  is 
specified),  then  the  meridional  velocity  increment  to  satisfy  continuity 
is  estimated  from 


Wsfreofitd  ~ Wc*ici«Ut«j 
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(99) 


The  velocity  increment  is  limited  to  ±10%  of  the  previous  m id- 
streamline  value.  Also,  if  the  flow  rate  is  less  than  the  desired  value, 
and  a valid  solution  to  the  momentum  equation  was  previously  computed, 
the  mid-streamline  meridional  velocity  is  noted,  as  it  is  a value  below 
which  the  correct  value  cannot  ' e for  the  subsonic  case,  or  above 
which  the  correct  value  cannot  lie  for  the  supersonic  case. 


If  the  solution  type  is  not  that  desired,  the  velocity  increment  is 
arbitrarily  set  to  plus  or  minus  10%  of  the  mid-streamline  value  in  an 
attempt  to  switch  to  the  supersonic  or  subsonic  case,  respectively.  Also, 
if  a valid  solution  to  the  momentum  equation  was  previously  computed, 
the  mid-streamline  velocity  is  noted  because,  for  example,  when  a sub- 
sonic solution  is  desired,  no  meridional  velocity  should  be  evaluated  which 
is  known  to  produce  a supersonic  solution. 

When  an  analysis-type  momentum  equation  is  being  used  only,  a 
check  is  now  made  on  the  velocity  increment  that  has  been  determined 
by  either  of  the  methods  described  above.  Suppose  that  a subsonic  solution 
is  desired.  (The  logic  is  "syn.mctr ical"  for  the  case  of  a desired  super- 
sonic solution.)  The  mid-strean  line  meridional  velocity  that  would  be 
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given  by  modifying  the  previous  value  by  the  calculated  increment  is 
checked  to  see  if  it  would  be  higher  than  the  lowest  value  which  has  been 
shown  to  imply  a supersonic  solution,  V*h  . If  so,  a new  increment 

is  calculated.  The  new  increment  is  given  by 

AYh=  (Vus^,,  - Vnolj)/ 2 (100) 

Thus,  the  new  velocity  will  be  within  the  range  of  possible  values, 
and  will  generally  refine  the  previously  established  limits  for  the  range. 

A similar  check  is  made  with  respect  to  the  minimum  plausible 
velocity,  that  is  V*  , the  highest  velocity  shown  to  imply  a subsonic 

solution  at  a flow  rate  lower  than  desired. 

The  meridional  velocities  on  each  streamline  are  incremented  by 
the  amount  derived  as  described  above,  and  a lower  limit  of  1.0  is  imposed 
on  the  value  produced. 

If  Subroutine  UD0308  is  being  used  (input  variable  NEQN  = 0 or  1), 
an  upper  limit  equal  to  the  local  sound  speed  is  also  imposed  upon  the 
meridional  velocity  at  each  streamline.  In  this  case,  the  maximum  mid- 
sireamline  meridional  velocity  allowable  for  a supersonic  solution 
( Vm  ) 's  set  equal  to  (the  sound  speed  x 1.05)  initially. 

When  Subroutine  UD0326  is  being  used  (input  variable  NEQN  - 2 or  3) 
V is  initialized  to  Z500.  0.  In  both  cases,  the  minimum 
permitted  velocity  for  a subsonic  solution  is  initialized  ta  zero. 

A check  is  made  on  the  convergence  of  the  momentum  and  continuity 
equat.ons,  and,il  these  do  not  both  converge.  a new  iteration  is  started, 
•uibject  to  the  maximum  number  of  iterations  permitted  by  the  variable 
ITMAX.  If  a new  iteration  is  started,  tire  entropy  distribution  is  recalcu- 
lated if  it  is  a function  of  meridional  velocity,  and  then  the  momentum 
equation  integration  p'acedure  is  entered  as  described  above. 

■i.  STREAMLINE  RE LO CA  TION  1TE  RATION 

The  positions  of  the  streamline  are  determined  iteratively  using  the 
condition  that  a constant  fraction  of  the  total  flow  should  lie  in  each  stream- 
tube.  The  flow  distribution  used  to  revise  the  streamline  coordinate  is  a 
function  of  the  streamline  characteristics,  and,  therefore,  the  calculation 
is  potentially  unstable.  -A  relaxation  factor  is  derived  in  the  following 
analysis,  which  is  based  upon  a similar  analysis  performed  by  Wilkinson 
for  two-dimensional  flow,  Refeience  12. 


Consider  a swirling  flow  in  a cylindrical  annulus,  with  no  enthalpy 
or  entropy  gradients.  Computing  stations  are  uniformly  spaced  Ax  apart, 
and  are  radial.  Equation  20  may  be  applied,  and  simplified  to  give 
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By  assuming  orly  small  deviations  from  midradius  conditions,  we 
can  approximate  Equation  101  by 


dV* 

Tr  = 
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If  we  introduce  an  error  in  the  streamline  radius  a*  one  station  only, 
then  the  corresponding  radius  of  curvature  is 
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Combining  Equations  102  and  103  gives 
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or 
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Mass  flux  distribution  is  proportional  to 
iniate  the  distribution,  using 


and  again,  we  appiox- 
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and 


Q = (/ 


(107) 


For  a high  radius  ratio,  the  flow  between  radii  f,  and  T*  is  propor- 

rfl 


tional  to 


and  from  the  previous  equations 


^ l>\ M h 

(m  = 0 + 

( (*l  — - ( f,  — Y*teO  ) j 


(108) 


To  determine  the  effect  of  the  error,  we  compare  the  predicted 
location  of  the  screamline  that  divides  the  flow  into  two  equal  parts  with  its 
correction  location,  that  is  when  E=0  The  predicted  location,  , may 
be  found  from  the  following  equation 
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Using  Equation  108,  and  assuming 
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this  yields 
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where 
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The  correction  location,  Y9  , is  similarly  found  to  be 
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Thus,  the  shift  in  streamline  position  is 


t(l  + Vr*) 


(114) 


The  shift  to  correct  the  error  E should  be  - E , so  that  the  correspon- 
ding relaxation  factor  R , is 
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The  case  considered,  that  is  where  the  error  E is  in. posed  at  one 
station  only,  is  one  extreme.  It  produces  the  largest  streamline  curvatures, 
and,  without  a relaxation  factor,  the  largest  shifts  m streamline  radius  from 
one  pass  to  the  next.  The  other  extreme  is  when  the  error  is  imposed  at  all 
stations,  producing  no  curvature.  In  this  case,  any  damping  (that  is,  a 
relaxation  factor  of  less  than  unity),  will  cause  the  streamline  shift  to  be 
underestimated.  An  optimum  relaxation  factor  F is  defined  as  one  which 
reduces  the  error  to  the  same  fraction  for  both  cases.  Noting  the  sign 
change  occurring  between  the  two  cases,  we  obtain 
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so  that 


h = 
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The  remaining  tern.  Q is  determined  by  assuming  the  fluid  to  be  a perfect 
gas.  Then 
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For  the  case  when  the  relative  flow  angle  is  specified,  the  analysis  starts 
from  Equation  18,  and,  following  the  same  steps  as  for  the  case  above, 
the  final  result  is 
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In  the  program,  the  length  Aft  is  taken  as  the  computing  station 
length  and  the  length  Ax  is  taken  as  the  meridional  distance  at  mid-annulus 
to  the  nearer  of  the  two  adjacent  computing  stations.  The  constant  8 in 
Equations  117  and  120  has  been  found  fron  experience  to  often  need  reducing 
to  6 to  ensure  stability.  It  is  input  as  the  variable  RCONST.  The  Mach 
number  squared  that  appears  in  Equations  119  and  121  is  limited  in  the 
program  to  the  value  input  as  XMMAX,  and  0.  6 has  been  found  to  be 
generally  satisfactory.  The  relaxation  factor  finally  determined  is  applied 
to  the  streamline  coordinate  -£L  , thus 

Vuv,  = "£o(ol  + F (^C  JcJJJ  - e(j ) (122) 


5.  CONVERGENCE  CRITERIA 

The  preceding  four  subsections  have  all  described  iterative  procedures 
involved  in  the  overall  determination  of  a solution.  The  convergence  of  each 
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of  these  procedures  is  measured  by  n.eans  of  a tolerance,  applied  as 
described  below,  and  related  to  a single  input  variable,  TOLNCE.  The 
normal  value  for  this  is  0.001. 

When  the  meridional  velocity  at  a station  is  determined  by  inte- 
gration of  the  analysis  form  of  the  momentum  equation,  the  velocity  at  any 
streamline  is  iterati  ely  determined  in  a stepping  procedure  from  the 
adjacent  streamline  that  is  nearer  to  the  mid-streamline.  This  iteration  is 
considered  converged  if  the  calculated  velocity  repeats  to  within  a tolerance 
of  TOLNCE/ 5.0,  applied  nondin  ensionally.  Thus,  :f  TOLNCE  = 0.001.  two 
successive  meridional  velocities  must  be  within  one  part  in  5000. 

Following  adjustment  of  the  velocity  level  on  all  strean  lines,  a check 
is  made  to  see  if  both  continuity  and  momentum  are  satisfied.  The  toler- 
ance used  is  TOLNCE/5.0.  unless  ? loss  coefficient  re-estimation  is 
specified  at  the  station  (NEVAL  > 0),  when  the  tolerance  is  TOLNCE/10.0. 
The  continuity  check  is  that  the  calculated  flow  must  equal  the  specified 
value  to  within  the  stated  tolerance,  applied  nondin, ensionally.  Then  the 
n omentum  equation  solution  is  checked  by  requiring  that  the  velocities 
equal  those  estimated  to  within  the  same  tolerance,  applied  nondimensionally. 

If  loss  coefficient  re -estimation  is  specified  at  a computing  station, 
an  extra  iteration  for  the  loss  coefficients  is  involved.  Convergence  of  this 
iteration  is  measured  by  comparing  the  estimated  meridional  velocities  with 
those  calculated,  and  the  tolerance  applied  is  TOLNCE/5.0.  With  LITER 
set  to  say  5,  this  iteration  will  probably  not  converge  to  within  the  above 
stated  tolerance  during  the  first  few  passes  through  the  program. 

Overall  convergence  of  the  solution  is  determined  by  two  checks  that 
are  made  at  the  end  of  each  pass.  The  meridional  velocities  at  each  mesh 
point  n.ust  repeat  to  within  a tolerance  of  TOLNCE.  applied  nondin. ensionally. 
Also,  the  proportion  of  the  total  flow  that  passes  between  the  hub  and  each 
streamline  unit  be  witi  .n  TOLNCE  of  the  desired  value  at  every  computing 
station.  For  example,  if  some  streamline  should  have  0.  5 of  the  flow 
between  it  and  the  hub,  and  TOLNCE  = 0.  001,  the  acceptable  range  is  from 
0.  499  to  0.  501. 

6.  PRANDTL-MEYER  FUNCTION 

The  Prandtl -Meyer  function  is  given  by  Equation  55,  and  in  the  course 
of  re-estimating  the  loss  coefficients,  it  this  is  specified  to  be  done.it  is 
desired  to  determine  the  Mach  number  corresponding  to  a given  angle,  & 

This  is  done  iteratively  using  Newton's  method,  with  the  quantity  J l"Nz 
being  sought  (rather  than  M ).  An  initial  estimate.  Fj  , of  J )—  NL  is  made 

equal  to>/|-H,  +0.1  where  is  the  Mach  number  before  the  expansion. 

Then  the  new  estimate,  » i gi''en  by 
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This  is  repeated  until  Fj  and  F^are  within  0.  000001  (subject  to  a maximum 
number  of  iterations  of  20). 

7.  STAGNATION  POINT  PROVISION 

If  the  "hub”  radius  at  the  inlet  station  is  zero,  a stagnation  point  is 
incorporated  into  the  solution  at  the  last  station  for  which  the  "hub”  radius 
is  zero.  This  must  occur  upstream  of  any  blade  rows,  and  the  whirl 
velocity  at  zero  radius  must  be  zero.  The  stagnation  point  is  handled  by 
imposing  the  zero  meridional  velocity  at  that  point,  and  only  integrating 
the  momentum  equation  inwards  as  far  as  the  second  streamline.  The  con- 
tinuity summation  includes  the  innermost  streamline  in  the  normal  n.anner. 

A finite  curvature  will  be  shown  in  the  results  at  the  stagnation  point,  but 
this  is  of  no  consequence  because  the  momentum  equation  is  not  computed 
in  the  innermost  streamtube. 

8.  INTERPOLATION 

The  various  tabulated  data  input  to  the  program  is  interpolated  at  the 
streamline  by  the  method  specified  by  the  input  variable  NTERP.  Two 
alternatives  are  available.  If  NTERP  = 1,  linear  point-to-point  interpolation 
(or  extrapolation)  is  performed.  This  is  not  to  be  preferred  generally, 
assuming  the  data  points  define  a smooth  curve. 

If  NTERP  = 0,  spline-fit  interpolation  is  performed.  The  spline 
curve  consists  of  a series  of  algebraic  cubics,  one  equation  applying  between 
each  adjacent  pair  of  data-points.  The  coefficients  are  determined  such  that 
the  resultant  curve  passes  through  each  point,  and  has  continuous  first  and 
second  derivatives.  Various  boundary  conditions  may  be  imposed  at  each 
end  of  the  curve;  in  this  case  the  second  derivative  at  the  end  point  is  set 
to  zero. 

9.  DEVIATION  ANGLE  DETERMINATION  FOR 
ANALYTIC  MEANLINE  BLADE  SECTION 

The  determination  of  the  deviation  angles  for  analytic  meanline  blades 
is  performed  in  Subroutine  UD0319.  A description  of  the  procedure  used  is 
included  in  this  section.  The  subroutine  is  only  entered  if  the  input 
variable  NAERO  = 1. 
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The  trailing  edge  ieviation  is  calculated  from 


s $o(0  Ksi,  Kfcc  + — j^111  + 0'5[%-  0*5) j + 


(124) 


The  parameters  K.Sn  , ^/c.  . and  aAM  are  given  in  the  input  as  variables 

XKSHPE,  AC,  and  DEL  TAD,  respectively.  The  camber  ( © ) and  solidity 
( ) are  derived  from  the  basic  determination  of  the  blade  sections.  The 

remaining  parameters,  £olo  , Ktt  . b » and  nt  , are  found  by  inter- 
polation from  data  given  in  Figures  161,  172,  164,  and  195  of  Reference 
6,  respectively.  Equation  124  is  based  upon  Equations  269  and  271  of 
Reference  6,  except  that, when  a parameter  is  shown  in  Reference  6 to  be 
a function  of  relative  air  inlet  angle,  blade  inlet  angle  is  substituted  in 
this  program. 


For  locations,  that  is, computing  stations,  within  a blade  row,  the 
deviation  is  determined  from 


$ = S+e  ( ) 


(125) 


The  fraction  ° Is  interpolated  from  the  input  variables  DVFRAC  and 
DM  as  a function  of  the  location  on  the  meridional  chordline  and,  i* 

NRAD  > 1,  as  a function  of  streamsurface  radius  at  the  trailing  edge. 

The  "station  control  card",  blade  speed,  blade  geometry,  loss  values 
(interpolated  at  the  appropriate  radii),  and  deviation  angles  are  output  onto 
the  file  specified  by  unit  LCG5  for  each  computing  station  within  the  blade, 
or  at  the  blade  trailing  edge,  in  the  format  required  by  the  aerodynamic 
section.  The  station  control  card  has  NDATA  set  equal  to  the  number  of 
streamsurfaces  on  which  blade  sections  are  produced,  NWORK  = 6,  and 
NDIMEN  = 0.  The  remaining  variables  are  set  according  to  the  values  in 
the  input  data  for  that  station. 


At  the  blade  leading  edge,  no  deviation  calculation  is  required,  but  a 
"station  control  card",  a speed  of  zero,  and  blade  geometry  are  output  onto 
the  file  specified  by  unit  LOG5  if  the  input  variable  IFANGS  for  the  leading 
edge  station  is  1.  The  station  control  card  has  NDATA  set  equal  to  the 
number  of  streamsurfaces  on  which  blade  sections  are  produced,  and  all 
remaining  variables  on  the  control  card  are  set  to  zero,  except  NGUT1, 
NOUT2,  and  NOUT3, which  are  reproduced  as  read  in. 
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10. 


TURBULENT  MIXING  CALCULATION  PROCEDURE 


The  equations  derived  in  Section  V.  12  describe  the  changes  in  the 
distribution  of  conditions  along  "streamlines"  of  an  axi symmetric,  turbulent 
flow  in  a cylindrical  annulus.  In  Section  V.  1,  equations  we  re  derived  for 
the  stationwise  variation  in  meridional  velocity  of  an  axisymmetric,  inviscid 
flow  in  an  arbitrary  annulus.  These  show  the  meridional  velocity  gradient 
to  be  a function  of  the  distributions  of  entropy,  total  enthalpy,  and  angular 
momentum  (as  well  as  the  streamline  characteristics).  In  order  to  incor- 
porate the  mixing  calculations  into  the  overall  computations,  first  the 
turbulent  flow  equations  are  used  to  determine  modifications  to  the  station- 
wise  distributions  of  entropy,  total  enthalpy,  and  angular  momentum. 

Then,  the  modified  distributions  are  used  in  an  inviscid  momentum  equa- 
tion to  compute  the  meridional  velocity  profile.  This  method  is  used  for  the 
interval  between  each  successive  pair  of  computing  stations  from  the  first 
to  the  last,  for  each  pass  through  the  iterative  procedure  to  determine 
the  streamline  locations,  and  hence  the  final  solution.  On  a step-by-step 
basis,  the  computation  proceeds  thus: 

1)  At  Station  1,  an  inviscid  momentum  equation  gives  the  meridional 
velocity  distribution. 

2)  The  turbulent  flow  equations  are  applied  for  the  interval  Station 
1 to  Station  2 to  give  the  distributions  of  entropy,  total  enthalpy, 
and  angular  momentum  at  Station  2. 

3)  If  a blade-row  lies  between  Stations  1 and  2,  the  distributions 
of  entropy,  total  enthalpy,  and  angular  momentum  at  Station  2 
are  redetermined,  using  the  modified  distributions  determined 
in  step  2 as  the  inlet  conditions.  Thus  the  effects  of  the  blade - 
row  (if  one  exists)  are  superimposed  upon  the  mixing  calculation. 

4)  An  inviscid  momentum  equation  is  used  to  determine  the 
meridional  velocity  distribution  at  Station  2. 

5)  Steps  2,  3,  and  4 are  repeated  for  the  intervals  Station  3 to  4, 

4 to  5,  and  so  on,  to  the  last  computing  station. 

6)  A new  streamline  pattern  estimate  now  exists,  based  upon  the 
meridional  velocity  profiles  computed  for  each  station.  Unless 
the  overall  solution  is  converged,  control  returns  to  Step  1 for 
another  pass. 

The  application  of  the  turbulent  flow  equations,  Step  2 above,  is  made 
as  follows.  The  interval  cf  con. pressor  annulus  is  approximated  by  a 
cylindrical  annulus,  for  which  the  hub  and  tip  radii  are  the  mean  hub  and 
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tip  radii  of  the  segment  of  the  compressor  ai..:alus,  and  for  which  the 
length  is  the  n.ean  length  of  the  segment  of  compressor  annulus.  The 
evolution  of  conditions  at  the  downstream  boundary  from  those  at  the 
upstream  boundary  is  made  by  stepping  downstream  through  a number  of 
equal  increments.  The  method  by  m.eans  of  which  the  number  of 
increments  is  chosen  is  described  below.  Each  increment  may  be  con- 
sidered to  have  an  upstream  boundary  and  a downstream  boundary. 
Conditions  at  the  first  upstream  boundary  are  those  at  the  upstream 
boundary  of  the  whole  mixing  interval.  Conditions  at  the  downstream 
boundary  of  an  intermediate  step  becon.e  the  upstream  boundary  conditions 
for  the  next  step,  until  finally  for  the  last  step,  the  downstream  boundary 
conditions  are  those  sought  for  the  downstream  boundary  of  the  mixing 
interval.  Inputs  to  the  calculation  are,  at  the  upstream  boundary,  the  mid- 
streamline static  pressure,  the  mid-streamline  meridional  velocity  (both 
as  computed  previously  by  the  inviscid  equations),  the  radial  distributions 
of  entropy,  total  enthalpy,  and  angular  momentum,  and  the  eddy  viscosity. 
For  each  mixing  step,  the  calculation  is  made  as  follows: 

1)  The  radial  component  equation  (Equation  78)  is  integrated 
to  give  the  static  pressure  distribution  at  the  upstream 
boundary.  Equation  78  is  written  in  the  following  form: 


£ 

I 

i f 
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- it1' 


(126) 


This  equation  is  used  to  step  from  streamline  to  streamline, 
mean  values  of  the  quantities  involved  being  used  in  each 
radial  interval.  As  the  density  distribution  is  not  known, 
this  must  be  determined  itervatively  from  the  static  pressures 
and  the  entropy. 


2)  For  each  streamline,  the  static  pressure  and  entropy  give  the 
static  enthalpy.  The  total  and  static  enthalpies  give  the  total 
velocity.  The  total  velocity  and  tangential  velocity  component 
give  the  axial  velocity  component. 

3)  For  the  first  step  only,  a continuity  summation  is  made  to 
determine  the  flow  rate. 


4'  The  energy  equation  (Equation  82)  is  appliec  on  each  streamline 
to  give  the  total  enthalpy  distribution  at  the  downstream  boundary. 
Equation  82  is  written  as  follows: 
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The  method  of  obtaining  all  radial  derivatives  is  described 
below. 

The  tangential  component  equation  (Equation  79)  is  used  on 
each  streamline  to  give  the  downstream  tangential  velocity 
component  distribution.  Equation  79  is  written: 

k 

Vi  Y\V<r  v / <i28) 


The  mid-streamline,  downstream  static  pressure  is  estimated 
(as  being  equal  to  the  upstream  value). 

The  downstream  radial  distribution  of  static  pressure  is 
determined  as  was  the  upstream  distribution  in  Step  1. 

The  axial  component  equation  (Equation  SO)  is  applied  on  each 
streamline  for  the  downstream  axial  velocity  distribution. 
Equation  80  is  written: 


For  each  streamline,  total  enthalpy,  tangential  and  axial 
velocities  yield  static  enthalpy.  Static  enthalpy  and  static 
pressure  give  density.  A summation  is  made  for  the  flow 
rate  at  the  downstream  boundary. 

If  the  flow  rate  computed  in  Step  9 differs  from  that  established 
in  Step  3,  the  mid-streamline  static  pressure  is  re-estimated, 
and  control  returns  to  Step  7.  (The  method  of  re-estimating 
the  static  pressure  is  given  below.)  Otherwise,  control 
passes  to  Step  11. 


11)  Entropy  at  the  downstream  boundary  is  given  by  static 

enthalpy  and  static  pressure.  Downstream  boundary  conditions 
become  upstream  boundary  conditions  for  the  next  step,  and 
control  returns  to  Step  1. 


Radial  derivatives  of  the  form 


occur  in  Equations  127, 


128,  and  129.  For  streamlines  away  from  the  walls,  the  following  method 
is  used: 


(Vr4WFj)  (vvfj-W 
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Values  near  a wall  must  reflect  the  boundary  condition:  *5^.*  ® 


at 


the  wall.  If  the  streamlines  are  numbered  1,  2,  3, 
then 


starting  at  the  wall, 


(131) 


A value  for  streamline  1 is  obtained  by  linearly  extrapolating  the 
values  for  streamline  3 and  2 to  the  radius  at  the  wall.  The  vaiue  at  streamline 
3 is  given  by  the  general  formula.  Equation  130. 

The  number  of  steps  into  which  each  interval  must  be  divided  is 
controlled  by  the  stability  of  the  three  equations  used  to  step  in  the  axial 
direction.  Consider  the  energy  equation,  and  a linear  variation  of  enthalpy 
with  radius.  The  effect  of  mixing  will  be  to  reduce  the  enthalpy  gradient, 
but  not  to  reverse  the  sign  thereof.  If  an  excessively  large  axial  step  is 
made,  the  sign  of  the  gradient  will  change,  and  so  the  limit  for  stability 
is  when  the  gradient  is  just  eliminated.  Writing  the  energy  equation  in 
approximate  form: 

Alt  Al  r ^ 

^ = y,  032) 
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The  critical  streamlines  are  those  adjacent  to  the  walls.  If  the  stream- 
lines are  numbered  1,  2,  3,  ....  starting  at  the  wall,  and  equispaced,  then, 
from  Equations  130  and  131  (and  for  a linear  enthalpy  distribution): 


a«t*l  dtMci 
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(134) 


where  flih  is  the  change  in  enthalpy  between  stream  lines  and  af  is 
the  streamline  spacing. 

The  limit  for  stability  is  when: 

rr  i)x  + 

^ axial,-  = ti\\ 

and  from  Equations  133  and  134 


(135) 


The  same  analysis  may  be  applied  to  the  tangential  and  axial  component 
equations,  with  the  same  result.  In  practice,  the  limit  of  Az  has  been 
found  to  be  slightly  more  than  half  the  value  given  by  Equation  135.  The 
program  selects  the  number  of  steps  to  be  four  times  that  given  by 
Equation  135,  but  even  this  will  often  lead  to  only  one  step  being  required 
for  each  mixing  interval. 
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The  re -estimation  of  the  n. id-streamline  static  pressure  to  satisfy 
continuity  at  the  downstream  boundary  of  each  mixing  step  (Step  9 in  the 
procedures  above),  is  accomplished  by  determining  the  value  of  change  of 
flow  with  static  pressure  for  the  one -dimensional  flow  of  a perfect  gas. 
Then  the  required  static  pressure  change  is  found  to  be: 


»,  - Wc.u«u)  i> 

P " w 


(136) 


where 


hj  is  the  ratio  of  axial  velocity  to  local  sound  speed. 
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SECTION  VII 


PROGRAM  OPERATION  AND  STRUCTURE 
1.  IMPLEMENTATION 

The  program  is  basically  written  in  USA  Standard  FORTRAN  and 
should  therefore  be  compatible  with  all  current  medium -to -large  computing 
systems,  except  as  noted  below.  All  development  running  has  been  done  on 
a Control  Data  Corporation  6600  computer  at  Wright-Patterson  Air  Force 
Rase.  The  operating  system  was  changed  from  Scope  3.  3 to  Scope  3.4 
during  the  program  development  period.  The  operating  system  includes 
CALCOMP  software  and  hardware  for  on-line  precision  plotting. 

The  program  is  loaded  in  three  principal  overlays  by  a small  resident 
main  program,  and  the  demands  of  the  CDC  Scope  overlay  scheme  are 
incorporated  into  the  coding.  Under  Scope  3.4,  the  central  memory  require- 
ments for  the  three  overlays  are  approximately  as  follows:  the  aerodynamic 
section,  115K;  the  analytic  meanline  blade  section,  125K;  and  the  arbitrary 
meanline  blade  section,  105K.  These  requirements  would  be  reduced  about 
10K  under  the  more  widely  used  Scope  3.  3 system.  The  aerodynamic  section 
is  divided  into  a principal  loading  and  two  alternative  loadings.  The  alterna- 
tive loads  are  Subroutines  UD0308  and  UD0326,  which  handle  two  alternative 
momentum  equation  formulations  selected  by  the  input  variable  NEQN.  In 
order  to  compile  and  run  the  program  under  system s other  than  CDC  Scope, 
the  overlay -related  calls  will  need  to  be  modified.  These  are: 

1)  Card  $MN$.  24  in  Program  UD0300  (the  main  program)  which 
calls  the  analytic  meanline  section,  cards  $AN$.  2 and  $AN$.  3. 

2)  Card  $MN$.  27  in  the  main  program.which  calls  the  aerodynamic 
section,  cards  $AR$.  2 and  $AR$.  3. 

3)  Card  $MN$.  31  in  the  main  program,  which  calls  the  arbitrary 
meanline  section,  cards  $AB$.  2 and  $AB$.  3. 

4)  Card  $AR$.  229  in  Program  UD03AR,  which  calls  Program 
UD0308,  cards  $08$.  2 and  $08$.  3. 

5)  Card  $AR$.  231  in  Program  UD03AR,  which  calls  Program 
UD0326,  cards  $26$.  2 and  $26$.  3. 

With  the  exception  of  the  main  program,  each  of  the  programs  referred  to 
above  is,  in  conventional  FORTRAN  terms,  a Subroutine.  The  PROGRAM 
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statement  attached  to  the  main  program.  Card  $MN$.  2,  is  peculiar  to  the 
CDC  Scope  system,  and  defines  the  following  file  usage  by  the  program: 


FORTRAN 


CDC  File  Name 

Variable  Name 

Value 

Description 

INPUT 

LOG1 

1 

Card  reader 

OUTPUT 

LOG2 

2 

Line  printer 

PUNCH 

LOG3 

3 

Card  punch 

FILE2 

LOG5 

5 

Scratch  file 

FILE3 

LOG6 

6 

Scratch  file 

PLOT 

- 

- 

CALCOMP  plot  file 

The  numerical  values  are  assigned  to  the  FORTRAN  variables  in  the  main 
program  starting  with  card  $MN$.  8.  (LOG4  is  redundant. ) LOG1  and 
LOG2  are  also  defined  in  Subroutine  UDG1  on  cards  $GN$.  4 and  $GN$.  5. 

Calls  are  incorporated  into  the  coding  to  subroutines  that  are  assumed 
to  be  available  as  part  of  a CALCOMP  plotting  package.  All  the  plotting  is 
only  performed  optionally  so  that  the  program  may  be  used  on  a system  that 
does  not  include  these  routines  in  its  library.  However,  it  may  be  necessary 
to  add  a dummy  subroutine  with  the  entry-point  names  PLOT, SYMBOL, 
NUMBER,  AXIS,  and  LINE,  as  some  systems  will  not  permit  program 
execution  with  unsatisfied  external  references. 

2.  PROGRAM  LOGIC 

The  complete  program  comprises  of  42  FORTRAN  programs  and 
subprograms.  They  are  distributed  among  the  three  principal  overlays  as 
follows: 

a.  Resident 

Program  UD0300.  The  main  program,  controls  entry  to  each 

section. 

b.  Aerodynamic  Section 

Program  UD03AR.  The  "main  program"  for  the  aerodynamic 
section,  controls  logic  within  the  section. 

Subroutine  UD0301.  Spline -fit  or  linear  interpolation  and  gradient 

determination. 

Subroutine  UD0302.  Reads  input  data. 
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Subroutine  UD0303. 
Subroutine  UD0304. 
Subroutine  UD0329. 

Subroutine  U DO 30  5. 

Subroutine  UD0306. 

Subroutine  UD0330. 

Subroutine  UD0307. 

Subroutine  UDO 309- 
Subroutine  UD0310. 
Subroutine  UD03 1 1 . 
Subroutine  UD0312. 
Subroutine  UDG1- 
Function  UDG2. 
Function  UDG3. 
Function  UDG4. 
Function  UDG5. 
Function  UDG6. 
Function  UDG7. 
Function  UDG8. 
Function  UDG9. 


Checks  lines  printed  per  page. 

Axisymmetric  turbulent  mixing  calculations. 

Determines  rr-i^al  derivatives  for  mixing 
calculation. 

Handles  interpolation  of  input  data  at  a particular 
computing  station. 

Modification  of  loss  coefficients  to  account  for 
radial  blade-wake  transfers. 

Solves  simultaneous  linear  equations  for  least- 
square  fit. 

Determines  entropy,  rothalpy  or  enthalpy,  and 
whirl  velocity  at  a computing  station. 

Re-estimates  loss  coefficients. 

Annulus  wall  boundary  layer  calculation. 

Regular  printed  output. 

Precision-plot  output. 

Reads  data  describing  fluid. 

Enthalpy  = f (entropy,  pressure) 

Entropy  = f (pressure,  enthalpy) 

Pressure  = f (enthalpy,  entropy) 

Specific  weight  = f (enthalpy,  entropy) 

Enthalpy  = f (pressure,  temperature) 

Temperature  = f (enthalpy,  entropy) 

Ratio  of  specific  heats  - f (enthalpy,  entropy) 

Square  of  Mach  number  - f (static  enthalpy, 
entropy,  square  of  velocity) 
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Subroutine  UDC  125. 


Creates  line-printer  plots. 


■ 

> 


Subroutine  UD0308. 
Subroutine  UD0326 


Determines  simultaneous  solution  to  momentum 
and  continuity  equations  when  NEQN  r 0 or  1. 

Determines  simultaneous  solution  to  momentum 
and  continuity  equations  when  NEQN  = 2 or  3. 


c.  Analytic  Meanline  Blade  Section 


Program  UD03AN.  The  "main  program"  for  the  analytic  mean 
line  blade  section.  Corresponds  to  Program  BLADE,  Ref.  3. 

Subroutine  UD0313.  Corresponds  to  Subroutine  BQ,  Ref.  3. 


Subroutine  UD0314.  Corresponds  to  Subroutine  CQ,  Ref.  3. 


Subroutine  UD0315.  Corresponds  to  Subroutine  Dl,  Ref.  3. 

Subroutine  UD0316.  Corresponds  to  Subroutine  EQ,  Ref.  3. 

Subroutine  UD0317.  Corresponds  to  Subroutine  FQ,  Ref.  3. 

Subroutine  UD0318.  Corresponds  to  Subroutine  GQ,  Ref.  3. 

Subroutine  UD0319.  Determines  deviation  angles  and  formats  data  for 

aerodynamic  section. 

d.  Arbitrary  Meanline  Blade  Section 

Program  UD03AB,  The  "main  program"  for  the  arbitrary  mean- 
line blade  section.  Corresponds  to  Program  ARBITR,  Ref.  4. 

Subroutine  UD0320.  Corresponds  to  Subroutine  BQ,  Ref.  4. 

Subroutine  UDQ321.  Interfaces  with  Subroutine  UD0327. 

Subroutine  UD0322.  Interfaces  with  Subroutine  UD0  327. 

Subroutine  UD0323.  Corresponds  to  Subroutine  EQ,  Ref.  4. 

Subroutine  UD0324.  Corresponds  to  Subroutine  FQ,  Ret.  4. 

Subroutine  UD0327.  Spline-fit  rout.ne  for  interpolation,  gradients  and 

integrals,  also  linear  interpolation. 
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The  overall  program  logic  is  illustrated  by  the  flow  chart  in 
Figure  2.  This  shows  all  the  significant  actions  in  Program  UD0300, 
expressed  by  means  of  the  FORTRAN  variables  and  statement  labels. 
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DO  120  K=  l,  NANAL 

CALL  Analytic  n.eanline  section 

L 

If  (NAERO.  EQ.O) 

* 

CALL  Aerodynamic  section 

1 

If  (NARBIT.  EQ.O) 

i 

DO  150  K=  1.  NARBIT 

CALL  Arbitrary  n.eanline  section 

! 

r 

i 

t 

Stop 


Figure  2.  Overall  Program  Logic 


Details  of  the  operation  of  the  analytic  meanline  section  may 
be  found  in  References  3 and  5.  The  only  significant  modification  made  in 
incorporating  the  program  described  therein  into  Program  UD0300  is  the 
addition  of  Subroutine  UD0319.  This  subroutine  is  entered  once  as  a final 
operation  if  data  from  the  analytic  meanline  blade  section  are  required  as 
input  for  the  aerodynamic  section. 

The  arbitrary  meanline  blade  section  operates  basically  as 
described  in  References  4 and  5,  but  a new  method  of  generating  the  mean- 
line has  been  developed.  Previously,  a number  of  meanlines  were 
examined,  based  upon  a range  of  boundary  conditions  at  the  leading  edge. 
The  n.eanline  was  then  selected  from  those  obtained  on  the  basis  of  the 
minimum  number  of  inflection  points  and  the  largest  minimum  radius  of 
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curvature  on  the  meanline.  The  segment  of  meanline  between  each  point 
where  the  blade  angle  was  specified  (a  computing  station  in  the  aerodynamic 
analysis)  was  an  algebraic  cubic,  so  that  the  resulting  meanline  was  a spline- 
curve  restrained  to  give  the  desired  blade  angle  at  each  computing  station. 

At  the  leading  edge,  the  radius  of  curvature  of  the  meanline  was  defined. 

The  new  procedure  is  to  fit  a spline-curve  through  the  specified  angles  at 
each  computing  station.  Then  the  analytic  integral  of  the  spline  curve  gives 
the  meanline  which  is  therefore  a piece-wise  quartic.  Boundary  conditions 
used  for  the  spline-fit  are  that  the  second  derivative  be  zero  at  each  end. 

The  new  procedure  has  been  found  to  give  meanlines  with  generally  larger 
minimum  radii  of  curvature  than  the  old  procedure,  and,  in  some  instances, 
to  give  a satisfactory  meanline  when  none  could  be  found  using  the  old 
procedure. 


The  logic  of  the  aerodynamic  section  is  illustrated  by  the  flow 
chart  shown  in  Figure  3 , which  is  essentially  a flow  chart  for  Program 
UD03AR.  All  significant  actions  and  calls  to  subroutines  (other  than  to 
Subroutine  UD0301  for  interpolation)  in  the  program  are  shown. 

3.  PROGRAM  AND  SUBPROGRAM  DESCRIPTIONS 

a.  Program  UD0300.  Cards  identified  by  $MN$.  This  is  the  main 
program  and  controls  entry  to  each  section,  as  shown  in  Figure  2 . 

b.  Program  UD03AR.  Cards  identified  by  $AR$.  This  is  the  main 
routine  of  the  aerodynamic  section,  and  it  controls  the  logic  of  the  section, 
as  shown  in  Figure  3.  In  terms  of  conventional  FORTRAN  it  is  a sub- 
routine, but  is  shown  as  a program  because  of  the  use  of  CDC  Scope  ovt  y. 
Some  significant  points  are  distributed  in  the  program  as  follows. 

Operation  First  Related  Card 


Call  UD0302  to  read  input. 

$AR$.  43 

Determine  if  streamline  estimate  is  required. 

$AR$.  51 

Begin  streamline  estimation. 

$AR$,  54 

Determine  geometric  constant  in  relaxation  factor. 

$AR$.  152 

Interpolate  inlet  conditions 

$AR$.  169 

Estimate  inlet  station  velocity. 

$AR$.  191 

Set  integers  to  begin  a pass. 

$AR$.  196 
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Figure  3.  Logic  of  Aerodynamic  Section. 
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Optional  call  to  UD0304  for  turbulent  mixing 
calculations. 

$AR$.  213 

Call  UD0305  to  interpolate  data  at  a bladed  station. 

$AR$.  215 

Optional  call  to  UD0306  for  blade  wake  transfer 
calculations. 

$AR$.  216 

Call  UDQ307  to  determine  entropy,  rothalpy,  or 
enthalpy  and  whirl  velocity. 

$AR$.  218 

Set  entropy,  enthalpy,  and  angular  momentum 
following  blade  free  space. 

$AR$.  220 

Call  UD0308  to  solve  momentum  and  continuity 
equations,  if  NEQN  = 0 or  1. 

$AR$.  229 

Call  UD0^26  to  solve  momentum  and  continuity 
equations,  if  NEQN  = 2 or  3. 

$AR$.  231 

Optional  call  to  UD0309  to  re-estimate  loss 
coefficients. 

$AR$.  234 

Check  convergence  of  meridional  velocity  in  loss 
re -estimation  loop. 

$AR$.  236 

Optional  call  to  UD0310  to  compute  annulus  wall 
boundary  layers. 

$AR$.  258 

Check  on  overall  convergence. 

$AR$.  263 

Re-estimate  streamline  locations. 

$AR$.  267 

Call  UD0311  for  printed  output. 

$AR$.  304 

Call  UD0309  for  printout  of  re-estimated  loss 
coefficients  if  specified. 

$AR$.  309 

Optional  call  to  UD0312  for  static  pressure  and 
streamline  plots. 

$AR$.  311 

Punch  streamline  pattern  if  specified. 

$AR$.  314 

Punch  strean  line  radii  for  analytic  meanline 
section  if  specified. 

$AR$.  318 
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Output  data  for  arbitrary  meanline  section 
if  specified. 

Optional  call  to  UD0325  to  plot  streamline 
relocation  shifts. 

c.  Subroutine  UD0301.  Cards  identified  by  $01$.  Performs 
linear  and  splinefit  interpolation  and  gradient  determination. 

d.  Subroutine  UD0302.  Cards  identified  by  $02$.  Reads  and 
prints  (aerodynamic  section)  input  data. 

e.  Subroutine  UD0303.  Cards  identified  by  $03$.  Checks  and 
updates  number  of  lines  printed  per  page. 

f.  Subroutine  UD0304.  Cards  identified  by  $04$.  Turbulent 

mixing  calculations. 

g.  Subroutine  UD0329,  Cards  identified  by  $29$.  Determines 

radial  derivatives  for  mixing  calculations. 

h.  Subroutine  UD0305.  Cards  identified  by  $05$.  This  subroutine 

handles  the  interpolation  of  data  prescribing  conditions  at  a computing 
station.  If  relative  flow  angle  is  not  specified,  this  is  derived  from  the 
data. 


$AR$.  323 
$AR$.  348 


i.  Subroutine  UD0306.  Cards  identified  by  $06$.  Loss  coefficients 
are  modified  to  account  for  radial  transfer  of  blade  wakes. 

j.  Subroutine  UD0330.  Cards  identified  by  $30$.  Solves 
simultaneous  linear  equations. 

k.  Subroutine  UD0307.  Cards  identified  by  $07$.  If  absolute 
angular  momentum  is  specified  (or  implied)  at  a computing  station,  this 
subroutine  determines  total  enthalpy,  entropy,  and  whirl  velocity.  If 
relative  flow  angle  is  specified,  this  was  determined  in  Subroutine  DD0305, 
and  here  the  rothalpy  and  entropy  are  determined.  Some  significant  points 
in  the  program  are  distributed  as  follows. 

Operation  First  Related  Card 

Determine  enthalpy  and  entropy  given  total  $07$.  33 

pressure  and  loss  coefficient  based  on  local 
dynamic  head. 
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Determine  enthalpy  and  entropy  given  total 
pressure  and  loss  coefficient  not  based  on 
local  dynamic  head. 

Determine  enthalpy  and  entropy  given  total 
pressure  and  isentropic  efficiency. 

Determine  enthalpy  and  entropy  given  total 
pressure  and  entropy  rise. 

Determine  whirl  velocity  from  enthalpy. 

Determine  whirl  velocity  when  enthalpy 
is  given. 

Determine  whirl  velocity  when  angular 
momentum  is  given. 

Set  whirl  velocity  when  given  directly. 

Determine  enthalpy  from  whirl  velocity. 

Determine  entropy  given  loss  coefficient  based 
on  local  dynamic  head  and  enthalpy,  directly  or 
derived. 

Determine  entropy  given  loss  coefficient  not  based 
on  local  dynamic  head,  and  enthalpy,  directly  or 
derived. 

Determine  entropy  given  isentropic  efficiency 
and  enthalpy,  directly  or  derived. 

Determine  entropy  given  entropy  rise  and 
enthalpy,  directly  or  derived. 

Determine  rothalpy,  given  relative  flow  angle. 

Determine  entropy  given  loss  coefficient  based 
on  local  dynamic  head  and  relative  flow  angle. 

Determine  entropy  given  loss  coefficient  not  based 
on  local  dynamic  head  and  relative  flow  angle. 


$07$. 46 

$07$. 62 

$07$. 66 

$07$.  69 
$07$. 72 

$07$.  76 

$07$. 79 
$07$.  81 
$07$. 85 

$07$.  97 

$07$.  112 

$07$. 115 

$07$.  118 
$07$.  122 

$07$.  134 
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Determine  entropy  given  isentropic  efficiency 
and  relative  Dow  angle. 

Determine  entropy  given  entropy  rise  and 
relative  Dow  angle. 


$07$. 146 
$07$. 158 


1.  Subroutine  UD0309.  Cards  identified  by  $09$.  This  subroutine 
'-estimates  the  relative  total  pressure  loss  coefficients.  If  IPRINT  = 0, 
“ new  values  are  placed  in  array  DATA2;  if  IPRINT  = 1,  the  new  values 
^ . . printed  out.  Some  significant  points  in  the  subroutine  are  distributed 
as  follows. 


Operation 

A maximum  value  for  the  calculated 
loss  coefficient  of  0.  7 is  set.  (This  is 
subsequenDy  only  imposed  when  the  value 
is  incorporated  into  the  calculation,  not 
when  it  is  printed  out. ) 

Heading  depending  upon  use  of  calculated 
loss  coefficients  is  printed. 

Cascade  solidity  is  interpolated. 

DiDusion  factor  is  determined. 

The  loss  parameters  are  interpolated  from  the 
10%,  50%,  and  90%  of  blade -height  curves  as 
a function  of  diffusion  factor. 

The  loss  parameters  are  interpolated  as  a 
function  of  radial  location. 

The  diffusion  con. pone nt  of  the  loss  coefficient 
is  determined. 

The  shock  component  of  the  loss  coefficient  is 
omitted  if  NDEL#  = 0. 

The  Prandtl -Meyer  expansion  angle  in  the  input 
data  is  interpolated. 

If  the  blade  geometry  at  the  leading  edge  is  not 
given,  the  incidence  angle  is  not  determined. 


First  Related  Card 


$09$. 28 


$09$.  35 

$09$. 61 
$09$.  77 
$09$.  82 

$09$. 93 
$09$.  101 
$09$.  103 
$09$.  105 
$09$. 106 
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The  expansion  angle  is  the  input  value  plus  the 
incidence  angle  (or  zero). 

The  inlet  relative  Mach  number  is  determined. 


$09$.  132 
$09$.  139 


The  Prandtl- Meyer  angle  following  the  expansion  $09$.  142 

is  determined. 

If  the  Prandtl -Meyer  angle  is  less  than  zero,  $09$.  144 

indicating  compression  to  a Mach  number  of 

less  than  unity,  the  suction  surface  Mach 

number  at  shockwave  impingement  is  set  to 

unity. 

The  Mach  number  at  the  shockwave  impingement  $09$.  149 

is  determined  iteratively. 
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The  mean  Mach  number  at  the  shock  wave  is  the 
average  of  the  relative  inlet  Mach  number  (or 
unity,  if  the  inlet  value  is  less  than  unity)  and 
the  suction  surface  impingement -point  Mach  number. 

$09$. 163 

The  mean  Mach  number  is  multiplied  by  the 
relative  inlet  Mach  number  if  this  was  less 
than  unity. 

$09$. 164 

No  shock-loss  is  computed  if  the  mean  Mach 
number  is  less  than  unity. 

$09$. 165 

The  shock-loss  component  of  the  loss  coefficient 
is  determined. 

$09$. 167 

The  total  loss  coefficient  is  determined. 

$09$. 170 

If  the  calculated  results  are  to  be  printed,  they 
are  not  incorporated  into  the  overall  computation. 

$09$.  171 

The  limiting  value  of  loss  coefficient  is  imposed, 
if  appropriate. 

$09$. 188 

m.  Subroutine  UD0310.  Cards  identified  by  $10$.  This  subroutine 

computes  the  blockage  due  to  the  annulus  wall  boundary  layers. 

n.  Subroutine  UD0311.  Cards  identified  by  $11$.  This  subroutine 
prints  the  regular  output. 
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o.  Subroutine  UD0312.  Cards  identified  by  $12$.  This  subroutine 
makes  plots  of  static  pressure  and  file  streamline  pattern.  All  calls  to 
CAJLCOMP  software  from  the  aerodynamic  section  occur  in  this  subroutine. 

p.  Subroutine  UDG1.  Cards  for  UDG1  and  the  associated  functions 
identified  by  $GN$.  The  constants  describing  the  fluid  are  read  in  by 
Subroutine  UDG1,  and  communicated  to  the  eight  functions  for  thermo* 
dynamic  properties  by  means  of  Common  storage  area  GAS.  This  sub- 
routine and  the  eight  associated  functions  may  be  modified  to  change  the  type 
of  fluid  from  a perfect  gas. 

q.  Functions  UDG2  through  UDG9.  These  eight  functions  give 
thermodynamic  properties  as  listed  in  the  previous  subsection. 

r.  Subroutine  UD0325.  Cards  identified  by  $25$.  This  subroutine 
makes  line-printer  plots  of  the  change  in  £ -coordinates  of  the  mid -stream- 
line  from  pass  to  pass,  for  each  station,  if  requested. 

s.  Subroutine  UD0308.  Cards  identified  by  $08$.  This  is  the  most 
important  subroutine  in  the  aero  section.  It  solves  the  momentum  and 
continuity  equations  when  NEQN  = 0 or  1.  Some  significant  points  in  file 
subroutine  are  distributed  as  follows. 


Operation 

First  Related  Card 

Determine  streamline  slope  angle. 

$08$.  64 

Determine  curvature  of  streamlines. 

$08$.  65 

Determine  streamwise  gradient  of  blade 
blockage. 

$08$. 68 

Determine  streamwise  gradient  of  entropy. 

$08$. 70 

Determine  streamwise  gradient  of  angular 
momentum. 

$08$. 71 

Determine  streamline  slope  angle. 

$08$. 104 
$08$. 107 

Determine  tangent  of  station  lean  angle. 

$08$. 131 

Determine  gradient  along  computing  station  of 
entropy  and  streamline  slope  angle. 

$08$.  140 
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Determine  gradient  along  computing  station 
of  angular  momentmn  and  enthalpy. 

$08$. 143 

Determine  local  sound  speed  for  design  case. 

$08$. 148 

Determine  gradient  along  computing  station 
of  rlsi^  and  rothalpy. 

$08$.  161 

Determine  local  sound  speed  for  analysis  case. 

$08$.  165 

Select  analysis  case  momentum  equation 
solution. 

$08$.  184 

Start  design  case  momentum  equation  solution. 
Determine  static  temperature. 

$08$.  185 

Determine  parameter  h in  momentum  equation. 

$08$. 220 

Determine  parameter  6 in  momentum  equation. 

$08$. 224 

Integrate  momentum  equation. 

$08$. 241 

Impose  a maximum  limit  on  any  meridional 
velocity  of  2. 0 times  the  mid-streamline  value. 

$08$. 271 

Impose  a minimum  limit  on  any  meridional 
velocity  of  1.0. 

$08$.  280 

New  estimate  of  velocity  is  stored  temporarily. 

$08$. 292 

Meridional  velocity  is  restra^ed  to  be  not 
greater  than  sound  speed. 

$08$.  293 

Relaxation  factor  is  determined  and  applied 
to  give  new  velocity  estimate. 

$08$.  303 

Call  UD0307  if  entropy  is  a function  of  velocity. 

$08$. 312 

Determine  static  enthalpy  and  square  of 
meridional  Mach  number. 

$08$. 313 

Start  analytic  case  momentum  equation  solution. 

$08$.  330 

Determine  static  temperature. 

$08$.  342. 
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Determine  meridional  velocity  gradient. 


$08$.  381 


Determine  meridional  velocity  on  adjacent 
streamline. 

Impose  a maximum  limit  on  the  meridional 
velocity  of  3.0  times  the  mid- streamline 
value. 

Impose  a minimum  limit  on  the  meridional 
velocity  of  1.0. 

Impose  a maximum  limit  on  the  meridional 
velocity  of  the  local  sound  speed. 

Check  convergence  of  meridional  velocity. 

If  entropy  is  a function  of  velocity,  call  UD0307. 

Determine  static  enthalpy  and  square  of 
relative  Mach  number. 

Start  continuity  equation  solution.  Integrate 
continuity  equation  and  equation  for  gradient 
of  flow  with  velocity. 

Normalize  fractions  of  flow  in  each  streamline 
for  streamline  re -estimation. 

Check  for  supersonic  solution. 

Determine  meridional  velocity  increment. 

Impose  limits  on  velocity  increment. 

Add  increment  to  meridional  velocities, 
impose  minimum  value  of  1.0,  maximum 
value  of  sound  speed. 

Check  convergence  of  momentum  and  continuity 
equations. 

Call  UD0307  if  entropy  is  a function  of  velocity. 
Determine  relative  flow  angle  for  design  cases. 


$08$. 397 
$08$.  398 

$08$. 410 

$08$. 424 

$08$. 429 
$08$.  452 
$08$. 454 

$08$. 478 

$08$. 506 

$08$. 508 
$08$. 511 
$08$. 512 
$08$. 519 

$08$. 555 

$08$. 565 
$08$. 592 
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$08$. 597 


Determine  enthalpy  and  whirl  velocity  for  analysis 
cases. 


t.  Subroutine  UD0326.  Cards  identified  by  $26$.  This  subroutine 
is  very  similar  to  Subroutine  UD0308,  and  it  solves  the  momentum  and 
continuity  equations  when  NEQN  = 1 or  2.  It  was  created  by  modification 
of  Subroutine  UD0308,  and  so  is  structured  identically.  Subroutine  UD0326 
does  not  have  the  local  sound  speed  as  an  upper  limit  on  meridional 
velocity,  and  the  A and  B parameters  of  the  momentum  equation  reflect 
the  alternative  formulations. 


u.  Program  UD03AN.  Cards  identified  by  $AN$.  This  is  the  main 
routine  of  the  analytic  meanline  blade  section  of  the  program.  In  terms  of 
conventional  FORTRAN  it  is  a subroutine,  but  appears  as  a program  here 
because  of  the  use  of  CDC  Scope  overlay.  See  Program  BLADE,  Reference 
3 for  details. 


v.  Subroutine  UD0313.  Cards  identified  by  $13$.  See  Subroutine 
BQ,  Reference  3 for  details. 

w.  Subroutine  UD0314.  Cards  identified  by  $14$.  See  Subroutine 
CO,  Reference  3 for  details. 

x.  Subroutine  UD0315.  Cards  identified  by  $15$.  See  Subroutine 
Dl,  Reference  3 for  details. 

y.  Subroutine  UD0316.  Cards  identified  by  $16$.  See  Subroutine 
EQ,  Reference  3 for  details. 

z.  Subroutine  UD0317.  Cards  identified  by  $17$.  See  Subroutine 
PQ,  Reference  3 for  details. 


aa.  Subroutine  UD0318.  Cards  identified  by  $18$.  See  Subroutine 
GQ,  Reference  3 for  details. 


ab.  Subroutine  UD0319.  Cards  identified  by  $19$.  This  subroutine 
determines  deviation  angles  and  hence  relative  flow  angles,  reads  in  other 
data  regarding  the  blade,  and  puts  onto  a scratch  file,  (LOG5),  data  for 
the  aerodynamic  section. 


ac.  Program  UD03AB.  Cards  identified  by  $AB$.  This  is  the  main 
routine  of  the  arbitrary  meanline  blade  section  of  the  program.  In  terms 
of  conventional  FORTRAN  it  is  a subroutine,  but  appears  as  a program  here 
because  of  the  use  of  CDC  Scope  overlay.  See  Program  ARB1TR,  Reference 
4 for  details,  except  that  some  simplification  occurs  through  use  of  the  new 
meanline  generation  procedure. 


Ill 


a 


ad.  Subroutine  UDO320.  Cards  identified  by  $20$.  See  Subroutine 

BQ,  Reference  4 for  details. 

ae.  Subroutine  UD0321.  Cards  identified  by  $21$.  Calls  Subroutine 
UD0327. 

af.  Subroutine  UD0322.  Cards  identified  by  $22$.  Calls  Subroutine 

UD0327. 

ag.  Subroutine  UD0323.  Cards  identified  by  $23$.  See  Subroutine 
EQ.  Reference  4 for  details. 

ah.  Subroutine  UD0324.  Cards  identified  by  $24$.  See  Subroutine 
IQ,  Reference  4 for  details. 

al.  Subroutine  UD0327.  Cards  identified  by  $27$.  This  subroutine 
provides  interpolation,  derivatives  and  integrals  based  upon  spline  fitting, 
and  also  linear  interpolation. 
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LIST  OF  SYMBOLS 


A 

a/c 

B 

B 

\> 

c 

b 

E 

F 

F. 

F 

3 

H 

H 

k 

I 

J 

Kfe 

JL 

n 

Km 

K. 

JV 

n> 

P 

P' 

Pr 


Parameter  in  momentum  equation 

Ratio  of  distance  to  maximum  camber  point  to  blade  chord 
Parameter  in  momentum  equation 

Blockage  due  to  annulus  wall  boundary  layers  and  blade  wakes 
Exponent  in  deviation  rule,  from  Figure  164,  Reference  6 
Blade  chord 
Diffusion  factor 

Arbitrary  error  in  stxeamline  location 
Optimum  streamline  relocation  relaxation  factor 
Radial  force  on  element  of  boundary  layer 
Radial  movement  of  boundary  layer  material 
Acceleration  due  to  gravity 
Total  enthalpy 

Boundary  layer  shape  factor,  ratio  of  displacement  to  momentum 
thicknesses 

Static  enthalpy 

Rothalpy,  defined  by  Equation  48 

Joules  equivalent 

Shape  factor  in  deviation  rule 

Thickness  factor  in  deviation  rule,  from  Figure  172,  Reference  6 
Distance  along  computing  station  outwards  from  hub 
Mach  number 

Ratio  of  meridional  velocity  component  to  speed  of  sound 
Ratio  of  tangential  velocity  component  to  speed  of  sound 
Ratio  of  relative  velocity  to  speed  of  sound 
Meridional  distance  in  direction  of  flow 
Total  pressure 

Ideal  (isentropic)  total  pressure 
Relative  total  pressure 
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LIST  OF  SYMBOLS  (Continued) 


|>  Static  pressure 

Q Rate  of  change  of  flow  with  meridional  velocity  at  mid-radius 

R Streamline  relocation  relaxation  factor 

r Radius 

Radius  of  curvature  of  stream  surface  projected  onto  meridional 
plane 

S Entropy 

X Total  temperature 

£ Static  temperature 

U « Free-stream  velocity 

y Velocity 

Radial  component  of  velocity 
Axial  component  of  velocity 
V*  Meridional  component  of  velocity 

V#  Tangential  or  whirl  component  of  velocity,  positive  in  direction 

of  rotor  rotation 

W Flow  rate 

w Relative  velocity 

Wo  Relative  tangential  velocity  component 

l n ) Specific  weight 

(t>|.  Specific  weight  based  on  total  conditions 

X Axial  coordinate 

<?(  Absolute  whirl  angle,  defined  by  Equation  1 

^ Relative  flow  angle,  defined  by 

Ratio  of  specific  heats 

$ Station  lean  angle,  defined  by  ta,r  = Vi>e 

£ Deviation  angle 

£ Reference  deviation  from  Figure  161,  Reference  6 

0,0  K 

^ Blade  lean  angle  on  a computing  station,  defined  by  1aw$ Pff 

with  8 increasing  in  the  direction  of  rotor  rotation  T bx 
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LIST  OF  SYMBOLS  (Continued) 


€ 

a 

6 

© 

<9 

X 

\? 

C 

eo 


Blade  lean  angle  projected  onto  constant-x  plane,  defined  by 
Equation  61 

Eddy  viscosity 

Isentropic  efficiency 

Camber  angJ 

Circumferential  coordinate 

Boundary  layer  momentum  thickness 

Area  coefficient,  defined  (1-  Blade  blockage) 

Kinematic  viscosity 
Density 

Cascade  solidity 

Boundary  layer  shear  stress  ». 

-L  / Or 

Streamline  slope  angle,  defined  by  f — ^ 

Speed  of  rotation  of  blading  and  coordinate  system 
Relative  total  pressure  loss  coefficient 
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